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AN  INVESTIGATION  OF  FLAME  STABILITY 
IN  A  COAXIAL  DUMP  COMBUSTOR 


by 

Edward  T,  Curran,  Ph.D. 

Dr.  Harold  E.  Wright,  Advisor 


An  experimental  investigation  of  the  flame  stability  characteristics 
of  a  coaxial  dump  combustion  chamber  was  made.  A  s^x-inon  diameter  com¬ 
bustor  was  studied,  and  the  dump  plane  inlet  diameter  was  varied  from 
2.5>h.to  5.0Wc^.Tests  were  carried  out  using  pre-mi xed  ethylene/air  and 
JP-4/air,  To  avoid  vitiation  effects,  a  clean-air  heater  was  used  to 
vary  the  temperature  of  the  inlet  air.  The  temperature  levels  were 
generally  in  the  range  1 000-1 250^R;  the  corresponding  inlet  pressure 
was  varied  between  approximately  10  and  40  psia.  Tnlet  Mach  numbers 
were  typically  in  the  range  0.2  to  0.9, 

The  flame  stability  data  were  not  collapsed  by  a  conventional  dump 
combustor  correlation.  Howeve  •,  a  correlating  parameter  was  derived  by 
modelling  the  annular  recirculation  zone  as  an  adiabatic  stirred  reactor 
using  a  one-step  chemical  reaction.  This  derived  correlating  parameter 
was  used  to  obtain  successful  correlations.  The  dominant  variable  in 
flame  stability  was  the  inlet  temperature;  the  subsidiary  variables  were 
the  step  height,  inlet  pressure,  and  velocity.  The  effect  of  velocity 
on  the  lean  extinction  limit  was  very  smaT'K  With  OP-4  fuel  the  exist¬ 
ence  of  complete  stabilityjopps  was  demonstrated  by  operating  at  very 
low  combustor  pressures. 

Some  tests  were  performed  with  additional  fuel  Injection  directly 
into  the  recirculation  zone.  Lean  extinction  limits  were  determined 
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with  varying  fractions  of  the  total  fuel  flow  being  injected  into  the 
recirculation  zone.  These  tests  yielded  an  estimate  of  the  air  toss 
flow  fraction  entrained  into  the  recirculation  zone. 

Finally,  tests  were  performed  using  a  transparent  quartz  combustor; 
it  was  observed  that  the  combustion  process  was  an  oscillatory  phenomenon 
with  the  reaction  zone  moving  rapidly  to  and  fro  along  the  combustor. 

A  comprehensive  review  of  the  literature  concerning  the  flame  stability 
and  aerodynamics  of  coaxial  dump  combustors  was  presented. 
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AN  INVESTIGATION  OF  FLAME  STABILITY 
IN  A  COAXIAL  DUMP  COMBUSTOR 


by 

Edward  T.  Curran,  Ph.D. 

Dr.  Harold  E.  Wright,  Advisor 

An  experimental  investigation  of  the  flame  stability  characteristics 
of  a  coaxial  dump  combustion  chamber  was  made.  A  six-inch  diameter  com¬ 
bustor  was  studied,  and  the  dump  plane  inlet  diameter  was  varied  from 
2.5"  to  5.0".  Tests  were  carried  out  using  pre-mixed  ethylene/air  and 
JP-4/air.  To  avoid  vitiation  effects,  a  clean-air  heater  was  used  to 
vary  the  temperature  of  the  inlet  air.  The  temperature  levels  were 
generally  in  the  range  1000-1250°R;  the  corresponding  inlet  pressure 
was  varied  between  approximately  10  and  40  psia.  Inlet  Mach  numbers 
were  typically  in  the  range  0.2  to  0,9. 

The  flame  stability  data  were  not  collapsed  by  a  conventional  dump 
combustor  correlation.  However,  a  correlating  parameter  was  derived  by 
modelling  the  annular  recirculation  zone  as  an  adiabatic  stirred  reactor 
using  a  one-step  chemical  reaction.  This  derived  correlating  parameter 
was  used  to  obtain  successful  correlations.  The  dominant  variable  in 
flame  stability  was  the  inlet  temperature;  the  subsidiary  variables  were 
the  step  height,  inlet  pressure,  and  velocity.  The  effect  of  velocity 
on  the  lean  extinction  limit  was  very  small.  With  JP-4  fuel  the  exist¬ 
ence  of  complete  stability  loops  was  demonstrated  by  operating  at  very 
low  combustor  pressures. 

Some  tests  were  performed  with  additional  fuel  injection  directly 
into  the  recirculation  zone.  Lean  extinction  limits  were  determined 


with  varying  fractions  of  the  total  fuel  flow  being  injected  into  the 
recirculation  zone.  These  tests  yielded  an  estimate  of  the  air  mass 
flow  fraction  entrained  into  the  recirculation  zone. 

Finally,  tests  were  performed  using  a  transparent  quartz  combustor; 
it  was  observed  that  the  combustion  process  was  an  oscillatory  phenomenon 
with  the  reaction  zone  moving  rapidly  to  and  fro  along  the  combustor. 

A  comprehensive  review  of  the  literature  concerning  the  flame  stability 
and  aerodynamics  of  coaxial  dump  combustors  was  presented. 
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I  INTRODUCTION 


The  subject  of  this  study  is  the  "sudden-expansion"  combustion 
system  which  has  recently  been  utilized  in  a  new  generation  of  missile 
propulsion  engines,  known  as  integral  rocket-ramjet  (IRR)  engines. 

This  engine  design  was  conceived  as  a  method  of  reducing  missile  vol¬ 
ume  by  packaging  the  required  rocket  booster  within  the  ramjet  engine 
combustor.  The  contrast  between  the  conventional  ramjet  engine  and 
the  IRR  engine  is  shown  in  Fig  1.  The  "sudden-expansion"  (SUE)  combus¬ 
tor  is  the  simplest  version  of  a  class  of  combustors  generically  des¬ 
cribed  as  "dump"  combustors,  examples  of  which  are  illustrated  in 
Fig  2.  It  will  be  noted  that  the  SUE  combustor  is  essentially  of 
coaxial  geometry,  and  it  is  created  simply  by  an  abrupt  increase  in 
cross-sectional  area.  This  abrupt  change  forces  a  separation  of  flow 
and  a  corresponding  annular  recirculation  zone  (RZ)  is  generated.  The 
flow  structure  created  provides  a  region  in  which  an  initial  flame  can 
be  stabilized  and  from  which  flame  can  be  propagated  into  the  entering 
fuel -air  mixture. 

The  introduction  of  "dump"  combustors  into  ramjet  engine  systems 
occurred  quite  rapidly  without  the  undergirding  of  a  large  data  base, 
and  the  initial  development  test  effort  was  largely  concentrated  on 
determination  of  the  overall  factors  which  influence  combustion  effi¬ 
ciency,  pressure  loss,  and  chamber  integrity.  Such  experimental 
testing  was  essentially  parametric  in  nature  and  yielded  urgently 
needed  performance  data  (Ref  1).  Furthermore,  there  were  little  data 
available  from  related  fields  of  research,  although  some  data  of 
relevance  were  subsequently  found  in  industrial  furnace  practice. 
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This  situation  stood  in  marked  contrast  to  that  existing  during 
the  intensive  period  of  development  of  conventional  ramjet  engines 
(circa  1950-1960).  At  that  time,  a  vast  amount  of  very  fruitful 
research  was  stimulated,  which  laid  an  adequate  foundation  for  the 
parallel  development  of  combustion  systems  for  gas  turbines,  after¬ 
burners,  and  ramjets.  A  representative  selection  of  work  reported 
during  this  period  may  be  found  in  the  technical  literature  (Refs  2 
through  8). 

In  view  of  the  paucity  of  data  concerning  dump  combustors  and  the 
general  lack  of  understanding  of  these  devices,  a  number  of  long  term 
investigations  has  been  initiated  at  various  research  organizations, 
and  examples  of  recent  progress  have  been  published  (Refs  9  and  10). 
However,  an  important  and  relatively  neglected  area  of  investigation 
which  is  vital  to  predicting  dump  combustor  performance  is  that  of 
flame  stability. 

The  term  "flame  stability"  is  used  to  describe  the  retention  of 
effective  combustion  within  the  chamber.  A  more  precise  term  is  thus 
"flame-holding,"  and  the  devices  which  serve  to  anchor  the  flame  sys¬ 
tem  within  the  combustor  of  conventional  ramjet  engines  were  oft-g-' 
termed  "flame-holders."  However,  current  usage  is  to  refer  to  such 
devices  as  "flame-stabilizers,"  and  the  process  of  flame  retention  as 
"flame  stabilization."  Generally,  flame  can  be  held  within  a  combustor 
over  a  range  of  fuel /air  ratios  varying  from  a  weak  mixture  limit  to  a 
rich  mixture  limit.  For  a  given  combustor  geometry,  the  range  of 
operation  is  primarily  dependent  on  combustor  entry  conditions  such  as 
the  velocity  (U),  temperature  (T),  and  pressure  (P).  For  a  single 
flame-stabilizing  element,  such  as  a  cone,  the  range  of  stable  operation 


may  be  displayed  in  terms  of  a  parameter  of  the  form  u“/P^D^TS  where 

Df  is  a  representative  dimension  of  the  flame-holder.  A  typical  flame 
stability  loop  for  such  a  bluff  body  is  shown  in  Fig  3,  and  in  this 
case  a  =  0  =  y  =  1,5=  1.5.  Thus,  for  a  given  size  of  flame- 
holder,  the  region  of  stable  operation  is  enhanced  by  low  velocity, 
high  temperature,  and  high  pressure.  For  an  integral  rocket-ramjet 
engine  operating  in  the  flight  regime  from  Mach  2.5  to  Mach  4.0,  the 
typical  range  of  variation  for  these  parameters  is:  velocity,  900  to 
1200  fps;  temperature,  1100  to  1 600°R ;  and  pressure,  10  to  100  psi. 

In  general,  it  will  be  impossible  to  maintain  flame  at  extreme  condi¬ 
tions  of  high  velocity  or  at  too  low  ^  temperature  or  pressure.  In 
practice,  the  weak  limit  is  of  more  importance  than  the  rich  limit 
because,  for  efficient  cruising  flight,  it  is  necessary  to  operate  the 
engine  at  very  lean  fuel-air  ratios.  The  rich  limit  is  usually 
achieved  at  fuel-air  ratios  greater  than  stoichiometric,  a  region  in 
which  it  is  not  normally  desired  to  erate  the  engine.  Additionally, 
as  fuel-air  ratios  approach  the  rich  limit,  experience  shows  that  very 
unstable  operating  conditions  are  encountered  which  result  in  severe 
thermal  and  mechanical  stresses  on  the  engine. 

It  is  thus  apparent  that  the  phenomenon  of  flame  stabilization  is, 
in  effect,  a  switch  which  controls  the  "on-off"  operation  of  the  com¬ 
bustion  system.  As  such,  it  limits  the  operating  regime  of  the  combus¬ 
tor  and  of  the  engine.  Furthermore,  if  extinction  of  combustion  occurs 
with  an  unmanned  vehicle  during  flight,  it  is  usually  impractical  to 
automatically  re-establish  combustion.  Consequently,  the  vehicle  flight 
is  terminated.  Thus,  for  all  classes  of  ramjet  engines,  the  maintenance 
of  flame  stability  is  vital. 
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Figure  3.  Typical  Flame  Stability  Loop  for  a  Bluff  Body 


In  the  case  of  the  IRR  combustor,  the  situation  is  more  severe 
than  for  conventional  ramjet  combustors.  The  latter  may  employ  mech¬ 
anical  flame-stabilizers  (e.g.,  gutters),  but  in  the  case  of  the  dump 
combustor,  no  stabilizers  can  be  inserted  into  the  combustor  because 
of  the  presence  of  the  rocket  grain.  Consequently,  flame  stability 
must  be  achieved  by  utilizing  the  flow  structure  created  by  the  abrupt 
change  in  area;  this  primary  means  of  stabilization  may  be  augmented 
by  low-blockage  mechanical  elements  located  at  the  dump  plane.  How¬ 
ever,  these  latter  elements  are  highly  undesirable  because  of  their 
attendant  pressure  losses.  As  noted  earlier,  little  data  are  available 
for  the  flame  stability  of  the  axi symmetric  step  geometries  associated 
with  coaxial  dump  combustors,  and  no  proven  parameter  exists  for 
correlating  the  data  that  are  available. 

Consequently,  the  first  goal  of  this  exploratory  investigation 
was  to  determine  the  flame  stability  characteristics  of  coaxial  dump 
geometries  of  various  diameter  ratios  (d/D),  over  a  range  of  inlet 
operating  conditions.  A  second  goal  was  to  construct  a  theoretical 
model  of  the  flame  stabilizing  mechanism  and  to  determine  a  suitable 
parameter  that  would  correlate  available  experimental  data. 

The  experimental  investigation  of  flame  stability  limits  was 
carried  out  in  the  test  facilities  of  the  Air  Force  Aero  Propulsion 
Laboratory.  A  baseline  combustor  diameter  of  6.0"  was  chosen  to 
match  the  facility  capability;  the  corresponding  length  was  30.0".  The 
diameter  of  the  inlet  to  the  combustor  was  varied  from  2.5"  to  5.0"  and 
thus  the  corresponding  step  heights  varied  from  1.75"  to  0.5".  Six 
basic  combustor  geometries  were  investigated.  Furthermore,  the  exhaust 
nozzle  area  was  also  varied  to  provide  various  inlet  velocity  levels 


to  the  chamber.  Typically,  inlet  Mach  numbers  were  in  the  range  0.2 
to  0.9.  No  flameholders  were  employed.  To  avoid  the  vitiation  effects 
associated  with  combustion  pre-heaters,  a  clean  air  heater  was  uti¬ 
lized  in  this  investigation  to  vary  the  temperature  of  the  inlet  air. 
The  temperature  levels  were  generally  in  the  range  1  COO-1 250°R;  the 
corresponding  inlet  pressure  was  varied  between  10  and  40  psia.  These 
test  conditions  were  selected  in  order  to  simulate  relatively  severe 
operating  conditions.  All  tests  were  carried  out  using  premixed  fuel 
and  air.  Both  ethylene  and  JP-4  were  used  as  fuels.  In  addition  to 
the  primary  tests  to  investigate  flame  stability,  tests  were  carried 
out  to  investigate  the  effect  of  additional  fuel  injection  into  the 
recirculation  zone;  also  a  transparent  quartz  combustor  was  briefly 
installed  to  permit  visual  study  of  the  internal  flame  pattern. 

This  report  basically  consists  of  four  main  divisions.  Initially, 
a  brief  review  of  the  flame  stability  characteristics  of  bluff  bodies 
is  given:  both  homogeneous  reactor  and  characteristic  time  theories 
of  flame  stabilization  are  reviewed.  This  section  is  concluded  by  a 
survey  of  the  relatively  meager  data  concerning  dump  combustors. 
Secondly,  an  analysis  of  the  aerodynamics  and  flame  stability  of 
coaxial  dump  combustors  is  presented;  this  is  followed  by  a  descrip¬ 
tion  of  the  experimental  program.  Finally,  the  results  of  this 
investigation  are  presented,  and  appropriate  conclusions  are  drawn. 

In  the  next  chapter,  a  brief  review  will  initially  be  given  of 
earlier  work  concerning  bluff-body  flame  stabilization,  followed  by  a 
survey  of  both  the  available  SUE-combustor  results  and  data  from 
similar  geometries  such  as  steps  and  recesses. 
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II  FLAME  STABILITY  AND  COMBUSTOR  OPERATION 


The  primary  function  of  the  combustion  chamber  is  to  maximize  the 
available  energy  released  by  the  combustion  process.  This  requires 
that  the  combustor  operate  at  a  high  combustion  efficiency  (t)c) ,  and 
with  minimum  pressure  loss  (APt)-  for  a  given  level  of  available 
energy  there  will  exist  a  region  of  trade-off  between  efficiency  and 
pressure  loss  which  permits  some  flexibility  in  the  engineering  design 
of  the  combustion  system.  The  combustion  process  must  be  maintained 
over  a  given  range  of  engine  operation.  Also,  the  process  must  usually 
be  accomplished  in  a  small  volume,  i.e.,  at  high  intensity,  and  the 
structural  integrity  of  the  combustor  must  be  assured. 

The  combustion  process  itself  obviously  requires  an  intimate 
contact  of  fuel  and  oxidant  elements  at  a  point  in  their  respective 
flow  histories  where  rapid  reaction  is  possible  and,  for  effective 
operation,  the  corresponding  energy  release  must  occur  before  the 
products  leave  the  combustor.  In  addition  to  the  turbulent  mixing 
and  molecular  diffusion  processes  in  the  combustor  which  bring  about 
intimate  contact  of  the  reactants,  it  should  be  borne  in  mind  that 
larger  scale  convective  flow  patterns  exist  which  transfer  energy  from 
the  combustion  products  to  the  entering  reactants  to  sustain  the 
process.  These  "backmixing"  flows  are  of  particular  importance  in 
dump  combustors.  The  flow  field  in  the  combustor  is  obviously  made 
much  more  complex  by  the  concurrent  and,  in  many  cases,  coupled  pro¬ 
cesses  of  transport  and  reaction.  Also  for  heterogeneous  processes, 
say  combustion  of  a  liquid  fuel,  the  processes  of  fuel  injection, 
atomization,  evaporation,  and  diffusion  can  be  identified,  without 
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exhausting  the  manifold  mechanisms  contributing  to  the  overall 
conversion  of  reactants  into  products.  The  introduction  of  simpli¬ 
fied  concepts  to  describe  the  complex  combustion  processes  occurring 
in  practical  combustors  was  one  of  the  key  contributions  of  the  work 
performed  in  the  period  1950-1960.  One  such  concept  was  that  of  the 
homogeneous  reaction  zone,  later  termed  the  Perfectly  Stirred  Reactor 
(PSR).  The  essential  features  of  this  concept  are  twofold;  namely, 
that  the  reactants  introduced  into  the  reactor  are  instantaneously  mixed 
throughout  the  reactor  volume,  and  secondly,  that  reaction  is  irrmediate. 
It  follows  then  that  the  exit  products  have  the  same  composition  and 
temperature  as  the  gases  within  the  reactor.  Although  this  is  indeed 
an  idealized  concept,  it  was  found  in  practice  that  certain  zones  in 
practical  combustors  were  so  intensively  mixed  (or  well-stirred)  that 
they  approached  the  ideal  PSR  behavior.  Of  course,  it  should  be  noted 
that  the  assumption  of  PSR  operation  permits  relatively  simple  calcula¬ 
tion  of  reactor  behavior.  Initially  the  stirred  reactor  concept  was 
applied  quite  widely  to  the  analysis  of  the  flow  regions  behind  bluff- 
body  stabilizers,  and  this  approach  permitted  a  rational  prediction  of 
flame  stability  limits.  The  concept  also  permitted  estimation  of  the 
upper  limit  of  combustion  intensity  when  mixing  limitations  were 
removed.  It  should  be  noted,  however,  that  stable  operation  of  a 
PSR  is  only  possible  when  the  combustion  efficiency  is  less  than  100%. 
This  latter  comment  serves  to  introduce  another  significant  contribu¬ 
tion  which  was  due  to  Bragg  (Ref  11):  he  proposed  that  for  maximum 
combustion  intensity,  the  ideal  combustion  chamber  should  be  one  in 
which  a  PSR  section  was  followed  immediately  by  a  Plug  Flow  Reactor 
(PFR),  as  illustrated  in  Fig  4.  This  latter  reactor  is  characterized 
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as  one  in  which  all  the  elements  move  through  the  reactor  with  constant 
and  equal  velocity  with  no  longitudinal  mixing  of  the  fluid  along  the 
flow  path.  The  function  of  the  PFR  is  to  permit  the  unburned  reactants 
from  the  PSR  section  to  be  completely  consumed. 

The  utility  of  the  PSR-PFR  concept  is  that  in  many  practical 
combustion  systems,  a  stirred  section  is  essential  to  provide  flame 
stability,  and  such  a  section  can  be  physically  identified  within  the 
combustor,  together  with  a  section  which  enables  the  combustion  reaction 
to  be  effectively  completed.  A  convincing  experimental  demonstration 
of  this  concept  is  given  by  Beer  and  Lee  (Ref  12).  Furthermore,  the 
basic  PSR-PFR  concept  can  be  generalized  so  that  complex  combustors  can 
be  represented  by  an  appropriate  network  of  PSR  and  PFR  elements.  Ulti¬ 
mately,  as  increasing  numbers  of  elemental  reactors  are  used,  one 
approaches  a  finite-difference  solution  of  the  combustor  flow  field. 

The  elaboratory  of  this  reactor-module  approach  is  dealt  with  in  both 
the  text  due  to  Vulis  (Ref  13)  and  that  of  Wen  and  Fan  (Ref  14).  Atten¬ 
tion  is  also  drawn  to  the  contributions  of  Essenhigh  (Refs  15,  16,  and 
17)  and  Swithenbank  (Refs  18  and  19)  to  this  literature. 

At  this  point,  it  is  appropriate  to  note  that  the  assumption  of 
a  perfectly  stirred  reaction  zone  determines  an  age  distribution  func¬ 
tion  for  particles  leaving  the  zone.  It  can  be  shown  that  this 
"residence  time  distribution"  (RTD)  function  is  of  an  exponential 
decay  form  (Ref  20).  Similarly,  for  a  plug  flow  reactor,  where  the 
particles  entering  the  reactor  at  a  given  time  all  subsequently  exit 
together,  the  RTD  function  is  a  simple  pulse,  or  delta  function, 
located  at  the  given  residence  time.  The  RDT  functions  are  also  illus¬ 
trated  in  Fig  4.  The  relevance  of  this  observation  is  that  zones  may 


be  identified  within  a  combustor  which  approximate  PSR  or  PFR  behavior 
by  observing  the  corresponding  residence  time  distributions.  However, 
such  observations  are  usually  carried  in  cold  flow  models,  particularly 
in  water  tunnel  models  of  combustion  systems.  Furthermore,  if  a  com¬ 
bustor  flow  field  is  modelled  by  a  network  of  PSR-PFR  elements,  the 
overall  RTD  of  the  combustor  model  may  be  experimentally  determined  to 
test  the  suitability  of  the  assumed  network.  Unfortunately,  such 
determinations  do  not  necessarily  confirm  the  uniqueness  of  the  model 
configuration. 

In  summary,  the  modelling  of  a  combustor  flow  by  a  reactor  network 
offers  a  physically  realistic  approach  which,  in  varying  degrees,  is 
capable  of  experimental  verification.  The  network  may  be  constructed 
at  various  levels  of  sophistication  while  maintaining  a  structured 
approach  to  the  representation  of  complex  reactor  flows.  In  this 
current  investigation,  the  representation  of  the  flame  stabilizing 
region  by  a  stirred  reactor  is  of  immediate  relevance,  and  this  dis¬ 
cussion  will  be  engaged  in  the  next  section. 

Stirred  Reactor  Characteristics 

In  this  section,  the  stable  operating  range  of  a  stirred  reactor 
is  discussed,  and  the  application  of  this  concept  to  the  prediction  of 
flame  stability  limits  is  introduced.  Following  this  initial  dis¬ 
cussion,  a  review  of  experimental  data  is  given,  and  a  discussion  of 
various  theories  of  flame  stabilization  is  presented. 

The  adiabatic  operation  of  a  stirred  reactor  is  now  considered. 

For  the  homogeneous  combustion  of  fuel  and  oxidizer,  various  reaction 
mechanisms  may  be  postulated.  However,  it  will  become  apparent  later 
in  this  discussion  that  a  simple,  one-step  reaction  model  is  quite 


adequate  to  describe  the  phenomenological  behavior  of  practical  flame 
stabilizing  systems.  For  purposes  of  illustration,  an  overall  reaction 
of  order  two  is  now  considered.  In  particular,  the  reaction  is  assumed 
to  be  of  order  one  for  both  fuel  and  oxidizer. 

The  volumetric  reaction  rate  (r)  may  then  be  written  as 

r  =  k,j\/Tr  e  'E/RT''||c0  Cfj  (1) 

where  CQ,  denote  the  concentrations  of  oxygen  and  fuel,  and  Tr  is 
the  reaction  temperature.  It  is  important  to  note  in  this  equation 
that  in  the  initial  stages  of  the  reaction,  the  reactants  are  consumed 
with  little  rise  in  temperature.  Subsequently,  the  exponential  tempera¬ 
ture  term  soon  gives  rise  to  a  rapid  increase  in  reaction  rate  until 
depletion  of  the  reactants  causes  the  rate  of  decrease.  The  typical 
progression  of  this  type  of  adiabatic  reaction  is  shown  in  Fig  5. 

Using  Eq ( 1 ) ,  the  rate  of  oxygen  consumption  by  the  reactor  may  be 
expressed  as: 


Typically  for  a  weak  mixture  of  hydrocarbon  and  air,  the  stoichiometry 
of  the  combustion  reaction  may  be  written  as  follows: 

<t>  Cx  Hy  +  (x  +  y/4)  02  +  3.76  (x  +  y/4)  N2 

=  x  «  <fi  CO2  +  y/2  <  *  H2O  +  (1-0  <t>  Cx  Hy  +  (x  +  y/4) 

(1  -  <0  02  +  3.76  (x  +  y/4)  N2  (3) 
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Temperature 


Figure  5.  Progress  of  Reaction 
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where 


0 


equivalence  ratio 

fractional  oxygen  consumption  efficiency 


t  = 


Utilizing  Eq(3),  the  mole  fractions  of  fuel  and  oxygen  may  be  written 

_ (1-t)0  ■ _ 


Xf  =  4.76  (x  +  y/4)  +  0  +  «0  (y/4 -1) 

(x  +  y/4)  (1  -  t0) _ 

X°  =  4.76  (x  +  y/4)  +  0  +  e0  (y/4 -1} 

Substitution  of  Eqs(4a)  and  (4b)  into  Eq(2)  yields 


(4a) 

(4b) 


VP  2 


h.  Je,:E/RTrl  A _ (1 -€)_(_! -<0)  _ 1 

r2  (Jr3/2  J  (5  (4 . 76 ( x+y/4 )  +  0  +  1 0  (y/4  -  1 ))  2J 


(5) 

If  the  contributions  of  the  smaller  terms  0  and  t<t>  (y/4  -  1)  appear¬ 
ing  in  the  denominator  are  neglected,  one  may  write 


m„ 


e-E/RTr  (i_<)  (i_t0) 


VPv 


3/2 


(6) 


r  -r 

For  rich  mixtures,  a  similar  expression  may  be  obtained  as  follows 


ny  e-E/RT 

V*  & 


VPr^ 


Tr3/2 


r  0  (1-Q‘ 

t 


(7) 


Thus  a  general  expression  is 

mr  e~E/RTr  0  (1-e)  (l-<j) 

VPr2  Tr  3X2 

where,  for  0  <  1 ,  j  =  0 

and  for  0  > 1 ,  j  =  1 


(8) 


The  term  — ^  has  come  to  be  known  as  the  air  loading  parameter  0 
VPr2 

of  the  reactor;  for  a  given  fuel,  this  term  is  a  function  of  the 
variables  0,  E,  and  Tr.  The  reaction  temperature  Tr,  in  turn, 


is  determined  by  T,  0 ,  and  t.  Herbert  (Refs  21  and  22)  has  exten¬ 
sively  characterized  the  performance  of  stirred  reactors.  In  Fig  6, 

mr 

the  various  relations  calculated  by  Herbert  between  ■  -*  and  <t> ,  <, 

VPf 

and  T  are  shown  schematically.  Thus  in  Fig  6(a),  it  is  shown  that  for 
a  given  0  and  T  the  loading  parameter  <t>  rises  to  a  local  peak 
value  0max  with  increasing  c  and  then  decreases  to  zero  as 
approaches  1.0.  For  the  condition  shown,  0max  occurs  at  a  high  value 
of  efficiency  (opt  (*opt  >0.80  approximately).  Furthermore,  it  can 
be  shown  from  heat  balance  considerations  that  stable  operation  of  the 
reactor  can  only  be  obtained  in  the  region  t*  <  (  <  1 ,  where  t*  =  topt- 
The  significance  of  0max  to  the  phenomenon  of  flame  stabilization 
is  that,  for  a  given  reactor  (flameholder)  configuration,  operating  at 
fixed  inlet  conditions  (T,  0,  and  P),  0max  determines  the  maximum  mass 
flow,  and  correspondingly,  the  maximum  flow  velocity  that  the  reactor 
can  sustain  before  extinction.  It  will  also  be  noted  that  for  a  given 
inlet  temperature,  0,^  increases  with  equivalence  ratio  in  the  weak 
mixture  regime,  and  decreases  with  equivalence  ratio  in  the  rich  mixture 
regime.  The  peak  value  of  0max  occurs  near  the  stoichiometric  equiva¬ 
lence  ratio.  Thus,  reactor  operation  is  constrained  between  a  weak  and 
rich  "blow-out"  limit  as  shown  in  Fig  6(b).  Also,  it  will  be  noted  that 
as  the  reactor  loading  is  increased,  the  range  of  equivalence  ratios 
over  which  stable  operation  can  be  maintained  is  diminished  with  both 
rich  and  lean  limit  equivalence  ratios  approaching  the  stoichiometric 
value,  <t>  =  1.  In  Fig  6(c)  it  will  also  be  observed  that  for  a 

given  0,  an  increase  in  inlet  temperature  (T)  causes  a  significant 
increase  in  the  level  of  0  and  correspondingly  of  the  value  of 
^max* 
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E 


S>(lb/sec  ft  atm  ) 

(a)  Relation  Between  <t> ,  and  <t>  (Constant  T) 


Rich  Blow-out 


->>n  qn  O^S^. 


.ean  Blow-out 


$(lb/sec  atm  ft  ) 

(b)  Rich  and  Lean  Blow-out  Limits  (Constant  T) 


=  1.0 


0  <t>(lb/sec  ft3  atm2) 200 

(c)  Relation  Between  <t>,  t,  and  T  (Constant  <t>) 

Fiqure  6.  Characteristic  Curves  of  Stirred  Reactor 
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It  is  convenient  at  this  point  to  introduce  a  modified  form  of  air 
loading  parameter.  The  blow-out  air  loading  0max  is  dependent  on  e, 
4>,  and  T.  Bragg  (Ref  11)  proposed  that  the  dependence  of  0max  on 
T  could  be  represented  by  a  relation  of  the  form  0max  <*  eT/b,  where 
b  depends  on  4> .  Subsequently,  Herbert  (Ref  22)  proposed  the  follow¬ 
ing  empirical  equation  for  b,  for  weak  mixture  conditions. 

b  =  24 5  (1.39  +  Ln  0/1.03),  0<1.O3  (9) 

This  relation  for  b  corresponds  to  temperature  (T)  expressed  in 
degrees  Kelvin. 

A  modified  air  loading  parameter  (&)  is  now  defined  as 

(10) 

The  validity  of  the  stirred  reactor  analysis  was  soon  confirmed  by 
the  testing  of  idealized  research  reactors.  Typically,  these  experi¬ 
mental  reactors  were  of  spherical  design  to  provide  minimum  surface  to 
volume  ratio,  of  small  scale,  and  were  operated  at  near-adiabatic  condi¬ 
tions.  The  dominant  design  feature,  illustrated  in  Fig  7(a),  was  the 
use  of  a  large  number  of  jets  to  create  intense  mixing  in  the  reaction 
zone.  Obviously,  a  high  pressure  drop  was  required  to  produce  these 
jets,  and  such  pressure  losses  could  not  be  tolerated  in  practical 
combustion  systems.  The  principal  experimental  aim  of  these  researches 
was  to  determine  the  "blow-out"  conditions  for  both  lean  and  rich  mix¬ 
tures  at  various  initial  conditions.  By  comparing  experimental  blow-out 
data  with  theoretical  estimates  based  on  equations  simulated  to  Eq(5), 
it  was  possible  to  deduce  an  overall  reaction  order  (n)  for  the  "global" 


r  - - - * - 

reaction  assumed,  for  appropriate  choices  of  rate  constant  (k-| )  and 
the  apparent  activation  energy  tern  (E).  The  classic  experimental  work 
in  this  area  was  performed  by  Longwell  and  Weiss  and  is  reported  in 
Ref  23:  Fig  7(b)  is  taken  from  their  data  and  confirms  the  excellent 
agreement  between  experimental  data  and  theoretical  prediction  which 
can  be  obtained  by  appropriate  choice  of  n,  E,  and  k-j .  It  is 
interesting  to  note  that  for  the  high  temperature  combustion  of  hydro¬ 
carbons  in  air,  an  overall  reaction  order  of  1.8  with  an  apparent 
activation  energy  of  42  KCal/Mol  was  suggested. 

Although  much  work  has  been  done  to  elucidate  the  detailed 

kinetic  and  mixing  mechanisms  associated  with  hydrocarbon  combustion 

in  stirred  reactors,  the  basic  one-step  model  has  continued  to  be  of 

utility,  and  some  useful  extensions  of  this  model  have  been  made  in 

recent  years.  One  such  refinement  has  been  an  improved  estimate  of  the 

rhv% 

'  pressure  exponent  n  occurring  in  the  loading  term  — ^  From  the 

Vr 

■  time  of  the  early  work  of  Longwell  and  Weiss  (Ref  23),  n  has  usually 

been  assumed  to  lie  in  the  range  1.75  to  2.00.  Indeed  Greenhough  and 
Lefebvre  (Ref  6)  utilized  n  =  1.75  for  kerosene-air  mixtures  with 

very  satisfactory  results.  However,  Hottel  et  al.  (Ref  24),  in  subse¬ 
quent  reactor  studies,  argued  that  n  varied  with  equivalence  ratio; 
they  presented  data  which  were  compatible  with  n  varying  from  1.7 
at  stoichiometric  conditions  down  to  1.0  at  an  equivalence  ratio  of  0.47. 
An  average  value  of  1.3  was  reported  for  these  experiments.  It  should 
also  be  noted  that  this  work  was  performed  using  a  non-standard  fuel, 

the  empirical  formula  for  which  was  C0  H  .  0  .  rtc. 

2.04  4.62  0.85 

It  is  appropriate  at  this  point  to  draw  attention  to  the  extended 
I  experimental  efforts  of  Hottel  et  al.  (Ref  25)  and  of  Clarke  and  his 
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associates  (Refs  26  and  27)  in  the  area  of  stirred  reactors.  These 
investigators  refined  and  extended  the  foundational  work  of  Longwell 
and  Weiss  (Ref  23). 

Subsequently,  Kretschmer  and  Odgers  (Ref  28)  carried  out  an 
extensive  study  of  data  obtained  from  propane  fueled  stirred  reactors; 
they  postulated  an  empirical  reaction  rate  equation  based  on  a  one- 
step  model.  In  this  model  the  value  of  the  exponent  n  varied  with 
equivalence  ratio  in  the  following  manner: 

Weak  mixtures  n  =  2 0 
Rich  mixtures  n  =  2/0 

and  the  individual  reaction  orders  for  the  fuel  and  oxygen  were  assumed 
to  be  the  same  and  thus  equal  to  n/2.  In  a  subsequent  paper,  Odgers 
and  Carrier  (Ref  29)  modified  the  above-recommended  exponents  as 
follows: 

Range  of  0  Recommended  Exponent 

<t>  <  0.5  1,0 

0.5  <  0  <  1.0  2  0 

1.0  <  0  <  2.0  2/0 

0  >  2  1.0 

Thus  in  no  case  does  the  value  of  n  become  less  than  unity. 

The  basic  mass  balance  assumed  (Ref  28)  was 

0  C3Hg  +  50 2  +  5m  N2  -v  3  je  C02  +  4  je  H20  + 

(0 -  jO  C^Hg  +  5(l-jt)  02  +  5mN2  (11) 

and  the  corresponding  loading  equation  becomes 
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(12) 


E 

\  k1  (m+1)  [5(l-j<)]n/2  [0-jt]n/2e  RTr 
VPrn  Rn  U  [5(m+l)  +  +  j<]n  [Tr](n"°-5^ 

where,  for  4>  <  1,  j  =  <t> 
and  for  <t>  >  1,  j  =  1 

However,  it  should  be  noted  that  the  use  of  a  semi-empirical 
reaction  rate  equation  led  to  substantial  variations  in  the  "effective" 
activation  energy  term  E/R.  Despite  these  substantial  variations  in 
activation  energy,  the  model  proposed  by  Kretschmer  and  Odgers  does 
harmonize  a  large  amount  of  data  concerning  propane  fueled  systems. 

To  apply  a  quotation  from  Levy  and  Weinberg  (Ref  30:252),  "...  this 
may  be  considered  rather  fortunate  as  it  tends  to  justify  both  the 
chemists'  contempt  for,  and  the  engineers'  use  of  the  global  reaction 
concept  each  within  their  respective  fields  of  application." 

In  a  related  development,  a  recent  paper  by  Ball al  and  Lefebvre 
(Ref  31)  modelled  the  weak  extinction  limits  of  propane  flames  sta¬ 
bilized  on  hollow-cone  flameholders  using  a  simple  one-step  reaction 
approach.  In  formulating  the  analysis,  the  effect  of  equivalence  ratio 
<t>  on  the  exponent  n  was  acknowledged;  however,  the  value  of 
n  =  2<t>  previously  recommended  (Refs  28  and  29)  was  not  used.  Instead, 
a  value  of  n  =  1.25  was  assumed  which,  together  with  a  fixed  activa¬ 
tion  energy  of  54,000  cals/g  mole,  yielded  an  excellent  prediction  of 
weak  blow-out  limits  over  equivalence  ratios  ranging  from  0.4  to  0.8. 

Thus  the  general  trend  of  reducing  n  at  lower  equivalence  ratios  has 
been  found  useful  in  practical  flameholder  testing  utilizing  propane 
as  fuel.  The  application  of  basic  stirred  reactor  theory  to  the  predic¬ 
tion  of  weak  extinction  limits  of  flameholders  will  now  be  considered  in 


more  detail . 
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Application  of  Stirred  Reactor  Model  to  Flame  Stabilization 


Early  designs  of  ramjet  engine  combustors  often  incorporated  a 
single  stabilizing  element  such  as  a  cone.  The  flame  pattern  corre¬ 
sponding  to  this  design  is  illustrated  in  Fig  8(a);  for  this  case,  it 
is  possible  to  separate  the  combustion  process  into  two  elements:  flame 
stabilization  and  flame  propagation  (or  flame  spreading).  Much  early 
research  was  therefore  devoted  to  the  flame  stabilizing  capabilities 
bluff  bodies  under  various  boundary  conditions  such  as  free-jet  (con¬ 
stant  pressure)  operation  or  confined  tube  (constant  area)  operation. 
Examples  of  such  early  work  are  well  documented  (Refs  32  and  33).  In  a 
similar  fashion,  flame  propagation  became  another  area  of  study. 

Obviously  the  configuration  shown  in  Fig  8(a)  represents  an 
inefficient  use  of  available  volume.  Consequently,  combustor  designs 
evolved  using  annular  stabilizing  elements  and,  for  more  controlled 
flow  distribution,  the  "can"  burner  concept  was  introduced.  These  con¬ 
cepts  are  shown  in  Figs  8(b)  and  8(c). 

It  is  apparent  that  in  contrast  to  the  system  shown  in  Fig  8(a), 
where  a  simple  description  of  the  flame  spreading  process  could  be 
given  in  terms  of  spatial  coordinates  and  a  turbulent  flame  propagation 
rate,  the  systems  shown  in  Figs  8(b)  and  8(c)  essentially  require  a 
"volumetric"  description.  This  volumetric  description  could  be  sub¬ 
divided  into  a  stabilizing  or  primary  zone  plus  a  secondary  zone  where 
combustion  is  effectively  completed.  The  point  to  be  made  is  that  the 
description  of  the  process  must  be  related  to  the  configuration  of  the 
combustor. 

For  any  of  the  geometries  shown  in  Figs  8(a),  (b),  and  (c),  the 
stabilizing  zone  is  a  region  of  intense  mixing  which  may  be  modelled, 
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at  various  levels  of  sophistication,  as  a  stirred  reactor.  Thus  in 
Fig  9,  the  recirculation  zone  (RZ)  behind  a  bluff  body  is  depicted. 
If  this  RZ  is  considered  as  a  stirred  reactor,  then  Eq(10)  could  be 
applied  to  this  zone  and  the  corresponding  blow-out  limits  described 
by  a  relation  of  the  form 


*LB0  =  f  (*) 


(13) 


In  order  to  describe  the  RZ  in  these  terms,  two  key  quantities  must  be 
known;  namely,  the  air  mass  flow  into  the  RZ,  and  V,  the  volume  of  the 
RZ.  For  a  typical  bluff  body  immersed  in  a  flowing  homogeneous  fuel- 
air  mixture,  it  is  generally  assured  that  mr  is  proportional  to  the 
density  of  the  stream  (p),  the  surface  area  of  the  RZ  (S),  and  the 
flow  velocity  at  the  edge  of  the  body  (U),  The  further  assumptions  are 

Hence 

I 

2, 


?  3 

made  that  S«Df  and  that  1/  a  Df 


m  pDf  U 

?  a  3  o  a 

V  p/  DfV 


U 


DfPT 


(14) 


where  Pr  has  been  put  equal  to  P, 
From  Eq(10),  at  blow-out  conditions 

m_ 


= 


VP.. 


1  „  u 

J7F  “  d 7 


T  e 


1 

T7F 


Thus  for  a  fixed  inlet  temperature 


I If  a 


u 


DfP 


(15) 


(16) 


A  number  of  assumptions  have  been  made  in  order  to  estimate  the 
quantities  mr  and  V.  The  validity  of  these  assumptions  will  not 
be  discussed  at  this  point,  but  will  be  examined  later  in  the  context 
of  coaxial  dump  burners.  It  should  be  noted,  however,  that  in  a 


26 


Figure  9.  Bluff  Body  Recirculation  Zone 


foundational  investigation  of  the  flame  stability  characteristics  of 
circular  discs,  published  in  1955,  DeZubay  (Ref  32)  was  able  to 
correlate  the  experimental  data  using  a  loading  parameter  of  the  form 


\Jj  a 


U 

p0.95 


07) 


Since  that  time  a  vast  number  of  investigations  of  flame  stability  have 

I  I  01 

successfully  utilized  parameters  of  the  form  u  _  to  correlate 

P^Df  y 

isolated  bluff  body  data  for  constant  inlet  temperature  conditions. 

The  effect  of  varying  inlet  temperature  on  stability  has  also  been 
introduced  in  a  direct  exponent  form  rather  than  in  the  factor 
T  eT/,b  of  Eq(  1 5) ;  thus  another  form  of  loading  parameter  is 


a 


p/»Df>T<5 


08) 


The  temperature  exponent  d  cannot  be  expected  to  be  strictly 
constant  because  the  appropriate  temperature  variation  is  contained  in 
the  factor  T  eT^b.  Thus  <J  can  be  shown  to  depend  on  both  the 
equivalence  ratio  {<f>)  and  the  range  of  initial  temperatures  considered. 
However,  over  limited  ranges  of  <j>  and  temperature,  the  use  of  an 
exponent  <5  is  acceptable  for  practical  correlations.  In  a  typical 
correlation  quoted  by  Stull  et  al .  (Ref  1),  a  value  of  4  equal  to 


*  Note  that  1)  and  P  in  this  expression  correspond  to  the  values  at 
flame  blow-off:  these  should  strictly  be  denoted  by  Ugg  and  Pgg.  How¬ 
ever,  to  avoid  complex  notation  this  distinction  will  not  be  observed 
when  the  context  is  clearly  addressing  blow-off  conditions. 
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1.5  is  given  and,  in  general,  values  of  6  fall  in  the  range  from 
1.0  to  2.0. 

Up  to  this  point  the  stability  of  the  recirculating  flow  region 
has  been  discussed  solely  in  terms  of  homogeneous  reactor  chemistry. 
The  related  influence  of  the  aerodynamics  of  the  flow  on  the  stabiliz¬ 
ing  process  has  been  considered  only  in  the  determination  of  simple 
estimates  jf  the  mass  flow  entering  the  recirculation  zone  and  the 
corresponding  zone  volume.  However,  it  will  be  apparent  that  the 
fundamental  processes  determining  flame  stability  are  the  transport 
phenomena  occurring  in  the  complex  flow  field  associated  with  the 
flameholder.  These  phenomena  are  dependent  on  both  the  overall  flow 
field  generated  by  the  flameholder  and  the  finer  details  of  the  flow 
field.  A  key  factor  in  determining  the  overall  flow  field  is,  of 
course,  the  overall  Reynolds  number  of  the  flow.  The  states  of  the 
boundary  layer/wake  flows  associated  with  the  flameholder  are  also 
important.  In  particular,  the  flame  stability  characteristics  will 
be  affected  as  the  state  of  the  wake  changes  from  laminar,  through 
transitional,  to  fully  turbulent.  For  most  practical  applications, 
the  flows  will  be  fully  turbulent. 

The  interrelation  between  the  chemical  and  aerodynamic  factors 
has  generally  been  acknowledged  by  assuming  that  for  a  given  fuel/air 
ratio  the  performance  of  a  combustion  system  is  dependent  on  the 
loading  (#)  and  the  associated  Reynolds  number  (Rg).  In  particular, 
Herbert  (Ref  34)  has  assumed  that 

Xlf  a  (Re)-f  (19) 

where  for  fully  turbulent  flows  f  has  a  value  in  the  region  of  0.05 
to  0.10. 
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One  of  the  finer  details  of  the  flow  field  is  the  state  of  the 
approach  stream  turbulence;  this  factor  can  influence  the  flame  sta¬ 
bility  to  varying  degrees  depending  on  the  geometry  and  flow  conditions 
of  the  flameholder.  The  effect  of  free  stream  turbulence  characteristics 
on  the  flow  field  associated  with  a  bluff  body  is  an  emerging  field  of 
study  (Ref  35),  and  no  generalized  comments  can  be  usefully  summarized 
at  this  time.  However,  in  the  recent  investigation  by  Ballal  and 
Lefebvre  (Ref  31)  of  the  stability  of  hollow  cone  stabilizers,  it  was 
shown  that  the  effect  of  free  stream  turbulence  intensity  could  be 
accounted  for  by  simply  modifying  the  local  velocity  term. 

Another  aerodynamic  phenomenon  of  considerable  significance  is 
that  of  interaction  of  the  flameholder  flow  field  with  the  boundary  of 
the  surrounding  flow.  The  blockage  produced  by  the  presence  of  the 
flameholder  in  a  pipe  results  in  an  acceleration  of  the  flow  as  it 
approaches  the  edge  of  the  flameholder  and  simultaneously  affects  the 
flow  field  structure  downstream  of  the  flameholder.  Hence,  the  degree 
of  blockage  is  a  factor  of  considerable  importance.  Thus  Wright  (Ref  36) 
found  that,  for  a  flat  plate  flameholder,  the  length  of  the  recirculation 
zone  (L)  varied  inversely  as  the  square  root  of  the  blockage  ratio 
(flameholder  area/duct  area).  Interestingly,  it  was  also  found  that 
the  width  of  the  flameholder  wake  (W)  varied  inversely  as  the  square  root 
of  the  blockage  ratio.  Thus  the  ratio  L/W  was  found  to  remain  approxi¬ 
mately  constant  even  though  L  and  W  varied  markedly  with  blockage. 

It  follows  from  the  above  that  the  blockage  ratio  is  an  important  aero¬ 
dynamic  factor  affecting  the  flame  stability  of  bluff  body  flameholders. 
It  should  also  be  noted  that  because  of  the  flow  deflection  produced  by 
some  bluff  body  shapes,  the  effective  aerodynamic  blockage  may  be 

significantly  greater  than  the  corresponding  geometric  blockage  (Ref  37). 
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To  recapitulate,  it  will  be  appreciated  that  the  analysis  of 
flame  stability  phenomena  is  a  most  complex  task.  The  description  of 
flameholder  performance  by  a  parameter  of  the  form  shown  in  Eq ( 1 8 )  is 
an  attempt  to  account  for  the  major  controlling  factors  in  a  simpli¬ 
fied  manner.  Consequently,  the  set  of  exponents  a,  0,  y,  6  will 
generally  assume  different,  but  nevertheless  bounded  values  which 
depend  on  the  specific  flameholder  configuration,  fuel,  and  approach 
flow  conditions. 

It  is  convenient  at  this  point  to  draw  attention  to  a  recent 
attempt  to  correlate  flame  stability  limits  by  Plee  and  Mellor  (Ref  38), 
utilizing  a  parameter  denoted  as  the  "inverse  theta  parameter"  (0^). 
The  origin  of  the  basic  9 L  parameter  may  be  found  in  the  early  work 
of  Greenhough  and  Lefebvre  (Ref  6),  who  were  seeking  a  parameter  to 
describe  the  efficiency  of  gas  turbine  combustors.  A  full  derivation 
of  this  parameter  was  given  subsequently  by  Lefebvre  (Ref  39),  and  some 
early  correlations  were  displayed  by  Lefebvre  and  Halls  (Ref  40'.  In 
principle,  the  combustion  process  is  modelled  as  a  turbulent  flame 
brush.  An  energy  balance  then  leads  to  the  result  that  combustion 
efficiency  is  proportional  to  the  ratio  of  turbulent  flame  speed  (Sy) 
to  a  reference  velocity  (Urey)  in  the  combustion  chamber. 


Damkohler's  relation  (Ref  41)  between  turbulent  flame  speed 
laminar  flame  speed  (Su),  and  Reynolds  number,  Sy  a  Su 
then  employed,  yielding 


(20) 

(sT), 

ReC’  is 


S  R  c 
u  e 


'ref 


(21) 
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A  second  order  reaction  process  is  then  assumed,  which  implies  that 


Su  is  independent  of  pressure.  The  effect  of  initial  temperature  is 
neglected  at  this  point,  leading  to  a  semi -empirical  relation  which  is 
essential  1 y 


^  =  f 


pl-75  A  D  0.75 
r  Mref  uref 


(22) 


where  the  terms  A  ,,  and  D  are  a  reference  combustor  area  and 
ret  ref 

diameter,  respectively. 

An  initial  temperature  correction  is  now  introduced  into  the 
expression,  yielding 


=  f  ( 6  L)  where 


P1-75*  .0  f°-75eT/l’ 

ref  ref 


(23) 

Thus  for  a  given  combustor,  the  combustion  efficiency  is  related  to  the 
pressure  (P)  and  temperature  (T)  conditions  at  entry  and  the  corre¬ 
sponding  mass  flow  (m).  The  "theta"  parameter  has  been  very  success¬ 
fully  utilized  for  the  correlation  of  combustion  efficiency  data  and  for 
the  initial  design  of  gas  turbine  combustors.  The  application  of  the 
theta  parameter  to  the  correlation  of  flame  stability  is  an  attempt  to 
extend  its  utility.  In  this  case,  the  reciprocal  of  the  theta  parameter, 
termed  the  "inverse  theta  parameter,"  is  a  more  suitable  form,  and  for 
a  given  combustor  it  may  be  expressed  as 


1  a  p0.75UTeT/b  (24) 

The  inverse  theta  parameter  can  also  be  expressed  as  «  U^/P^T4, 

the  form  conventionally  used  for  correlating  flame  stability  data. 


parameter 


However,  there  was  little  direct  evidence  that  the  6^ 
would  be  effective  in  correlating  flame  stability  data,  although  partial 
success  was  achieved  (Ref  6).  Plee  and  Mellor  (Ref  38)  subsequently 
found  that  use  of  the  0^  parameter  did  indicate  stability  trends 
(i.e.,  the  weak  extinction  limit  ^BO  increasec)  with  but  good 

correlations  were  not  generally  obtained.  However,  it  would  be  pre¬ 
mature  to  discard  the  inverse  theta  parameter  at  this  time  pending 
further  experimental  work. 

Up  to  this  point  the  phenomenon  of  flame  blow-off  has  been 
approached  solely  from  the  point  of  view  of  stirred  reactor  theory. 
However,  other  theories  of  flame  stabilization  have  been  proposed,  and 
some  are  of  technological  interest.  In  the  next  section  a  brief  review 
of  such  theories  will  be  given. 

Alternative  Criteria  for  Flame  Stabilization 

Historically,  flame  stability  criteria  were  evolved  to  describe 
the  stable  operating  range  of  laboratory- type  burners.  These  criteria 
were  sought  employing  such  relatively  fundamental  combustion  concepts 
as  the  laminar  and  turbulent  flame  speeds  associated  with  pre-mixed 
combustibles.  A  significant  starting  point  for  this  discussion  is  the 
work  of  Putnam  and  Jensen  (Ref  42),  in  which  was  derived  a  Peclet  num¬ 
ber  criterion  for  predicting  "flash  back"  of  laboratory-type  burners. 
These  investigators  used  a  conventional  velocity  gradient  approach 
which  basically  assumes  that  stable  operation  occurs  when  flow  velocity 
(U)  and  laminar  flame  propagation  velocity  (Su)  are  matched:  their 
criteria  for  equilibrium  was  that 

Peu  -  k2  (Pes)2  (25) 
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where  Peu  is  the  Peclet  number  based  on  average  flow  velocity  (-— ) , 

SuDf 

P  is  the  Peclet  number  based  on  laminar  flame  speed  (— — ), 
es  3 

0 ^  is  a  characteristic  dimension  of  flameholder,  and 
a  is  the  thermal  diffusivity. 

This  flash-back  criteria  was  also  applied  to  describe  blow-off  data  of 
small  axial  rod  flameholders;  a  good  correlation  was  obtained  with  an 
expression  of  the  form  given  in  Eq(25).  An  important  observation  made 
by  Putnam  and  Jensen  was  that  Eq(25)  could  be  expressed  as 


k  S 
k2  bu 


In  the  above  expression,  the  term  may  be  regarded  as  a 

a 

characteristic  flow  time  (  rdu)  and  the  term  can  be  regarded  as 

a  characteristic  reaction  time  (t  ).  Thus  the  Peclet  number  cri- 

aS 


terion  can  also  be  written  as  r. 


i£  ras ’  and  if  Tdu  is  less 


than  Tas,  blow-off  of  the  flame  will  occur.  In  general,  it  would 
be  anticipated  that  t  c  would  be  a  function  of  fuel-type,  initial 
temperature  and  pressure,  and  fuel-air  ratio.  The  concept  inherent  in 
Eq (26 )  is  a  balance  between  a  time  associated  with  the  flow  (  rd(j)  and 
a  time  characteristic  of  the  reaction  (r  ):  this  concept  will  appear 

aS 

repeatedly  in  this  review. 

The  basic  Peclet  number  criterion  was  placed  on  a  more  general 
theoretical  basis  by  Spalding  and  Tall  (Ref  43);  for  high  values  of 
Peclet  number,  their  criterion  for  flame  stabilization  was  written 

i  i  n  re-  2 
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which  is  of  the  same  form  as  Eq ( 25) .  However,  it  must  be  noted  that 
these  authors  introduced  the  laminar  flame  speed  into  their  analysis  as 

O 

a  measure  of  the  reaction  rate  (<*  S  )  and  not  as  an  index  of  flame 
propagation  rate.  The  form  of  correlation  of  data  anticipated  by 
Eq ( 27 )  is  shown  schematically  in  Fig  10.  Acceptable  data  correlations 
have  been  published  (Ref  44);  however,  additional  data  acquired  since 
the  publication  of  this  reference  have  not  been  satisfactorily  corre¬ 
lated  as  shown  by  Shchetinkov  (Ref  45).  In  order  to  utilize  the 
criterion  of  Eq(27),  the  Peclet  number,  based  on  the  laminar  flame 
speed,  has  to  be  evaluated.  Unfortunately,  neither  Su  nor  thermal 
diffusivity  values  are  generally  available  for  gas  mixtures  at  condi¬ 
tions  of  technological  interest;  consequently,  this  particular  criterion 
has  not  found  wide  application. 

Another  flame  stability  criterion  which  deserves  brief  mention  is 
that  of  extinction  due  to  "flame-stretch."  This  concept  also  arose  in 
connection  with  laboratory  studies  of  turbulent  flames  and  was  originated 
by  Karlovitz  et  al.  (Ref  46).  It  was  observed  that  flame  propagation 
was  arrested  in  regions  of  flow  where  steep  velocity  gradients  existed. 
The  corresponding  increase  in  flame  surface  area,  or  "flame-stretch"  in 
the  steep  gradient  area,  was  quantified  using  a  factor  termed  the 
"Karlovitz  number  (K)"  defined  as 


dp  dU 
U  d£ 


(28) 


where  U  (£)  is  the  local  velocity  of  unburned  gas,  ^7-  the  velocity 
gradient,  and  dQ  a  representative  flame  thickness  defined  as 

^o  =  a/%* 
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The  utility  of  the  Karlovitz  number  is  that  a  critical  value  of  K 
may  be  determined  above  which  flame  extinction  occurs.  Lewis  and  von  Elbe 
(Ref  47)  have  applied  the  Karlovitz  criterion  to  bluff-body  flame  sta¬ 
bilization  with  modest  success.  However,  the  validity  of  the  flame- 
stretch  concept  is  still  a  subject  of  controversy  as  exemplified  in  the 
discussion  of  Melvin  and  Moss  (Ref  48).  Also,  like  the  Peclet  number 
criterion,  the  Karlovitz  number  has  found  no  general  application  to 
combustion  processes  of  technological  interest. 

Meanwhile,  in  the  early  Soviet  literature  (Ref  49)  on  bluff-body 
flame  stability,  the  work  of  Dunsky  is  referenced  as  the  origin  of  the 
so-called  Mikhel'son  criterion  (Mi),  defined  as 


Mi  = 


rz 


cz 


(29) 


where  r rz  is  the  mean  particle  residence  time  in  the  recirculation 
zone  (RZ)  of  the  stabilizer  and  r ^  is  a  combustion  reaction  time. 
In  this  case,  r  is  taken  as  tne  characteristic  time  associated 

6  o 

with  the  laminar  flame  speed,  namely  <= — ,  where  <5  is  the  flame 

u  a 

thickness.  Also  <5  Q  is  taken  as  proportional  to  so  that  rcz 
is  defined  as  r_  =  .  Dunsky  used  a  characteristic  dimension 


cz 


>u 


2X  so  that  the  Mikhel'son  criteria  is  defined  as 

T, 
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Shchetinkov  (Ref  45)  observes  that  Mi  displays  a  constant  value 
(  *  0.45)  only  for  limited  experiments  with  small  conical  stabilizers. 
A  related  but  distinct  theory  put  forward  by  Zukoski  and  Marble 


(Ref  50)  in  the  U.S.A.  has  also  been  successful  in  predicting  blow-off 
limits  of  certain  combustor  configurations  and  has  found  wide  application. 


A  suitable  configuration  by  which  to  introduce  the  elements  of 
this  theory  is  the  isolated  stabilizer  shown  in  Fig  9.  In  this  sketch 
the  recirculation  zone  is  seen  to  be  enveloped  by  a  shear  layer  which 
begins  as  an  extension  of  the  boundary  layer  and  becomes  progressively 
thicker  in  the  downstream  direction.  Mass  and  energy  exchange  occurs 
between  the  outer  flow  and  the  recirculation  zone  flow  through  the 
shear  layer.  Under  stable  operating  conditions,  a  flame  is  initiated 
in  this  mixing  layer  which  subsequently  propagates  into  the  surrounding 
flow.  Now,  according  to  stirred  reactor  theory,  this  initial  flame 
would  be  extinguished  when  the  stirred  reactor  attained  blow-out.  How¬ 
ever,  Zukoski  and  Marble  postulated  that  the  only  active  combustion 
associated  with  the  RZ  is  limited  to  the  shear  layer  at  the  boundary, 
and  that  a  portion  of  the  burned  gas  from  this  flame  zone  is  recirculated 
upstream  to  replenish  the  RZ.  The  criterion  of  flame  stability  is  that 
sufficient  reactant  gas  is  ignited  during  its  residence  time  in  the 
shear  layer  to  propagate  a  flame.  Thus,  if  U  is  the  local  flow 
velocity  and  L  a  length  representative  of  the  RZ  length,  then 

=  L/U  is  a  measure  of  the  exposure  of  the  fresh  gas  to  the 
hot  combustion  products  which  serve  as  an  ignition  source.  In  par¬ 
ticular,  a  critical  time  r  .  is  defined  for  ignition  of  the  flow, 
thus 

Tc.  =  4  =  n,  (3,) 

where  Ugg  is  the  measured  blow-off  velocity. 

The  time  t  .  was  denoted  as  the  "characteristic  ignition  time." 
Typical  experimental  values  of  t  .  derived  from  Zukoski  and  Marble 
(Ref  50),  and  from  the  works  of  Wright  (Ref  36)  and  Broman  (Ref  51)  are 
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shown  in  Fig  11.  The  variation  of  t  with  0  and  initial 
temperature  is  shown  for  the  specific  fuels  utilized;  the  effect  of 
pressure  is  not  documented.  Two  aspects  of  the  Zukoski-Marble  model 
are  now  familiar.  Firstly,  the  fluid  dynamics  of  the  problem  are 
decoupled  from  the  chemistry.  depends  solely  on  the  fluid 

dynamics  and  rci  on  the  chemistry.  Secondly,  the  blow-off 
criterion,  r/u/ TC1-  =  can  regarded  as  a  ratio  of  residence 
time  to  chemical  reaction  time  which  is  in  the  form  of  Damkohler's 


first  dimensionless  group  for  chemical  reactors  (Ref  52).  This  resi¬ 
dence  time/reaction  time  approach  has  appeared  in  recent  years  in  many 
forms.  Thus  Mironenko  (Refs  53  and  54)  introduces  a  flow  time  defined 

as  — ,  where  L  is  the  lenath  of  the  RZ  zone  and  Ure] 

rur  ure-| 

is  the  relative  flow  velocity.  Thus  II  ,  =  U  +  l)  where  Um  is 

J  re l  m  p  m 

the  downstream  velocity  of  the  flow  past  the  stabilizer  and  Up  is 

the  upstream  velocity  in  the  RZ  behind  the  flameholder.  Also,  a 

characteristic  reaction  time  7*  is  introduced  related  to  the  tem- 

ce 

perature  rise  required  to  increase  the  associated  reaction  rate  by  the 
exponential  factor  (e).  Using  these  revised  definitions,  t^u  and 
rce,  Mironenko  displays  good  agreement  with  practical  stabilizer 
configurations,  both  in  the  above  references  and  in  a  more  recent  work 
co-authored  with  Solokhin  (Ref  55), 

The  introduction  of  the  recirculation  zone  length  (L)  instead 
of  a  flameholder  dimension  (D^)  as  a  characteristic  length  for  the 
flow  time  has  led  to  increased  complexity.  Typically,  for  turbulent 
flow  past  a  flameholder,  Lisa  rapidly  fluctuating  quantity  and 
experimental  determination  of  a  time-mean  value  of  L  is  difficult. 
Extensive  investigations  of  the  flow  structure  and  associated  L 


values  behind  bluff  body  flameholders  have  been  carried  out  by 
Wright  (Ref  36),  Bovina  (Ref  49),  and  Winterfeld  (Refs  56  and  57). 

These  investigations  have  shown  that  L  is  affected  by  flameholder 
shape,  flow  velocity  and  turbulence,  duct  blockage,  and  heat  release. 
Thus  the  apparent  decoupling  of  flow  and  chemical  phenomena  achieved 
with  the  rdu/  rci  criteria  is  lost  when  L  is  used  as  a  character¬ 
istic  length.  Nevertheless,  Baev  and  Tret'yakov  (Ref  58)  have  success¬ 
fully  correlated  data  using  a  non-dimensional  RZ  length  (U,)  which 
is  corrected  for  the  influence  of  flameholder  shape  and  blockage.  This 
work  is  of  particular  interest  because  the  chemical  time  utilized  is  a 
combustion  time  based  on  experimental  measurements  (Refs  59  and  60)  and 
corrected  for  pressure  variation  by  use  of  the  expression 

rc(p)  =  Tc(p=l)  (32) 

where  rc(p)  is  the  combustion  time  corresponding  to  pressure  P 
and  n  is  the  order  of  the  reaction  (quoted  as  1.6  for  gasoline-air 
mixtures  and  1.2-1. 3  for  propane-air  mixtures).  The  characteristic 
combustion  times  for  gasoline  utilized  by  Baev  and  Tret'yakov  (Ref  59) 
are  illustrated  in  Fig  12  and  compared  to  the  ignition  times  obtained 
experimentally  by  Zukoski  and  Marble  (Ref  50).  Obviously,  both 
sources  yield  curves  of  similar  shape  but  of  differing  magnitudes. 

In  deriving  their  correlation,  Baev  and  Tret'yakov  (Ref  58) 
utilized  a  quenching  distance  criterion  where  the  flame  thickness  at 
blow-off  was  related  to  the  transverse  flameholder  size.  The  final 
correlating  parameter  was 
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As  indicated  earlier,  Winterfeld  (Ref  57)  carried  out  a 

thorough  study  of  flow  phenomena  associated  with  recirculation  zones 

behind  various  shapes  of  bluff  body.  He  determined  both  the  size  k-mi 

shape  of  the  recirculation  zones,  the  associated  velocity  field,  and 

the  corresponding  mean  residence  times  for  both  cold  and  burning  flows. 

He  was  unable  to  relate  blow-off  limits  to  the  mean  residence  time  in 

the  zone  and  concluded  that  flameholding  was  controlled  by  a  flow 

time/reaction  time  criterion.  Winterfeld  considered  that  flameholding 

would  only  be  obtained  if  those  flow  elements  which,  entering  the 

downstream  end  of  the  RZ  and  flowing  upstream,  attained  ignition  before 

they  emerged  at  the  trailing  edge  of  the  flameholder.  In  this  case  the 

characteristic  residence  time  could  be  written  rn  ,  where  llr7 

rz 

is  a  mean  upstream  velocity  along  the  axis  of  RZ.  This  latter  velocity 
was  assumed  proportional  to  the  maximum  upstream  velocity  along  the  RZ 
axis,  Urzmax. 

Finally,  the  appropriate  residence  time  rluTZ  was  expressed  as 


TL 


rz 


(34) 


The  term  U/Urz  is  of  the  order  0.3-0. 5  for  typical  bluff  body 

I,  id  A 

shapes. 

In  a  subsequent  work,  however,  Winterfeld  (Ref  61)  returned  to 
a  simpler  criterion  and  employed  the  ratio  of  a  characteristic  stay 
time  T£u  to  a  characteristic  reaction  time  r  to  describe  sta¬ 
bility.  Significantly,  he  has  pointed  out  that  such  times  cannot  be 
clearly  definable  when  so  many  physio-chemical  processes  are  occurring 
concurrently.  However,  the  stay  time  T^u  was  broadly  defined  as 
the  ratio  of  an  appropriate  length  and  a  related  velocity:  typically 


J 
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T^u  =  L/U.  In  selecting  a  characteristic  time  for  the  reaction, 
Winterfeld  effectively  considered  the  turbulent  flame  as  a  wrinkled 
laminar  flame  with  the  corresponding  reaction  time  taken  as 


T 


as 


a 


(35) 


where  <5q  is  equal  to  a  representative  flame  thickness.  It  was 
further  assumed  that  <5  Q  a  a/S  ,  where  a  is  an  effective  thermal 
diffusivity.  Thus,  for  the  characteristic  reaction  time  r  ,  one 

aS 

has 


T 


as 


a 


(36) 


and  the  ratio  t n  /  -r  can  be  written 

LU  aS 


L  _u_ 
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(37) 


as  “a 

Winterfeld  has  termed  the  right-hand  side  of  Eq(37)  the  characteristic 


factor  KpL  for  flame  stabilization. 
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which  is  a  useful  parameter  for  comparing  the  flameholding  performance 
of  different  fuels.  Winterfeld  further  assumes  that  Su  «  P’°’^ 

(which  implies  a  reaction  order  of  1.8)  and  that  the  thermal  diffusivity 
is  proportional  to  pressure;  these  assumptions  yield  an  overall  stability 
criterion  given  by 
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Yet  another  approach,  which  at  first  sight  is  similar  to  the 
Zukoski  and  Marble  theory,  is  the  recently  introduced  characteristic 


time  method  of  analysis  due  to  Mellor  (Ref  62).  In  principle,  Mellor 
has  divided  the  typical  combustion  process  into  basic  phenomena  such 
as  hetereogeneous  processes,  fluid  mechanical  effects,  and  chemical 
reactions.  Each  of  these  elements  is  characterized  by  an  appropriate 
time  scale,  and  these  characteristic  times  are  applied  to  those  regions 
of  the  combustor  where  a  particular  phenomenon  is  dominant.  In 
particular,  for  the  case  of  the  flameholder,  the  stabilization  process 
is  viewed  as  a  competition  between  the  ignition  delay  time  of  the  fuel 
(  Thc)  and  the  turbulent  mixing  time  in  the  shear  layer  (  Ts^)  that 
exists  between  the  recirculating  burner  gases  and  the  surrounding 
free  stream.  The  stability  criterion  is  thus  Ts//Thc  =  ^ '  TVie 

characteristic  times  utilized  in  this  analysis  are  basically  order-of- 
magnitude  estimates  of  the  processes  which  they  characterize.  Never¬ 
theless,  their  previous  utilization  has  led  to  good  predictions  of 
gaseous  emissions,  and  more  recently  Plee  and  Mellor  (Ref  63)  have 
shown  that  blow-off  limits  of  various  combustor  designs  can  be  corre¬ 
lated  by  appropriate  selection  of  and  r^. 

The  turbulent  mixing  time  in  the  shear  layer  is  assumed  to  be 
equal  to  the  characteristic  decay  time  of  the  large  scale  eddies; 
this  time  is,  in  turn,  related  to  the  ratio  of  turbulence  integral 
length  scale  to  the  turbulent  rms  velocity.  The  assumption  is  then 
made  that  the  integral  scale  is  proportional  to  the  flameholder  width 
(Df)  and  that  the  turbulent  velocity  varies  with  an  "adjusted" 
velocity  (U^)  at  the  lip  of  the  flameholder.  Thus 

ts i "  4  (4°> 
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The  velocity  is  adjusted  to  allow  for  the  fact  that  the  shear 
layer  temperature  (T^)  is  higher  than  the  inlet  temperature  (T), 
using  the  following  equation 


(41) 


The  corresponding  characteristic  time  for  ignition  is  calculated  from 

a  simplified  overall  reaction  rate  expression  as 
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where  (f/a)s  is  the  stoichiometric  fuel-air  ratio.  Since  this  ratio 
is  relatively  constant  for  hydrocarbons  of  interest,  the  term  (f/a)s 
is  eliminated.  Thus  the  criterion  for  blow-off  becomes 
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Once  again  a  decoupling  of  the  fluid  mechanics  and  the  chemistry  is 
achieved.  As  a  concrete  example  of  the  chemical  time,  Plee  and 
Mellor  (Ref  63)  have  utilized  the  data  of  Ballal  and  Lefebvre  (Ref  31) 
which  related  to  the  stabilization  of  pre-mixed  propane-air  flames  to 
obtain 
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where  the  prime  indicates  the  presence  of  the  temperature  ratio  in 
the  right-hand  term.  Both  the  apparent  activation  energy  term  and 
the  pre-exponential  constant  (10~^)  were  determined  empirically.  By 


46 


employing  and  r^c,  a  very  convincing  correlation  of  the  data 

of  Ballal  and  Lefebvre  was  obtained,  described  by  the  expression 

*  2-n  Thc  -  °'46  (45> 

I 

where  r ^  and  T^c  are  in  msec.  This  correlation  is  illustrated 

in  Fig  13.  Furthermore,  successful  stability  correlations  were 
obtained  for  a  variety  of  combustion  systems.  However,  it  must  be 
noted  that  no  pressure  dependence  is  included  in  this  criterion  other 
than  the  minor  influence  of  pressure  on  flame  temperature. 

Up  to  this  point,  two  simplified  models  of  flame  stability  have 
basically  been  considered:  the  stirred  reactor  model  and  the  class  of 
models  described  as  flow  time/reaction  time  models.  In  the  stirred 
reactor  case,  the  stabilization  process  is  modelled  in  an  Eulerian 
manner,  the  homogeneous  reaction  zone  being  fixed  in  space.  For  the 
time-ratio  models,  attention  is  focused  on  the  time  course  of  the 
reaction. 

In  practice,  the  actual  description  of  the  RZ  flow  field  involves 
both  spatial  and  temporal  viewpoints.  Analysis  of  the  processes 
occurring  in  the  RZ  involves  a  number  of  complex  phenomena  embracing 
several  disciplines.  First,  the  state  of  the  approach  flow  needs 
description  as  a  turbulent  field  whose  structure  is  in  general  modified 
as  the  flow  accelerates  around  the  flameholder.  The  flow  which 
separates  from  the  flameholder  lip  entrains  hot  gases  from  the  RZ. 

In  fact,  "engulfment"  appears  to  be  a  better  description  than  entrain¬ 
ment  since  recent  studies  of  mixing  interfaces  have  revealed  large 
islands  or  parcels  of  ambient  fluid  being  drawn  into  the  jet-like  flow. 
These  islands  are  subsequently  "rolled  up"  or  "stretched"  into 
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Figure  13.  Characteristic  Time  Correlation  for  Lean  Blow-off  of 
Bluff  Body  Stabilizer 


’I 


interleaving  layers  of  fresh  mixture  and  hot  products,  resulting  in 
rapid  intermixing  and  reaction  of  these  gases.  Meanwhile,  these 
reacting  parcels  are  moving  downstream,  and  a  portion  of  them  will 
subsequently  be  captured  by  the  ^circulating  flow  and  will  emerge 
later  from  the  RZ  further  upstream  as  a  new  parcel  of  hot  gas  which 
will  in  turn  be  engulfed  by  the  surrounding  flow.  Thus  a  cycle  exists 
of  the  "birth"  of  parcels  of  reacting  mixture,  followed  by  their 
"death"  either  in  the  RZ  or  the  downstream  flow.  Successful  stabiliza¬ 
tion  is  ultimately  dependent  on  the  statistical  outcome  of  this 
"birth-death"  process.  There  is  correspondingly  a  complex  interrelation 
between  the  fluid  mechanics  and  the  chemical  aspects  of  this  process 
which  ultimately  must  be  handled  on  a  statistical  basis.  Obviously, 
treating  this  phenomenon  either  as  a  homogeneous  adiabatic  stirred 
reactor  or  as  a  time-ratio  model  is  a  gross  simplification.  Conversely, 
a  more  exact  analysis  is  a  formidable  task.  A  start  on  this  task  has 
recently  been  made  by  Spalding  (Ref  64),  who  originated  the  so-called 
"ESCIMO"  theory  of  turbulent  combustion.  ESCIMO  stands  for  the  sequen¬ 
tial  events  occurring  in  turbulent  combustion  described  as  Engulfment, 
Stretching,  Coherence,  Interdif fusion,  together  with  a  reminder  that 
the  process  of  interdiffusion  is  analyzed  in  the  reference  frame  of  a 
Moving  Observer  travelling  with  a  parcel  of  fluid.  The  term  coherence 
is  used  to  denote  that  the  parcels  of  intermixed  fluids  cannot  be 
unmixed,  but  maintain  their  "coherence"  as  they  flow  downstream.  The 
successful  evolution  of  the  ESCIMO  theory  should  provide  a  powerful 
tool  to  elucidate  the  mechanism  of  flame  stabilization. 

In  retrospect,  it  is  surprising  that  both  the  stirred  reactor 
approach  and  the  time-ratio  models  both  yield  good  correlations  for 
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specific  flameholders.  This  is  probably  due  to  the  fact  that  the 
parameters  involved  appear  in  the  form  of  ratios  rather  than  absolute 


values.  Thus  in  the  case  of  the  stirred  reactor,  the  air-loading 
mr 

term  _L_ involves  two  quantities,  mr  and  V,  whose  absolute  magni- 

1  mr 

tudes  are  difficult  to  assess.  However,  the  ratio  y—  may  be  a 

meaningful  statistic  describing  the  overall  entrained  flow  to  volume 


ratio  of  a  collection  of  parcels  of  reacting  gas  associated  with  the 


stabilizer.  Similarly,  in  time-ratio  theories  the  ratio  of  residence- 


to-reaction  time  rather  than  the  individual  absolute  values  will 


probably  be  more  meaningful.  The  question  as  to  whether  the  stirred 
reactor  or  the  time- ratio  model  better  represents  physical  reality 
has  been  previously  argued  (Refs  65  and  66)  on  the  basis  of  limited 
experimental  data.  As  noted  earlier,  both  approaches  are  gross 
simplifications,  and  this  discussion  will  not  be  enlarged;  however, 
it  is  appropriate  to  point  out  that  the  Zukoski-Marble  approach  could 
be  embraced  by  reactor  theory  by  considering  the  shear  layer  as  a 
"zone-reactor"  and  modelling  the  flow  behind  a  bluff  body  as  a  reactor 
network  rather  than  a  single  stirred  reactor.  In  this  way,  a  reactor- 
module  approach  to  combustion  modelling  could  be  preserved  as  a 
rational  framework  of  analysis.  Furthermore,  the  reactor-module 
approach  remains  pertinent  for  those  combustor  configurations  where  a 
spatial  description  is  not  readily  applicable  (see  Fig  8(c))  and  in 
which  controlling  shear  layer(s)  cannot  clearly  be  identified,  e.g., 
in  the  multi-inlet  dump  combustors  shown  in  Fig  2. 

In  the  next  section  the  discussion  of  flame  stabilization  will 


be  narrowed  to  consider  prior  work  which  is  relevant  to  dump  combustor 
configurations. 
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Flame  Stability  of  Dump  Combustors  and  of  Similar  Geometries 

As  noted  in  Chapter  I,  one  of  the  primary  incentives  for  this 
investigation  was  the  paucity  of  data  concerning  the  flame  stability 
of  coaxial  dump  combustors.  In  this  section,  the  meager  data  availa¬ 
ble  for  coaxial  combustors  will  be  discussed  first,  and  subsequently  a 
review  of  investigations  concerning  the  flame  stability  of  similar 
geometries  such  as  steps  and  wall  recesses  will  be  given. 

It  should  be  noted  that  coaxial  dump  combustors  have  been  utilized 
in  the  following  applications:  industrial  furnaces  (including  "tunnel" 
burners);  supersonic  combustion  ramjet  engines  (scramjets);  and  sub¬ 
sonic  combustion  IRR  engines.  Relevant  data  from  industrial  furnace 
studies  will  be  discussed  later.  The  coaxial  dump  combustor  has  found 
application  in  scramjet  engines  (Refs  67,  68,  and  79),  and  it  is  inter¬ 
esting  to  note  that  the  coaxial  configuration  is  one  that  can  be  used 
with  both  subsonic  and  supersonic  flows  at  the  dump  plane.  However, 
the  flow  and  combustion  phenomena  occurring  with  established  supersonic 
flows  differ  markedly  from  those  associated  with  subsonic  flows.  This 
current  discussion  will  be  restricted  to  subsonic  flows. 

Many  investigations  of  the  flame  stability  of  basic  coaxial 
combustors,  i.e.,  without  flameholders  at  the  entry,  have  been  reported. 
Initial  parametric  investigations  by  Barclay  (Ref  70)  of  a  number  of 
coaxial  combustor  geometries  with  liquid  fuel  (JP-4)  injection  yielded 
stability  data  with  no  clearly  discernable  correlating  parameter. 

Further  parametric  investigations  by  Stull  et  al .  (Ref  1)  using  liquid 
fuel  (JP-4)  injected  through  flush  wall  injectors  located  a  short  dis¬ 
tance  upstream  from  the  dump  plane  showed  that  the  fuel-air  ratios  at 
blow-off  could  be  correlated  on  the  basis  of  a  modified  loading  parameter 
given  by 
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where  =  pressure  in  atmospheres, 

=  temperature  (°R)/1000, 
h  =  step  height  =  (D-d)/2, 

D  =  combustor  diameter  (ft),  and 
d  =  inlet  diameter  at  the  dump  plane  (ft). 

This  data  is  shown  in  Fig  14.  It  is  significant  that  the  x 
parameter  apparently  collapses  data  from  geometrically  similar  combus¬ 
tors  ranging  in  diameter  from  2  to  5  inches.  Obviously,  this  x 
parameter  is  similar  to  the  form  of  correlating  parameter  Ua/P^D^7T^ 
used  to  describe  the  stability  of  bluff  body  flameholders ;  however, 
the  characteristic  dimension  of  the  flameholder  has  been  replaced  by 
the  step  height.  It  should  be  noted,  however,  that  such  correlations 
are  usually  applied  when  homogeneous  fuel-air  mixtures  are  used.  In 
this  case,  however,  the  correlation  appears  successful  even  with  wall 

injectors  close  to  the  dump  plane. 

The  investigation  by  Stull  et  al.  (Ref  1)  indicated  the  potential 
of  applying  bluff  body  correlating  parameters  to  coaxial  combustors. 
However,  a  basic  investigation  where  each  variable  was  independently 
varied  was  not  performed  in  this  series,  and  thus  the  functional  form 
of  the  correlation  could  not  be  confirmed.  Furthermore,  it  will  be 
observed  (Fig  14)  that  a  few  data  points  were  not  in  accord  with  the 
correlation. 

Results  of  further  investigations  performed  by  Craig  et  al.  (Ref  71) 
on  large-scale  dump  combustors  (12"  diameter)  are  shown  in  Fig  15. 

Once  again  fuel  was  injected  from  the  wall  a  short  distance  upstream 
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Figure  14.  Flame  Stability  Correlation:  Small  Coaxial 
Dump  Combustors 
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Figure  15.  Flame  Stability  Correlation:  12"  Coaxial  Dump  Combustor 


from  the  dump  plane.  In  this  case,  the  lean  extinction  limits  are 
similar  to  those  obtained  by  Stull  et  al.  (Refl).  However,  an  inter¬ 
esting  observation  was  that  decreasing  the  step  height  (h)  improved 
the  combustion  efficiency,  and  that  the  ratio  chamber  length  to  step 
height  (Lc/h)  was  more  significant  to  controlling  performance  than  the 
length  to  diameter  ratio  (Lc/D).  Indeed,  it  was  suggested  that  the 
optimum  step  size  would  be  the  smallest  one  capable  of  sustaining 
flame  over  the  operating  range  of  the  combustor.  A  further  investiga¬ 
tion  was  made  by  Chambers  (Ref  72)  of  two  dump  combustors  of  4"  and  5" 
internal  diameter  fed  by  a  2.5"  dia.  pipe.  In  this  case,  a  homogeneous 
premixed  fuel  (JP-4/air)  mixture  was  used.  Unfortunately,  only  30  data 
points  were  obtained  for  various  inlet  temperatures  (660°R  to  1248°R) 
and  pressures  (16  to  79  psi),  and  once  again  it  was  not  possible  to 
study  the  effect  of  varying  independently  the  inlet  pressure,  velocity, 
and  temperature.  The  reduced  data  (Ref  72)  is  displayed  in  Fig  16. 
Another  investigation  of  a  10"  diameter  coaxial  dump  combustor  fed  by 
a  3.35"  pipe  and  using  liquid  fuel  (JP-4)  flush-wall  injection  was 
reported  by  Schmotolocha  and  Economos  (Ref  73).  This  combustor  was  of 
relatively  low  d/D  ratio  and  had  an  overall  length  of  35",  which  was 
barely  long  enough  to  accommodate  the  recirculation  zone  generated  by 
the  step.  Data  from  this  investigation  is  shown  in  Fig  15. 

Another  recent  investigation  of  a  3"  diameter  coaxial  combustor 
with  a  2"  diameter  inlet  has  been  reported  by  Choudhury  (Ref  74),  and 
the  corresponding  lean  limits  are  shown  in  Fig  17.  This  set  of  data 
was  apparently  obtained  at  one  inlet  temperature  (590°R)  and  at  approxi¬ 
mately  atmospheric  pressure.  An  interesting  aspect  of  this  investigation 
was  the  use  of  different  wall  materials  which  provide  an  indication  of 


Figure  16.  Flame  Stability  Correlation  for  Dump  Combustors:  Pre-Mixed  Fuel-Air 


Stability  Oata  for  a  Coaxial  Dump  Combustor  and  for  Steps 


the  effect  of  wall  heat  losses  on  flame  stability.  In  this  case, 
pre-mixed  propane  was  used  as  a  fuel. 

All  the  above  investigations  were  performed  in  support  of  ramjet 
combustor  technology.  An  overall  observation  concerning  the  above  data 
is  that  in  no  case  is  the  existence  of  a  complete  flame  stability  loop 
unequivocably  confirmed  by  the  data,  although  some  investigators  have 
sketched  in  the  anticipated  curves.  It  is  also  interesting  to  note 
that  the  data  of  Choudhury  (Ref  74)  shows  only  a  slight  rise  in  the 
lean  extinction  fuel-air  ratio  as  the  velocity  increases.  Similarly, 
inspection  of  Figs  14  to  17  shows  a  relatively  constant  lean  limit  to 
exist  over  wide  ranges  of  operation.  It  will  also  be  noted  that  the 
majority  of  work  has  been  performed  with  heterogeneous  fuel -air  mix¬ 
tures,  and  little  data  have  been  acquired  with  homogeneous  mixtures. 

In  addition  to  the  above  combustor  studies,  two  basic  investigations 
are  noteworthy.  The  first  effort  by  Ross  (Ref  75)  was  a  study  of  the 
flame  stability  of  a  small  scale  (D  =  1",  d  =  0.5",  approximately) 
coaxial  combustor  fabricated  from  Vycor  glass  tubing.  This  investigation 
was  conducted  to  study  the  effects  of  oscillatory  combustion  on  the 
limits  of  stabilization.  A  variable  length  combustor  (16"-20")  was 
employed  to  vary  the  oscillation  frequency.  The  fuel  used  was  propane, 
and  the  data  obtained  are  shown  in  Fig  18.  Although  the  long  test  con¬ 
figurations  (Lc  =  16",  18",  and  20")  were  not  representative  of 

ramjet  combustor  practice,  a  number  of  relevant  observations  can  be 
made.  Because  of  the  transparent  combustor,  it  was  possible  to 
observe  a  distinct  annular  recirculation  zone  with  attached  flame 
brush.  Also,  as  the  lean  limit  was  approached,  the  flame  retreated 
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Figure  18.  Flame  Stability  Data  for  a  Small  Coaxial  Dump  Combustor 


towards  the  dump  plane,  appearing  as  a  faint  annular  brush  just  prior 
to  blow-out.  For  the  16"  and  18"  chambers,  the  lean  limit  equivalence 
ratio  remained  relatively  constant  over  a  wide  range  of  velocity 
(90-200  feet/sec),  but  for  the  20"  chamber  the  weak  limit  value 
increased  progressively  with  velocity  (Fig  13).  As  noted  earlier,  a 
longitudinal  resonant  oscillation  was  encountered  in  this  investigation, 
the  occurrence  of  which  completely  modified  the  flame  pattern.  This 
was  perhaps  the  first  warning  of  similar  oscillatory  phenomena  which 
have  been  manifested  in  recent  dump  burner  research  (Ref  1). 

The  other  basic  research  effort  is  that  of  Baev  and  Tret'yakov 
(Ref  58).  In  this  study,  a  symmetric, two-dimensional  dump  combustor 
with  transparent  side  walls  was  utilized.  Two  step  heights  (25  and 
17  mm)  were  studied,  and  both  rich  and  weak  limits  were  determined  as 
a  function  of  velocity.  Although  a  closed  loop  was  not  obtained, 
reduced  step  height  produced  a  marked  deterioration  in  both  weak  and 
rich  limits.  The  length  of  the  recirculation  zones  behind  the  steps 
was  determined  from  flame  photographs.  The  zone  length  was  not 
affected  by  fuel-air  ratio,  but  increased  markedly  with  combustion 
chamber  exit  velocity;  however,  no  exhaust  nozzle  appeared  to  be  used. 
Unfortunately,  the  reporting  of  this  important  investigation  was 
obscured  by  its  brevity. 

For  completeness,  mention  is  made  of  an  investigation  of  a  simple 
coaxial  combustor  reported  by  Altenkirch  and  Mel  1  or  (Ref  66).  In 
this  case,  the  flame  stabilizing  element  was  formed  by  a  simple  orifice 
in  a  circular  plate  spannvn  the  combustor.  The  combustor  geometry  is 
thus  similar  to  that  of  a  coaxial  dump  combustor.  Unfortunately,  the 
stability  data  reported  cannot  be  readily  expressed  in  terms  of 


blow-out  velocity.  So  far,  attention  has  been  confined  to  combustion 
systems  of  the  sudden-expansion  type.  Several  flame-stability  investi¬ 
gations  have  been  performed  on  related  geometries  such  as  two-dimension 
steps  and  wall  recesses.  These  tests  will  now  be  reviewed. 

A  most  recent  study  is  that  of  Kawamura  (Ref  76).  In  this 
investigation,  propane  fuel  was  ejected  as  a  thin  sheet  from  a  slot  at 
the  top  of  backward-facing  step.  The  fuel,  mixed  with  the  co-flowing 
airstream,  and  the  behavior  of  the  resulting  ignition  front  was  observed. 
Although  this  work  is  not  truly  comparable  to  one  with  homogeneous 
fuel-air  mixture,  two  interesting  observations  can  be  made.  The  first 
is  that  the  lean  and  rich  blow-off  limits  remained  apart,  and  a  closed 
stability  loop  was  never  obtained.  Second,  the  length  of  the  RZ  under 
burning  conditions  was  about  half  that  of  the  comparable  cold  flow 
value  (3.5  step  heights  compared  to  7  under  cold  flow  conditions). 
However,  the  RZ  length  exhibited  little  variation  with  fuel/air  ratio 
or  freestream  velocity. 

An  earlier  investigation  by  Snyder  (Ref  77)  concerned  the  flow  of 
pre-mixed  propane-air  over  a  step.  In  this  case,  the  boundary  layer 
thickness  at  the  dump  position  was  systematically  varied  by  using 
different  approach  lengths  upstream  of  the  step.  Unfortunately,  the 
step  height  (h)  was  not  reported,  and  it  is  thus  not  possible  to 
assess  the  relative  effect  of  boundary  layer  thickness.  A  general 
observation  was  that  increasing  boundary  layer  thickness  improved 
stability.  A  complete  stability  loop  was  obtained,  but  the  incomplete 
reporting  of  test  data  leaves  little  room  for  further  comment.  However, 
once  again  it  was  noted  that  the  RZ  length  was  found  to  remain  essen¬ 
tially  constant  above  a  Reynolds  number  of  2.5  x  10^. 
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For  completeness,  another  step-stabilization  investigation  should 
be  noted.  This  investigation  was  performed  by  DeZubay  (Ref  78),  with 
the  primary  goal  of  investigating  the  effect  of  wall  heat  transfer  with 
the  recirculating  zone;  some  data  on  rich  limits  are  given,  but  not 
introduced  here.  However,  it  should  be  noted  that  heat  loss  from  the 
RZ  is  more  likely  to  occur  with  dump  or  step  configurations  compared 
to  bluff  body  stabilizers. 

Another  configuration  related  to  the  dump  combustor  is  the  flame 
stabilizer  configuration  formed  by  machining  a  recess  in  the  chamber 
wall.  This  "recessed-wall"  configuration  is  attractive  for  practical 
combustors  (Ref  79  and  80)  because  it  offers  flame  stabilization  with 
no  mechanical  blockage  of  the  stream  and  thus  with  low  pressure  loss. 
However,  the  floor  of  the  recess  is  subject  to  high  heat  loads  due  to 
the  recirculating  flow  in  the  cavity. 

A  pioneering  investigation  of  recessed  wall  flameholders  was 
reported  by  Huellmantel  et  al.  (Ref  81).  In  this  case,  a  pre-mixed 
propane  air  stream  was  used,  and  the  test  section  was  a  symmetrical 
rectangular  channel  with  recesses  cut  on  the  top  and  bottom  walls  of 
the  channel.  To  determine  the  effect  of  recess  shape,  12  different 
recess  contours  were  tested.  This  study  found  that  recesses  of  short 
length  provided  inadequate  RZ  volume  for  stability  and  were  associated 
with  rough  burning.  The  influence  of  recess  depth  (i.e.,  effectively 
the  step  height)  was  also  examined;  reducing  the  step  height  from  the 
maximum  of  0.25  in.  resulted  in  reduced  stability  limits,  as  would  be 
anticipated. 

A  further  extension  of  this  initial  work  was  made  by  Povinelli 
(Ref  82),  the  main  achievement  of  which  was  to  correlate  the  data  of 
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Huellmantel  et  al.  on  the  basis  of  the  Zukoski-Marble  criterion.  An 
interesting  observation  in  this  case  was  that  the  length  of  the  RZ  was 
essentially  constant  despite  variations  in  flow  velocity  and  equivalence 
ratio. 

A  recent  and  comprehensive  investigation  of  recess  wall/step 
configurations  has  been  carried  out  by  Choudhury  (Ref  74).  Although 
the  main  thrust  of  this  work  was  the  enhancement  of  stability  and 
combustion  efficiency  by  the  use  of  a  system  of  small  air  jets  injected 
upstream  of  the  dump  plane,  much  of  the  experimental  work  covered  basic 
recessed  wall/step  geometries.  As  already  shown  in  Fig  17,  the  blow-off 
limits  for  both  two-dimensional  steps  and  an  axisymmetric  combustor  show 
little  change  over  the  range  of  velocity  explored.  The  effects  of  vary¬ 
ing  inlet  temperature  and  pressure  were  not  studied.  The  effect  of 
step  height  was,  however,  examined;  basically  no  effect  of  step  height 
on  stability  was  found  for  heights  between  0.25  and  0.75  inch.  However, 
reduction  in  step  height  below  0.25  in.  resulted  in  degradation  of 
stability.  The  length  of  the  RZ  was  found  to  remain  between  6  and  7 
step  heights  for  cold  flow  conditions  and  "similar  behavior"  was 
observed  during  burning.  A  major  service  performed  by  Choudhury  has 
been  to  draw  increased  attention  to  the  pressure  fluctuations  and  noisy 
operation  associated  with  dump  combustors.  Even  in  cold  flow  the 
reattachment  process  associated  with  the  free  shear  layer  gave  rise  to 
marked  pressure  fluctuations  and,  of  course,  the  related  phenomena  of 
"cavity  tones"  is  well  known.  However,  the  introduction  of  combustion 
into  this  environment  can  lead  to  severe  manifestations  of  rough 
burning.  Choudhury  has  documented  this  phenomenon  in  some  detail,  and 
no  further  discussion  of  this  combustion  instability  will  be  given. 
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The  limited  results  available  concerning  the  flame  stabilizing 
properties  of  steps/recesses  indicate  that  the  overall  RZ  length 
( L/h )  may  change  from  cold  flow  to  burning  conditions.  However,  no 
marked  variation  of  L/h  occurs  with  varying  flow  velocity  and  fuel- 
air  ratio  (Refs  76,  77,  and  82).  The  results  concerning  the  effect  of 
step  height  variation  on  stability  may  be  explained  by  postulating  a 
critical  step  height  of  about  0.25  inch.  Above  this  critical  height, 
the  step  dimension  does  not  appear  to  affect  the  stabilizing  limits 
(Ref  74).  Below  this  height,  the  stability  range  is  narrowed  by 
decreasing  step  height  (Refs  74,  81).  However,  the  results  presented 
earlier  (Ref  58)  concerning  the  opposing-step  combustor  configuration 
of  Baev  and  Tret'yakov  do  not  harmonize  with  the  above  observations 
concerning  the  effect  of  velocity  on  the  RZ  length  or  the  existence  of 
a  critical  step  height. 

The  limited  experimental  work  concerning  the  flame  stability  of 
coaxial  dump  combustors  and  related  geometries  has  been  reviewed.  It 
is  now  appropriate  to  turn  to  the  theoretical  analysis  of  the  dump 
combustor  with  the  object  of  modelling  the  flame  stabilizing  process. 
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Ill  THEORETICAL  ANALYSIS 


In  Chapter  II,  it  was  pointed  out  that  a  combustion  chamber  can  be 
modelled  as  a  network  of  reactor  modules.  To  preserve  this  framework 
of  analysis,  the  lean  extinction  limits  of  a  coaxial  dump  combustor 
will  be  estimated  by  assuming  the  recirculation  zone  to  be  a  stirred 
reactor.  In  order  to  evaluate  the  air  loading  expression,  Eq(10),  of 
the  stirred  reactor,  two  terms  are  needed  which  are  determined  by  the 
aerodynamics  of  the  combustor.  These  terms  are  the  recirculation  zone 
volume  (V)  and  the  mass  flow  (mr)  entrained  into  the  zone.  An  exten¬ 
sive  review  of  experimental  data  concerning  flow  phenomena  in  coaxial 
dump  configurations  was  therefore  undertaken  and  is  presented  in 
Appendix  A.  The  results  of  this  review  will  now  be  briefly  outlined. 
Aerodynamics  of  Coaxial  Dump  Configurations 

The  general  features  of  the  flow  field  under  non-burning  conditions 
are  shown  in  Fig  19.  Separation  of  the  flow  occurs  at  the  inlet  dump 
plane  and  is  followed  by  reattachment  to  the  wall  further  downstream. 

The  reverse  flow  velocity  in  the  annular  recirculation  zone  is  also 
depicted.  The  reattachment  process  at  the  wall  is  a  very  unsteady 
phenomenon.  Thus,  the  point  of  reattachment  was  observed  by 
DeRossett  and  Przirembel  (Ref  83)  to  oscillate  at  a  frequency  of  the 
order  of  several  hundred  Hertz.  Thus,  reattachment  occurs  over  a 
region  rather  than  on  a  line.  Furthermore,  the  width  of  this  region 
has  been  found  (Ref  83)  to  broaden  with  both  the  step  height  and  the 
initial  thickness  of  the  inlet  boundary  layer;  appreciable  pressure 
fluctuations  also  accompany  the  reattachment  process.  The  recircula¬ 
tion  zone  was  also  found  (Ref  83)  to  include  a  sub-zone  just  downstream 
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of  the  step,  and  this  is  also  shown  in  Fig  19.  Also,  several 
investigators  (Refs  84,  85,  and  86)  have  shown  evidence  that  a 
rotary  motion  about  the  chamber  axis  occurs  within  the  recirculation 
zone.  Despite  the  complexity  of  the  reattachment  process,  the  time- 
mean  length  (L)  of  the  recirculation  zone  may  be  determined  within 
reasonable  limits.  In  Fig  20,  a  composite  plot  of  the  values  of  L/h 
obtained  by  various  investigators  is  given.  It  will  be  noted  that  at 
high  Reynolds  numbers,  the  ratio  L/h  appears  to  become  relatively 
constant  and  independent  of  the  flow  Reynolds  number.  Typically,  the 
length  (L)  is  found  to  be  between  8  and  9  step  heights  (Rets  84  and  85). 
However,  both  Pennucci  (Ref  84)  and  Drewry  (Ref  85),  who  tested  coaxial 


chambers  fitted  with  terminal  convergent  nozzles,  observed  that  a  reduc¬ 
tion  in  nozzle  throat  size  appeared  to  reduce  the  reattachment  length 
(L)  slightly.  Limited  experimental  data  (Refs  87  and  88)  indicate 
that  L/h  may  also  be  affected  by  the  diameter  ratio  (d/D).  However, 
for  modelling  purposes,  one  may  assume  that  k  =  L/h  is  approximately 
constant  for  Reynolds  numbers  greater  than  about  3  x  105. 

The  boundary  of  the  recirculation  zone  determined  in  experimental 
studies  may  be  approximated  by  a  parabolic  curve  as  shown  by  the  data 
of  Chaturvedi  (Ref  89).  Thus,  with  the  length  and  shape  of  the 
recirculation  zone  known,  the  corresponding  volume  can  be  calculated. 

In  addition  to  the  recirculation  zone  volume,  the  mass  flow 
entering  into  the  recirculation  process  is  required  in  order  to  cal¬ 
culate  the  air  loading  of  the  reaction  zone.  Typically,  in  analyses 
of  bluff  body  flame  stabilizing  processes,  this  entrained  mass  flow  is 
usually  calculated  (Ref  21  and  31)  by  an  expression  of  the  form 

mr  a  pSU 
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where  S  is  an  effective  entrainment  surface  area,  U  is  the  local 
flow  velocity  at  the  trailing  edge  of  the  bluff  body,  and  P  is  an  appro¬ 
priate  density.* 

2 

Usually  S  is  assumed  proportional  to  ,  yielding  the  result  that 

mr  a  P  Df2  U  (48) 

However,  in  the  case  of  the  coaxial  dump  geometries,  a  number  of 
investigations  has  been  made,  in  connection  with  industrial  furnace 
development,  to  measure  the  maximum  fraction  of  the  inlet  mass  flow 
entering  into  the  recirculation  process.  This  extensive  literature  is 
also  reviewed  in  Appendix  A. 

Typically,  industrial  furnaces  are  fed  by  two  streams  consisting 
of  a  primary  jet  which  carries  the  fuel,  surrounded  by  a  secondary  flow 
of  air.  With  zero  secondary  air  flow,  the  furnace  operation  corresponds 
to  that  of  a  coaxial  dump  combustor.  This  case  is  documented  in  the 
literature  (Refs  90,  91,  and  92).  Unfortunately,  the  diameter  ratio 
(d/D)  of  a  typical  furnace  is  usually  much  smaller  than  that  of  a  ram¬ 
jet  coaxial  combustor.  However,  some  useful  data  have  been  obtained 
from  furnace  investigations.  The  results  of  various  furnace  studies  are 
given  in  Fig  21,  showing  recirculated  mass  flow  fraction  (mrc/m0)  as  a 
function  of  diameter  ratio  (d/D).  The  data  are  well  represented  by  the 
following  equation  due  to  Wingfield  (Ref  91). 

(49) 

where  m„  denotes  the  primary  jet  mass  flow  and  iti  is  the  measured 
o  rc 

recirculation  mass  flow.  At  this  point,  one  may  now  proceed  to 
evaluate  the  air  loading  parameter  of  the  recirculation  zone. 

*  See  discussion  in  Appendix  B. 
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Recirculated  Mass  Flow  Fraction  -  Diameter  Ratio 
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Loading  Parameter  of  Recirculation  Zone 


The  loading  parameter  of  interest  is  that  given  in  Eq(10),  and  for 


the  moment,  attention  is  drawn  to  the  quantity 


fen)  '  ™ 


s  may  be 


expressed  in  terms  of  a  mass  flow  fraction  (mr/mn)  as  follows: 


For  the  coaxial  dump  combustor,  one  may  set  Pr  equal  to  P0  with 
negligible  error  to  obtain 


Thus  for  given  inlet  conditions,  the  terms  (ror/mQ)  and  (A,o/V)  determine 
reactor  loading.  Alternatively,  the  quantity  (mr/VPrn)  may  be  evaluated 
using  the  conventional  expression  for  the  entrained  mass  flow  given  in 
Eq(47).  In  this  case 


where  kt  represents  an  entrainment  factor,  and  P  has  been  put 
equal  to  pQ. 

From  Eq(50)  and  Eq ( 52 )  one  has  the  equality 


In  this  case,  an  estimate  of  the  consumed  mass  fraction  (mr/m0)  may  be 
derived  from  Eo(53)  if  the  geometry  of  the  recirculation  zone  and  the 
effective  value  of  kt  are  known.  As  a  crude  approximation,  one  may 
assume  that  the  zone  is  bounded  by  a  straight  line  extending  from  the 
point  of  separation  at  the  dump  plane  to  the  point  of  reattachment  at 
the  chamber  wall.  It  is  further  assumed  that  the  point  of  reattachment 
is  located  at  a  distance  kh  from  the  dump  plane.  As  noted  previously, 
k  is  approximately  equal  to  8.5  for  cold  flows.  The  surface  area  of 
the  recirculation  zone  can  be  obtained  by  integration,  and  is  given  by 


=  rr  k  h  (D  -  h) 


(54) 


where  denotes  the  surface  area  based  on  a  straight  line  boundary. 
The  corresponding  volume  is 


Vi 


£-jp-  (3D  -  2h) 


(55) 


For  k  >  >  1,  Eq(54)  may  be  approximated  as 


S ^  =  n  kh  (D  -  h) 


(56) 


and  the  corresponding  recirculated  mass  flow  fraction  becomes 
m 


t—  =  k  t  k 

m0 


1 


smi 


1 


(57) 


which  for  high  diameter  ratios  (d/D),  and  with  appropriate  choice  of 
coefficients,  yields  results  similar  in  form  to  those  given  by  Eq(49). 
A  characteri sti c  dimension  may  also  be  derived  for  this  linear  RZ 
geometry  by  writing 

t-  v/  D  /2h\  ^3"~d) 

~  S/  6  \D  /  (z-^) 


(58) 
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Correspondingly,  Eq ( 52 )  may  be  expressed  as 


VPr" 

and  this  form  of  equation  is  relatively  insensitive  to  variations  in  the 
value  of  k.  Also  h  is  effectively  proportional  to 

However,  as  noted  earlier,  a  better  approximation  to  the  recirculation 
zone  shape  is  obtained  with  a  parabolic  curve.  In  this  case,  the  expres¬ 
sion  for  the  surface  area  (S  )  becomes  * 
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The  corresponding  volume  is  given  by 

2 


V_ 


ir  kh 


2  ,  ,  4  . 
3  d  +  5  h 


(61) 


Once  again  a  characteristic  dimension  (hp)  may  be  defined  expressed  as 


h  i  V 

0  D  Sp 


(62' 


The  ratio  hp/D  may  be  calculated  from  Eq ( 60 )  and  Eq (61),  using  the 
identity  D  =  d  +  2h.  A  curve  of  hp/D  versus  2h/D,  for  various 
values  of  k,  is  given  in  Fig  22.  For  k  =  8.5,  a  good  representa¬ 
tion  cf  this  relationship  is 

1.095 


•(£) 
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where  0.1  <  <  0.6. 


*  See  Appendix  A,  Eq(A-l)  and  Eq(A-2). 
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If  the  exponent  in  Eq(63)  were  equal  to  unity,  the  diameter  D  could 

be  eliminated  and  Fp  would  be  exactly  proportional  to  h. 

In  addition  to  the  theoretical  estimates  and  measured  values  of 

recirculated  mass  flow  fraction,  it  is  possible  in  principle  to  measure 

the  amount  of  air  entrained  into  the  recirculation  zone  and  consumed  under 

actual  operating  conditions  by  a  technique  introduced  by  Lefebvre  et  al. 

(Ref  93).  This  technique  was  used  to  study  the  air  entrainment  into 

the  recirculation  zone  of  a  bluff  body  and  consists  of  injecting  a  small 

fraction  of  the  total  fuel  flow  directly  into  the  recirculation  zone. 

•  • 

Both  the  mainstream  fuel  flow  (fQ)  and  the  fuel  flow  (f  )  injected  into 

the  recirculation  zone  are  measured  at  the  lean  extinction  condition. 

« 

The  total  fuel  flow  (ft)  entering  the  recirculation  zone  is  then 


where  mr  is  the  air  mass  flow  burned  in  RZ,  and  m„  denotes  the 
total  air  mass  flow  at  the  dump  plane. 

The  assumption  is  made  that  the  lean  extinction  occurs  at  the  same 
overall  fuel -air  ratio  in  the  recirculation  zone  irrespective  of  whether 
the  fuel  enters  the  zone  by  direct  injection  or  by  entrainment.  Thus, 
for  fixed  upstream  conditions,  the  total  fuel  flow  (ft)  should  remain 
constant  irrespective  of  the  division  of  fuel  between  the  mainstream 
and  the  recirculation  zone.  The  following  relationship  then  exists 
between  the  various  fuel  flows 
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Consequently,  for  given  upstream  conditions,  a  plot  of  f  against 
f  should  give  a  straight-line  relationship,  the  slope  of  which  yields 
the  desired  quantity  (mr/m0). 

This  technique  has  been  successfully  applied  to  a  gas  turbine 
combustor  configuration  by  Verduzio  and  Campanaro  (Ref  94). 

Chemical  Reaction  Model 

In  order  to  estimate  the  lean  stability  limit  of  the  coaxial  dump 

combustor,  the  annular  recirculation  zone  is  modelled  as  a  stirred 

reactor.  For  the  combustion  of  a  hydrocarbon  denoted  by  C  H  with 

x  y 

excess  air,  an  expression  for  the  air  loading  term  0  was  given  in 
Eq(5).  This  equation  was,  however,  for  a  reaction  of  order  two.  The 
more  general  expression  for  a  reaction  of  order  n  is  derived  in 
Appendix  C  and  is  given  in  the  following  equation: 

VPrn  Rn  Tn“0'5  €<t>  |(|TH-I)(x  +  y/4)  +<P  +  (y/4-l)Jn 

(66) 


It  is  assumed  that  an  effective  activation  energy  (E)  and  a  reaction 
order  (n)  can  be  assigned  for  the  combustion  reaction.  Then  for  a  given 
equivalence  ratio  ($)  and  initial  air  temperature  (T),  the  air  loading 
achieves  a  maximum  value  (#max)  at  a  certain  value  of  e.  This  maximum 
value  0max  is  the  highest  loading  the  reaction  zone  can  sustain  before 
extinction.  Values  of  t—  0mav  have  been  calculated  for  both  ethylene 

|T!<iX 

and  JP-4  fuels  burning  in  air,  for  inlet  temperatures  varying  from  600- 
1200°R,  and  for  equivalence  ratios  ranging  from  0.3  to  1.0.  These 
calculations  were  made  for  values  of  E  varying  from  30  to  54  KCal/Gm 
Mole,  and  for  values  of  n  from  1.25  to  2.00. 
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can 


Ballal  and  Lefebvre  (Ref  31)  have  shown  that  the  term  0 

IIWa 

be  approximately  represented  as  a  function  of  equivalence  ratio  and 
initial  temperature,  of  the  form 


where  b  and  t  are  constant  for  an  assumed  activation  energy  and 
reaction  order.  Alternatively,  one  may  express  the  variation  of 

^Wx  as 

*max  *  TS*‘  (68) 

Calculated  values  of  s  and  t  are  shown  in  Fig  23  for  ethylene; 
similar  data  for  JP-4  is  given  in  Appendix  C.  For  ethylene,  an  activa¬ 
tion  energy  of  about  42  KCal/Gm  Mole  would  be  anticipated.  If  a 
reaction  order  of  1.75  is  assumed,  then  the  corresponding  value  of  s 
is  2.86  and  the  value  of  t  is  5.3. 

From  Eq(68),  the  maximum  loading  may  be  expressed  as 

%  C  f 

-~r  »  h  T  *  (69) 

VPrn 

In  practice  the  actual  loading  at  which  a  combustor  blows  out  is 
significantly  less  than  the  theoretical  maximum  value  #max»  princi¬ 
pally  because  the  combustor  is  not  perfectly  mixed.  This  proportionate 
reduction  in  actual  loading  compared  to  the  ideal  loading  is  included 
in  the  constant  term  k3. 

From  Eq(59),  with  PQ  ~  Pr,  one  may  write 
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E  (KCal/Mole) 


and  t  as  a  Function  of  E  and 


Finally,  using  Eq(69),  and  assuming  h £<*  h,  one  obtains  the  following 
expression  for  the  lean  blow-out  limit 

1 

*uso  *  0  •  JTTn^TT  •  (70) 

Thus  for  the  previous  assumptions  made  for  ethylene  (E  s  42  KCal/Gm 
Mole,  n  =  1.75)  yielding  s  =  2.86,  t  *  5.3,  one  obtains 


<t> 
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(71) 


Thus  a  correlating  parameter  is  anticipated  of  the  following  form: 


^LBO  =  k4  p  b  (72) 

o  o 

In  practical  applications,  the  use  of  the  stagnation  temperature  T* 

Lo 

is  usually  preferred  rather  than  TQ.  For  modest  flow  velocities, 

Tt0  ~  To  and  direct  substitution  of  Tt0  in  Eq(72)  is  permissible. 
For  those  higher  inlet  velocities  where  compressibility  effects  are 
important,  Eq(72)  may  be  expressed  in  terms  of  M0  and  Tt0  rather 
than  U0  and  T0.  However,  the  use  of  U0  and  Tt0  is  a  preferred 
form,  and  the  small  errors  involved  in  assuming  Tt0  ~  T0  are  usually 
accepted. 

It  should  be  emphasized  that  the  exponents  given  in  Eq (71 )  are 
obtained  on  the  basis  of  a  one-step  reaction  model  described  by  a 
given  activation  energy  and  reaction  order.  A  perfectly  stirred 
adiabatic  reaction  zone  was  also  assumed.  Consequently,  the  values  of 
these  exponents  derived  from  practical  combustors  which  are  neither 
perfectly  stirred  nor  completely  adiabatic  may  deviate  substantially 
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from  the  theoretical  values.  Furthermore,  the  common  assumption  made 
in  Eq (52)  that  P  *  P0  will  also  contribute  to  such  deviations. 

It  is  now  appropriate  to  consider  the  experimental  program  planned 
both  to  explore  the  validity  of  the  form  of  correlation  given  in 
Eq(72)  and  to  estimate  the  value  of  the  recirculated  mass  flow  fraction. 


IV  EXPERIMENTAL  PROGRAM 


Aims  of  the  Experimental  Program 

The  first  goal  of  the  experimental  program  is  to  determine  the 
lean  extinction  limits  of  the  stabilizer  geometry  created  by  the  sudden 
expansion  section  of  a  coaxial  dump  combustor.  In  particular,  the 
effects  of  inlet  pressure,  temperature,  velocity,  and  geometry  on  the 
lean  extinction  limit  are  sought.  Ideally,  the  effect  of  combustor 
scale  would  also  be  included  in  this  investigation;  however,  in  order 
to  keep  the  scope  within  available  resources,  this  investigation  is 
restricted  to  a  6"  diameter  combustor.  Geometry  variation  is  accom¬ 
plished  by  varying  the  diameter  of  the  inlet  from  2,5"  to  5.0"  with 
corresponding  step  heights  varying  from  1.75"  to  0.5".  The  second 
goal  of  the  experimental  program  is  to  evaluate  the  utility  of  the 
form  of  correlating  parameter  given  in  Eq(72)  over  the  range  of  experi¬ 
mental  conditions  and  geometries  available. 

The  third  goal  is  to  determine  experimentally  values  of  the 
recirculated  mass  flow  fraction  by  employing  the  technique  of  split- 
fuel  injection  introduced  by  Lefebvre  et  al.  (Ref  93). 

Before  describing  the  design  of  the  experimental  hardware  and 
execution  of  the  experimental  program,  it  is  appropriate  to  highlight 
some  of  the  problems  peculiar  to  determining  flame  stability  limits. 
Determination  of  Stability  Limits 

The  experimental  determination  of  accurate  blow-off  limits  for  a 
dump  combustor  is  not  a  simple  task  because  of  the  complex  flow  inter¬ 
actions  which  occur  between  the  chemical,  aerodynamic,  and  acoustic 
phenomena  associated  with  combustion.  A  practical  dump  combustor 
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experiences  highly  turbulent  conditions  due  to  the  nature  of  the 
separation/reattachment  process  so  that  perturbations  continually 
exist  in  local  flow  conditions.  As  the  lean  blow-out  limit  (LBO)  is 
approached,  the  combustion  process  becomes  increasingly  vulnerable  to 
disturbances,  and  successive  determinations  of  the  blow-off  fuel-air 
ratios  will  exhibit  substantial  scatter,  even  under  apparently  similar 
operating  conditions. 

Detection  of  the  actual  point  of  blow-off  is  another  source  of 
difficulties.  A  typical  experimental  procedure  to  determine  a  lean 
blow-off  point  is  to  gradually  reduce  the  fuel  flow  until  extinction 
occurs.  Unfortunately,  reduction  in  fuel  flow  is  usually  accompanied 
by  a  corresponding  increase  in  air  flow,  and  the  fuel  and  air  flows  at 
the  exact  instant  of  blow-out  can  be  difficult  to  determine  with  pre¬ 
cision.  Additionally,  the  detection  of  the  point  of  blow-out  can  be 
difficult  if  the  cessation  of  combustion  cannot  be  visually  observed 
or  aurally  detected.  In  many  cases,  a  thermocouple  can  be  placed  in 
the  recirculation  zone  to  signal  blow-off,  but  the  resulting  lags  in 
detection  can  be  appreciable.  Additionally,  the  rate  at  which  the 
fuel  flow  is  changed  as  the  blow-off  limit  is  approached  can  signifi¬ 
cantly  affect  the  result;  this  is  not  a  trivial  matter  because  in  many 
cases  blow-off  points  have  to  be  determined  expeditiously  to  avoid 
overheating  of  the  combustor. 

In  the  case  of  isolated  bluff  body  stabilizers  or  simple 
derivatives  of  such  systems,  testing  can  be  accomplished  under  non- 
constrained  (constant  pressure)  or  constrained  (constant  area) 
conditions.  In  non-constrained  operation,  testing  can  be  carried 
out  in  a  "free-jet"  mode,  enabling  the  velocity  and  temperature  at  the 


82 


flameholder  to  be  varied  independently,  while  the  pressure  level 
remains  constant.  In  constrained  operation,  the  presence  of  a  choked 
exhaust  nozzle  constrains  the  inlet  velocity  to  vary  in  a  manner  which 
is  related  to  the  effective  heat  release  in  the  combustion  chamber. 
Typically,  in  most  experimental  tests  the  combustor  may  be  operated 
either  at  chosen  inlet  pressure  and  temperature  conditions,  with  the 
inlet  velocity  varying  with  fuel -air  ratio,  or  at  constant  mass  flow 
and  temperature  conditions  with  the  inlet  pressure  and  velocity  chang¬ 
ing  with  fuel -air  ratio.  If  the  flameholder  is  subjected  to  identical 
inlet  conditions  under  both  "constrained"  or  "non-constrained"  opera¬ 
tion,  then  the  corresponding  blow-off  fuel-air  ratios  should  be  the 
same.  In  practice,  significant  and  unpredictable  differences  can  exist 
due  to  the  dynamic  coupling  effects  which  occur  under  constrained 
operation.  Generally,  the  unconstrained  stable  operating  range 
decreases  when  a  choked  exhaust  nozzle  terminates  the  chamber,  and  also 
as  the  combustion  chamber  length  is  increased.  These  effects  could  be 
classified  as  changes  in  the  downstream  boundary  conditions  which 
affect  the  behavior  of  flow- induced  oscillations  and  thereby  the 
blow-off  condition. 

An  unexpected  bonus  from  the  use  of  a  choked  exhaust  nozzle,  however, 
is  that  the  lean  blow-off  is  usually  attained  crisply  and  can  be  deter¬ 
mined  with  increased  precision.  Furthermore,  the  test  chamber  is 
effectively  decoupled  from  the  vagaries  of  the  facility  exhaust  system. 
Another  source  of  dynamic  interaction  is  the  inlet  piping  system 
upstream  of  the  combustor.  Pressure  pulses  arising  from  the  combustion 
process  can  cause  significant  changes  in  inlet  flow  conditions,  particu¬ 
larly  if  the  inlet  duct  is  operating  near  an  Incipient  separation 
condition. 


The  above  comnents  are  made  in  the  context  of  determining  lean 
blow-out  limits.  The  determination  of  rich  blow-out  limits  is  even 
more  problematical  because  of  the  generally  violent,  and  sometimes  near 
explosive,  behavior  of  combustors  as  the  rich  limit  is  approached. 
Indeed,  the  discharge  of  fuel-rich  products  into  facility  exhaust  sys¬ 
tems  can  be  extremely  hazardous,  and  for  this  reason  such  determinations 
are  often  prohibited. 

From  the  above  discussion,  it  will  be  appreciated  that  the 
determination  of  accurate  flame  stability  limits  is  not  entirely 
straightforward.  However,  because  of  the  ability  to  test  bluff  body 
stabilizers  in  an  unconstrained  environment,  several  in-depth  studies 
or  recirculation  zone  processes  have  been  made  (Refs  49,  56,  and  57). 
Such  stabilizers  can  be  examined  with  comparative  ease,  using  either 
optical  or  probe  instruments. 

In  contradistinction  to  the  bluff  body  case,  the  coaxial  dump 
combustor  presents  a  much  more  challenging  task.  Firstly,  the  sta¬ 
bilizing  element  cannot  be  decoupled  from  the  combustor  chamber  proper 
because  it  is  inherently  part  of  the  chamber  geometry.  Secondly,  the 
hot, recirculating  gases  directly  scour  the  chamber  walls,  thus  making 
it  difficult  to  use  transparent  materials  and  usually  restricting 
operating  time  due  to  excessive  wall  heating.  It  is  also  very  diffi¬ 
cult  to  optically  probe  just  one  segment  of  the  recirculation  zone 
since  the  line  of  sight  will  generally  pass  through  two  zones. 

Thirdly,  it  is  undesirable  to  use  physical  probes  in  this  particular 
environment  because  the  presence  of  a  relatively  large,  cooled  probe 
leads  to  significant  disturbance  of  a  recirculating  flow  field. 

Finally,  the  basic  nature  of  the  separation-reattachment  process 
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yields  a  highly  turbulent  quasi-steady  flow  field  which,  given  the 
needed  instrumentation,  would  require  extensive  and  time-consuming 
measurements  for  detailed  characterization. 


The  foregoing  considerations  have  been  taken  into  account  in  the 
design  of  the  experimental  combustors  and  their  associated  instrumenta 
tion.  The  design  of  the  test  hardware  will  now  be  discussed. 

Combustor  Design 

The  basic  combustor  chamber  was  of  6"  diameter;  a  size  suited  to 
the  test  facilities  available.  A  given  test  combustor  was  made  up  of 
three  components,  as  shown  in  Fig  24.  These  components  consisted  of 
a  10"  long  convergent  inlet  section,  the  chamber  proper,  and  a  con¬ 
vergent  exhaust  nozzle.  The  convergent  inlet  section  was  selected 
rather  than  the  usual  inlet  pipe  in  order  to  obtain  a  more  uniform 
inlet  velocity  profile  and  a  reduced  level  of  turbulence.  The  com¬ 
bustor  chamber  consisted  of  one  master  section  to  which  additional 
sections  could  be  added  to  vary  overall  length.  The  master  section 
was  18"  long  and  instrumented  with  15  wall  static  pressure  taps.  The 
chamber  was  machined  internally  to  provide  a  truly  circular  cross- 
section.  For  some  later  tests,  11  external  wall  thermocouples  were 
added  along  the  length  of  this  master  section.  For  the  majority  of 
tests,  a  12"  uninstrumented  chamber  was  added  to  the  master  section, 
yielding  a  30"  overall  length  and  a  chamber  length: diameter  ratio  of  5 
This  value  was  selected  because  of  its  relevance  to  current  ramjet 
engine  designs,  and  also  it  provided  sufficient  length  to  avoid 
interference  with  the  anticipated  recirculation  zone  sizes.  A  prime 
geometric  variable  is  the  step  height.  Thus,  six  inlet  sections  were 
fabricated,  yielding  dump  plane  diameters  of  2.5,  3.0,  3.5,  4.0,  4.5, 
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and  5.0  inches;  correspondingly  the  step  height  could  be  varied  from 
1.75  to  0.5  inches  by  increments  of  0.25  inches. 

In  order  to  vary  the  entry  velocity  over  a  meaningful  range,  three 
convergent  exit  nozzles  of  area  ratios  A*/A  equal  to  0.40,  0.50,  and 
0.60  were  utilized;  thus  together  with  straight-pipe  operation 
(A*/A  =  1.0)  four  nozzle  sizes  were  available.  Subsequently,  opera¬ 
tion  of  the  combustor  in  a  straight  pipe  mode  yielded  very  scattered 
data,  and  two  more  nozzles  (A*/ A  =  0.20  and  0.80)  were  fabricated 
to  provide  a  wider  range  of  velocity  values. 

All  components  were  joined  together  by  flanges.  Considerable 
care  was  taken  during  assembly  to  ensure  the  concentricity  of  the 
inlet  and  combustor  at  the  dump  plane.  This  was  to  ensure  that  a 
uniform  step  height  was  maintained.  Previous  experience  has  shown 
that  lack  of  concentricity  could  result  in  partial  flame  extinction 
instead  of  a  precise,  complete  blow-off. 

The  inlet- to-combustor  joint  incorporated  a  specially  designed 
plate  to  permit  fuel  injection  directly  into  the  annular  RZ,  A  number 
of  tests  were  run  with  secondary  fuel  injection  in  order  to  assess  the 
air  flow  entrained  into  the  recirculation  zone.  Typical  injection 
plates  are  shown  in  Fig  25;  six  plates  were  manufactured,  one  for  each 
dump  size. 

The  available  test  hardware  configurations  are  summarized  in 
Table  I.  All  components  were  fabricated  from  steel. 

Test  Facility 

This  investigation  was  part  of  a  continuing  series  of  combustor 


studies  carried  out  in-house  by  the  Ramjet  Engine  Division  of  the  USAF 
Aero  Propulsion  Laboratory  (AFAPL).  Consequently,  the  AFAPL  test 


Figure  25.  Secondary  Injection  Plates 


facilities  were  utilized  for  this  investigation.  A  brief  discussion 
of  relevant  parts  of  this  facility  follows. 

The  combustor  test  rig  is  mounted  on  a  floating  thrust  stand 
which  is  suspended  from  four  flexures;  a  load  cell  measures  the  thrust 
forces  associated  with  the  combustor.  The  stand  is  illustrated  in 
Fig  26.  The  air  supply  to  the  test  rig  is  supplied  through  12  two-inch 
diameter  flexible  hoses.  The  thrust  stand  is  statically  calibrated 
periodically,  using  a  reference  load  cell.  Additionally,  prior  to 
actual  testing,  the  accuracy  of  the  thrust  reading  is  checked  under 
no- flow  and  flowing  conditions. 

The  air  supply  and  exhaust  systems  are  shown  schematically  in 
Fig  27.  The  air  is  obtained  from  a  300  psi  laboratory  supply  and  may 
be  heated  by  an  indirectly  fired  gas  furnace  or  by  an  in-line  combus¬ 
tion  heater.  The  latter  heater  is  a  modified  J-85  combustion  chamber 
and  is  operated  on  ethylene  with  an  oxygen  make-up  flow  added  upstream 
of  the  heater.  However,  to  avoid  the  effects  of  vitiation  by  combus¬ 
tion  products,  only  the  clean  air  heater  was  used  in  this  test  series. 
The  air  mass  flow  rate  to  the  combustor  is  measured,  using  flange-tap 
square-edge  orifice  plates,  and  two  orifice  installations  are  available, 
depending  on  the  mass  flow  level. 

After  leaving  the  combustor,  the  exhaust  gases  are  cooled  by  a 
water-spray  and  sucked  out  through  a  high  capacity  exhauster  system. 

The  exhaust  pressure  may  be  reduced  down  to  pressures  of  3  psia.  The 
only  hydrocarbon  fuel  readily  available  in  large  quantities  was 
commercial  grade  ethylene,  and  the  bulk  of  the  experimental  work  was 
performed  using  this  gas.  Subsequently,  supplies  of  this  gas  were 
interrupted,  and  several  tests  were  later  run  using  JP-4. 
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Figure  26.  Ramjet  Test  Stand 


For  some  tests  with  ethylene,  both  mainstream  injection  and 
secondary  injection  into  the  RZ  were  required.  Consequently,  the 
bottled  gas  supply  was  split  into  two  feed  lines  and  the  fuel  flow 
measured  by  Cox  turbine- type  flowmeters.  Once  again  calibration  of 
these  flowmeters  was  periodically  accomplished,  using  a  reference 
flowmeter.  When  JP-4  was  used,  only  one  feed  line  was  required;  JP-4 
was  pumped  from  a  main  storage  tank  and  the  flow  rate  also  measured 
by  a  Cox  flowmeter. 

The  data  acquisition  system  is  shown  schematically  in  Fig  28. 

Thus,  incoming  sensor  information  can  be  fed  from  a  patch  panel  to 
either  a  Hewlett-Packard  digital  data  system  or  a  Mod  Comp  II  computer 
controlled  data  system.  This  latter  system  samples  at  a  rate  of  about 
5000  channels  per  sec.  Each  data  point  is  the  average  of  about  45 
separate  scans  of  the  data  system  in  order  to  average  out  electrical 
noise. 

In  addition  to  acquiring,  storing,  and  processing  all  raw  data, 
the  Mod  Comp  II  system  displays  on  a  video  screen,  mounted  on  the  test 
console,  all  parameters  needed  to  run  the  combustor  at  the  required 
test  conditions.  Furthermore,  at  the  end  of  any  given  test  run, 
selected  data  is  displayed  on  a  Tektronix  graphics  display  terminal; 
hard  copies  of  this  data  are  readily  available,  and  selected  quantities 
may  be  cross-plotted  to  further  aid  analysis  of  the  data. 

The  gas  flow  path  through  the  test  rig  and  some  associated 
instrumentation  are  shown  in  Fig  29.  Initially,  a  flow  straightening 
honeycomb  and  several  flow  smoothing  grids  were  installed;  however, 
due  to  occasional  occurrences  of  "flashback,"  these  components 
received  substantial  damage  and  were  subsequently  removed  from  the 
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Figure  29.  Gas  Flow  Paths  Through  Test  Rig 


rig.  Fuel  was  injected  into  the  stream  some  seven  feet  upstream  of 
the  dump  station  through  a  32-point  injection  system  to  ensure  that  a 
uniform  fuel -air  mixture  was  obtained.  The  injector  consisted  of 
eight  radial  spokes,  each  carrying  four  holes,  with  injection  perpen¬ 
dicular  to  the  flow  direction. 

Experimental  Investigations 

Three  distinct  test  series  were  performed  as  follows: 

a.  Flame  stability  of  plain  dump  combustors. 

b.  Blow-out  limits  with  additional  fuel  injection  into  RZ. 

c.  Flame  visualization  studies  with  a  quartz  combustor. 

The  first  and  major  investigation  was  the  effect  of  geometric  variables 
on  flame  stability.  The  configurations  tested  are  shown  in  Table  I. 

The  baseline  inlet  configuration  was  the  3.5"  diameter  dump  size,  and 
the  effect  of  varying  this  inlet  size  from  2.5"  to  4.5"  was  of  prime 
interest,  A  standard  procedure  that  was  executed  prior  to  any  com¬ 
bustion  tests  was  to  check  the  functioning  of  the  test  rig  over  a  wide 
range  of  cold  mass  flows.  At  each  mass  flow,  a  read-out  of  key  data 
was  obtained,  thus  providing  a  final  check  on  the  integrity  of  line 
connections  and  data  acquisition  system  performance.  An  incidental 
benefit  was  the  accumulation  of  a  broad  data  base  on  the  fluid  mech¬ 
anics  of  dump  combustors. 

As  noted  earlier,  the  majority  of  tests  was  performed  using 
ethylene  (95%  pure),  supplies  of  which  were  readily  available,  Two  main 
problems  were  encountered  with  this  gas.  These  were:  the  violent  and 
explosive  combustion  conditions  encountered  in  the  determination  of  rich 
blow-out  limits  and  the  occurrence  of  "flashback"  phenomenon.  The 
determination  of  rich  blow-out  limits  was  difficult  because  of  the 
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increasingly  unsteady  nature  of  the  combustion  process  as  the  rich 
blow-off  limit  was  approached;  this  led  to  poor  repeatability  of  data. 
Furthermore,  the  accumulation  of  fuel-rich  gases  in  the  exhaust  ducting 
led  to  several  explosions  from  which  minor  mechanical  damage  resulted. 
Consequently,  rich-limit  testing  was  soon  eliminated.  The  second 
phenomenon  experienced  with  ethylene  was  combustion  upstream  of  the 
main  burner.  This  phenomenon  occurred  occasionally  during  ignition 
transients  or  when  adjusting  fuel  and  air  mass  flows  at  low  air  flow 
rates,  and  also  with  small  area  exit-nozzles.  The  mechanism  of  this 
phenomenon  was  not  explored  nor  its  brief  occurrence  documented;  when 
it  occurred  its  effects  were  very  damaging,  leading  to  destruction  of 
turbulence-reducing  screens,  pitot  probes,  and  even  part  of  the  5"  inlet 
converging  section.  Plee  and  Mellor  (Ref  95)  have  recently  reviewed 
this  phenomenon. 

Opportunity  was  taken  to  obtain  flame  stability  limits  with 
pre-mixed  JP-4  and  air.  In  this  case  the  flow  at  the  rich  blow-off 
limit  was  much  smoother,  and  repeatable  data  were  obtained.  Also,  the 
occurrence  of  flashback  was  much  reduced. 

For  either  fuel,  the  experimental  procedure  to  determine  blow-off 
limits  was  the  same.  At  a  given  inlet  temperature,  the  combustor  was 
tested  over  a  range  of  mass  flows.  For  a  given  mass  flow,  the  combustor 
was  ignited  at  a  high  fuel-air  ratio  (typically  0.05),  and  the  mass  flow 
readjusted  to  the  desired  condition.  The  lean  blow-off  limit  was 
determined  by  incrementally  reducing  the  fuel  flow  while  maintaining 
the  desired  air  mass  flow  level  until  extinction  occurred.  This  pro¬ 
cess  had  to  be  accomplished  expeditiously  to  avoid  overheating  of  the 
combustor,  but  slowly  enough  to  ensure  that  steady  flow  prevailed  after 
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each  incremental  change  of  fuel  flow.  Several  tests  were  performed 
at  the  came  flow  conditions  to  ensure  that  repeatable  data  were 
obtained.  A  Kistler  pressure  transducer  located  at  the  dump  plane  was 
used  to  give  a  precise  indication  of  blow-off.  The  visual  display  of 
the  transducer  output  clearly  indicated  the  cessation  of  combustion  at 
blow-off.  A  thermocouple  was  also  located  in  the  RZ,  but  the  Kistler 
display  was  a  much  superior  indicator  of  blow-off.  Rich  blow-off 
limits  were  determined  in  a  similar  fashion.  During  the  test  phase, 
two  experimental  constraints  were  generally  encountered.  Considerable 
care  was  required  both  to  maintain  a  constant  inlet  temperature  as  the 
air  mass  flow  varied  and  to  avoid  combustion  instability  effects.  The 
temperature  problem  was  due  to  the  thermal  inertia  associated  with  the 
extreme  length  of  the  air  supply  pipe  from  the  air  heater  to  the  rig 
(several  hundred  feet).  This  problem  was  overcome  by  running  test 
points  alternately  at  high  and  low  air  mass  flows.  Combustor  "screech" 
was  encountered  at  high  pressure,  high  mass  flow  conditions,  and  was 
generally  avoided  by  testing  at  reduced  pressures.  It  was  also  noted 
that  the  "screech"  generally  disappeared  as  the  lean  blow-off  limit 
was  approached.  However,  this  phenomenon  was  of  sufficient  interest 
that  a  separate  investigation  using  the  same  hardware  has  subsequently 
been  initiated. 

As  part  of  a  general  upgrading  of  facility  equipment,  the  clean 
air  heater  initially  used  in  this  investigation  was  subsequently 
replaced  by  a  more  energy-efficient  unit.  During  the  replacement 
phase,  some  of  the  flame  stability  tests  were  repeated,  using  a 
vitiated  heater  with  oxygen  make-up.  A  considerable  narrowing  of  the 
stability  range  was  immediately  noted,  and  no  further  testing  was 
accomplished  with  this  heater. 
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The  second  type  of  investigation  performed  was  to  study  the 
effect  of  additional  fuel  injection  into  the  RZ.  As  noted  previously, 
the  motivation  for  this  investigation  was  that  the  determination  of  the 
lean  blow-off  limits  of  a  split  injection  system  enables  an  estimate  to 
be  made  of  the  air  entrainment  rate  into  the  recirculation  zone.  The 
procedure  for  determining  the  lean  blow-off  limits  in  this  case  was 
similar  to  that  outlined  previously  except  that  the  fuel  injection  rate 
into  the  RZ  was  maintained  at  a  constant  value  while  the  main  fuel  flow 
was  varied  to  obtain  the  lean  blow-off  limit.  The  fuel  injection  plates 
utilized  for  secondary  fuel  injection  were  of  the  type  shown  in  Fig  25. 
Secondary  fuel  is  injected  axially  through  concentric  rings  of  holes. 

A  large  number  of  holes  was  utilized  to  keep  the  fuel  injection  velocity 
low  so  as  to  minimize  interference  with  the  main  recirculation  zone. 
Furthermore,  for  each  test  geometry  the  effect  of  secondary  fuel  injec¬ 
tion  on  the  recirculation  zone  was  studied  by  observing  the  effect  of 
fuel  flow  on  the  wall  static  pressure  distribution.  Gaseous  nitrogen 
was  substituted  for  fuel  in  these  tests  to  conserve  fuel.  In  general, 
no  significant  effect  of  gas  flow  on  the  structure  of  the  recirculation 
zone  was  observed  and,  indeed,  no  effect  was  anticipated  because  such 
fuel  flow  rates  were  very  small.  Two  interesting  effects  were  observed 
with  secondary  fuel  injection.  Firstly,  two  types  of  blow-off  were 
observed.  For  relatively  low  secondary  fuel  flows  a  normal  lean  blow- 
off  was  obtained  as  the  main  fuel  flow  was  reduced.  However,  at  higher 
secondary  fuel  flows,  as  the  lean  blow-off  was  approached,  a  change  in 
combustion  mode  was  observed,  and  the  combustor  would  not  extinguish 
until  the  secondary  flow  was  reduced.  This  latter  mode  had  to  be 
avoided  in  order  to  prevent  ambiguity  in  the  desired  test  data. 


Another  significant  observation  was  that  under  certain  conditions  the 
injection  of  small  quantities  of  secondary  fuel  led  to  a  much  quieter- 
running  combustor;  thus,  secondary  injection  may  suppress  some  insta¬ 
bility  mechanisms. 

A  final  test  series  was  performed  to  observe  the  blow-off 
phenomenon.  For  this  series,  an  annealed  quartz  combustor  was  utilized, 
and  the  installation  is  shown  in  Fig  30.  It  will  be  seen  that  the  com¬ 
bustor  is  secured  between  the  upstream  and  downstream  flanges  by  four 
spring-loaded  steel  rods.  Each  end  of  the  quartz  combustor  was  sealed 
with  a  ring  gasket  fabricated  from  "carborundum  fiber  frax"  material  to 
which  was  applied  a  silicone  rubber  coating.  After  curing  of  the  coat¬ 
ing,  an  effective  high  temperature  gasket  was  obtained.  This  method  of 
installing  a  quartz  chamber  functioned  quite  satisfactorily  under  com¬ 
bustion  conditions  provided  that  "screech"  was  avoided.  Flame  blow-off 
was  observed  with  three  inlet  sizes:  2.5",  3.5",  and  4.5".  Additionally, 
some  studies  were  made  of  flame  blow-off  phenomena  when  secondary  fuel 
injection  was  employed.  All  flame  visualization  studies  were  photo¬ 
graphically  recorded,  usually  at  a  frame  speed  of  500  frames  per  second. 

The  results  obtained  from  the  various  test  series  outlined  above 


will  now  be  described. 


V  RESULTS  AND  DISCUSSION 


It  will  be  recalled  from  the  preceeding  Chapter  IV  that  three 
distinct  test  series  were  performed  as  follows: 

a.  Determination  of  the  flame  stability  limits  of  plain  dump 
combustors,  using  both  ethylene  and  JP-4  as  fuels. 

b.  Determination  of  lean  extinction  limits  with  additional 
injection  into  the  recirculation  zone. 

c.  Flame  visualization  studies  using  a  transparent  quartz 

combustor. 

In  addition  to  the  three  basic  test  series,  a  cold-flow  test  was 
run  to  determine  the  velocity  profiles  associated  with  each  combustor 
inlet.  These  tests  were  run  over  a  range  of  inlet  Reynolds  and  Mach 
numbers.  The  results  are  shown  in  Appendix  D,  and  it  will  be  noted 
that  substantially  constant  velocity  profiles  were  obtained. 

The  results  of  each  basic  test  series  will  now  be  presented  and 
discussed. 


Flame  Stability  Tests 

In  order  to  investigate  the  effects  of  the  primary  variables  of 
temperature,  pressure,  step  height,  and  velocity  on  the  stable  operating 
range,  it  is  obviously  desirable  to  vary  these  quantities  independently 
of  one  another.  The  first  three  variables  are  easily  controlled;  how¬ 
ever,  for  a  given  combustor  configuration  operating  at  a  given  inlet 
temperature  and  with  a  choked  exit  nozzle,  the  Mach  number  at  the  dump 
plane  is  dependent  on  the  effective  heat  release  in  the  combustor. 

At  the  blow-out  condition,  the  dump  plane  velocity  will  typically 

U  “ 

satisfy  Eq(49),  i.e.,  ♦LB0  «  — TTTT  ’  and  will  simultaneously 

P_  n  T„ 
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correspond  to  a  value  determined  by  the  effective  heat  release  in  the 
combustor.  It  will  be  recalled  from  Chapter  IV  that  for  a  given  test, 
the  lean  blow-out  condition  was  determined  by  progressively  reducing 
the  fuel  flow  while  the  inlet  temperature  and  mass  flow  were  maintained 
constant.  For  successive  tests  at  a  given  inlet  temperature  but  with 
various  mass  flow  levels,  it  was  found  that  the  dump  plane  Mach  number 
at  blow-out  remained  practically  constant  over  the  relatively  limited 
range  of  variation  of  fuel-air  ratios  which  occurred  in  determining  the 
lean  extinction  points.  Thus,  for  a  given  inlet  temperature  the  dump 
plane  velocity  at  lean  extinction  is  approximately  constant  even  though 
the  mass  flow  (and  therefore  the  pressure)  is  varied.  Obviously,  for  a 
given  nozzle  throat  size  the  dump  plane  velocity  will  decrease  with 
increasing  dump  diameter,  and  conversely  for  a  given  dump  diameter  this 
velocity  will  rise  if  the  throat  area  is  increased.  With  this  introduc¬ 
tory  note  on  flow  conditions  concluded,  one  may  now  proceed  to  consider 
the  data.  The  results  obtained  with  ethylene  will  be  considered  first. 

Tests  with  Ethylene.  The  test  matrix  performed  with  ethylene  was 
presented  in  Table  I.  It  will  be  noted  that  all  dump  sizes  were 
successfully  run,  the  dump  plane  diameters  varying  from  2.5"  to  5.0". 

In  the  latter  case,  it  was  interesting  to  observe  that  successful 
flame  stabilization  and  combustor  operation  was  achieved  with  only  a 
1/2"  step.  Unfortunately,  only  limited  testing  was  possible  as  the 
5.0"  convergent  section  received  severe  damage  due  to  "flashback"  of 
the  flame  upstream  of  the  combustor. 

The  experimental  data  acquired  for  a  typical  dump  plane 
diameter  (3.5")  are  shown  in  Figs  31,  32,  and  33;  each  figure  corre¬ 
sponds  to  a  given  nozzle  area  ratio.  Only  lean  extinction  data  are 


TABLE  I  COMBUSTOR  GEOMETRIES  TESTED 


All  of  the  above  testing  was  with  clean  air. 
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31.  Lean  Extinction  Limits:  3.5"  Dump  Dia.,  0.40  Nozzle 
Area  Ratio;  Ethylene  Fuel 
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Figure  32,  Lean  Extinction  Limits:  3.5"  Dump  Dia.,  0.50  Nozzle 
Area  Ratio;  Ethylene  Fuel 


105 


shown  because  of  the  hazardous  nature  of  the  rich  limit  testing. 
Considering  Fig  32,  it  will  be  recalled  that  a  line  of  constant  inlet 
temperature  operation  also  corresponds  effectively  to  a  line  of  constant 
velocity.  It  is  immediately  apparent  that,  at  a  given  temperature  and 
velocity,  the  effect  of  increasing  pressure  is  to  reduce  the  lean 
extinction  fuel-air  ratio.  This  variation  is  also  in  accord  with  that 
anticipated  from  Eq(71). 

The  effect  of  velocity  on  the  lean  extinction  limit  was 
evaluated  by  two  methods.  The  first  was  to  compare  the  data  from  com¬ 
bustors  of  the  same  step  height,  operated  at  identical  inlet  pressure 
and  temperature  conditions,  but  with  different  exit  nozzle  throat  sizes. 
Such  a  comparison  is  shown  in  Fig  34,  and  it  is  apparent  that  inlet 
velocity  does  not  exert  a  strong  effect  on  lean  extinction  limits. 
Unfortunately,  only  one  such  comparison  could  be  made  because  of  the 
difficulty  of  expeditiously  setting-up  identical  temperature  conditions 
at  the  entrance  of  each  of  the  desired  combustor  configurations.  This 
was  due  to  the  thermal  inertia  of  the  heater-inlet  pipe  system,  A 
second  approach  was  therefore  adopted,  and  this  consisted  of  operating 
a  given  combustor  in  a  straight-pipe  (i.e.,  no  exit  nozzle)  unchoked 
mode.  Under  these  conditions,  the  inlet  pressure  and  temperature  were 
maintained  approximately  constant,  and  the  velocity  varied  by  changing 
the  mass  flow.  Typical  data  are  shown  in  Fig  35,  using  Mach  number  for 
convenience;  thus  the  lean  limit  is  relatively  insensitive  to  velocity, 

A  regression  analysis  of  the  data  showed  that  a  slight  increase  of  lean 
limit  fuel-air  ratio  with  velocity  is  present. 

Returning  to  Fig  31,  it  is  quite  clear  that  temperature 
exerts  a  major  effect  on  the  lean  limit  fuel -air  ratio.  Thus  with 
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of  Inlet  Mach  Number  on  Lean  Extinction  Limit;  Ethylene  Fuel 


increasing  temperature,  the  lean  extinction  fuel-air  ratio  decreases 
markedly.  This  observation  is  also  in  accord  with  Eq(71).  The  effect 
of  step  height  variation  on  the  lean  limit  is  illustrated  5n  Fig  36: 
in  this  case  the  lean  extinction  limits  of  two  combustors,  operated 
under  identical  inlet  conditions,  are  compared.  The  expected  beneficial 
effect  of  increased  step  height  is  clearly  seen,  and  this  observation 
also  accords  with  Eq( 71 ) . 

Following  the  conclusion  of  the  tests  with  ethylene,  two  sets  of 
tests  were  run  using  JP-4  as  fuel.  These  results  will  now  be  described. 

Tests  With  JP-4.  The  first  test  series  performed  with  JP-4  is 
presented  in  Table  I.  It  will  be  noted  that  once  again  three  nozzle 
area  ratios  were  available;  namely,  0.40,  0.50,  and  0.60.  A  second  set 
of  tests  was  performed  with  two  additional  nozzle  area  ratios  equal  to 
0.20  and  0.80,  respectively.  Some  difficulty  was  experienced  with  the 
smaller  nozzle  because  of  the  occurrence  of  "flashback"  at  the  low  dump 
plane  velocities  corresponding  to  this  nozzle  area  ratio,  and  only 
limited  data  were  obtained.  An  additional  feature  of  the  tests  with 
JP-4  was  that  rich  blow-out  limits  were  also  obtained.  In  the  case  of 
JP-4,  the  rich  limit  extinction  was  relatively  quiet  and  less  hazardous 
than  with  ethylene.  Also,  the  availability  of  a  larger  exhaust  nozzle 
permitted  full  stability  loops  to  be  demonstrated  for  the  first  time 
for  dump  configurations, 

A  typical  set  of  data  acquired  for  a  given  configuration  is 
shown  in  Fig  37,  and  both  lean  and  rich  extinction  limits  are  shown. 

The  lean  limit  variation  is  similar  to  that  obtained  previously  with 
ethylene.  It  will  be  noted  that  increasing  pressure  improves  the 
stable  operating  range.  The  effect  of  velocity  on  the  lean  limit 
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Figure  36.  Effect  of  Step  Height  on  Lean  Extinction  Limit; 
Ethylene  Fuel 
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fuel-air  ratio  was  only  studied  in  this  case  by  operating  combustors  of 
the  same  diameter  ratio  (d/D)  at  identical  inlet  pressure  and  temperature 
conditions  but  with  different  nozzle  throat  sizes.  A  comparison  is  shown 
in  Fig  38,  and  once  again  it  is  apparent  that  velocity  has  only  a  small 
effect  on  the  lean  extinction  limit.  Returning  to  Fig  37,  it  is  clear 
that  increasing  temperature  exerts  a  major  beneficial  effect  on  the  lean 
extinction  limit.  The  effect  of  step  height  variation  is  shown  in  Fig  39, 
and  once  again  increasing  step  height  reduces  the  lean  extinction  fuel- 
air  ratio.  Thus  the  effects  of  pressure,  temperature,  velocity,  and  step 
height  are  very  similar  to  those  observed  with  ethylene. 

One  noteworthy  effect  encountered  in  the  JP-4  test  series  was 
the  occurrence  of  combustor  "screech"  at  high  inlet  temperatures.  In 
one  instance  this  produced  considerable  scatter  in  the  lean-extinction 
data  as  shown  in  Fig  40.  The  "screech"  condition  was  eliminated  by 
moving  the  fuel  injectors  upstream  into  the  highly  turbulent  area  where 
air  enters  the  test  rig.  The  multi-hole  finger  injectors  were  also 
replaced  with  fan-spray  injector  elements.  The  improvement  of  the  data 
that  occurred  with  "screech"  elimination  is  also  shown  in  Fig  40. 

Another  test  excursion  was  to  investigate  the  effect  of 
combustor  length  to  diameter  ratio  0-c/D)  on  the  lean  stability  limit. 

It  is  well  known  that  for  conventional  stabilizers,  an  increase  in 
combustor  length  results  in  a  decreased  stable  operating  range.  The 
reason  for  this  decrease  is  usually  attributed  to  increased  flow 
perturbations  arising  from  acoustic  effects.  The  significance  of  this 
variable  was  assessed  by  operating  the  2.5"  dump  combustor  with 
combustor  lengths  of  18",  30",  and  36".  The  corresponding  results 
are  shown  in  Fig  41;  it  will  be  noted  that  only  a  minor  reduction  in 

stable  operating  range  occurred  in  this  case. 
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Figure  40.  Effect  of  Screech  on  Quality  of  Data; 
JP-4  Fuel 
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Another  brief  test  excursion  was  to  investigate  the  effect  of 
using  vitiated  air  on  the  stable  operating  range.  In  this  case  an 
ethylene- fueled  combustion  heater  was  used  with  oxygen  make-up. 

Despite  the  replenishment  of  the  oxygen,  a  distinct  reduction  in  stable 
operating  range  was  experienced  as  shown  in  Fig  42.  This  result  was 
rather  surprising  considering  that  the  incoming  air  was  only  heated 
from  600  to  1000°R. 

As  a  final  test  goal,  it  was  decided  to  determine  the 
complete  stability  loop  for  given  combustor  configurations.  This  was 
successfully  accomplished  by  operating  these  combustors  at  very  low 
pressures,  typically  at  less  than  10  psia.  The  results  obtained  are 
shown  in  Figs  43  and  44,  and  they  demonstrate  conclusively  that  sta¬ 
bility  loops  do  exist  for  the  sudden  expansion  configuration.  This 
demonstration  concluded  the  flame  stability  tests  of  the  plain  dump 
combustor.  No  general  discussion  has  been  given  of  data  obtained  with 
the  3",  4",  and  5"  diameter  inlets.  This  was  due  firstly  to  the  fact 
the  5"  inlet  was  destroyed  due  to  "flashback"  early  in  the  test  program. 
Secondly,  the  3"  and  4"  inlets  were  tested  only  with  ethylene,  and  the 
limited  data  obtained  simply  confirmed  the  results  obtained  with  the 
2.5",  3.5",  and  4.5"  combustors.  The  second  test  series  which  utilized 
direct  injection  into  the  recirculation  zone  will  now  be  described. 

Lean  Extinction  Tests  With  Split  Flow  Injection 

The  purpose  of  these  tests  was  to  determine  the  amount  of  air 
entrained  and  consumed  in  the  recirculation  zone.  It  will  be  recalled 
that  in  order  to  estimate  the  fraction  of  air  entrained,  it  is  nec¬ 
essary  to  use  a  split  fuel  injection  system.  Thus  fuel  is  injected 
both  into  the  main  air  flow  and  directly  into  the  recirculation  zone. 
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Figure  42.  Effect  of  Vitiation  on  Stable  Range;  JP-4  Fuel 


Fuel /Air  Ratio 


Figure  43.  Complete  Stability  Loop;  JP-4  Fuel 


Combustor 

Pressure 


The  fuel  injection  plates  used  for  secondary  fuel  admission  are 
illustrated  in  Fig  25.  Tests  were  performed  with  these  inlet  sizes: 


2.5",  3.5",  and  4.5".  The  data  for  the  2.5"  inlet  is  shown  in  Fig  45. 

It  will  be  observed  that  the  tests  were  performed  for  various  air  mass 
flow  and  temperature  levels.  The  data  scatter  increases  slightly  with 
the  more  turbulent  conditions  experienced  at  higher  mass  flows.  Atten¬ 
tion  is  also  drawn  to  the  very  small  fuel  flow  injected  into  the  recir¬ 
culation  zone.  This  required  precise  fuel  flow  metering  for  accurate 
results.  It  will  also  be  seen  that  all  test  conditions  yielded  parallel 
straight  line  relationships  between  the  primary  and  RZ  fuel  flows  at 
lean  extinction,  which  implies  that  the  recirculated  mass  flow  fraction 
does  not  vary  with  the  mass  flow  rate. 

Testing  at  the  smaller  step  sizes  corresponding  to  the  3.5"  and 
4.5"  inlet  diameter  had  to  be  accomplished  deftly  and  expeditiously  to 
prevent  the  combustor  wall  from  overheating.  The  data  for  the  4.5" 
diameter  inlets  are  shown  in  Fig  46.  The  data  for  the  3.5"  inlet  had 
to  be  discarded  because  of  erroneous  secondary-fuel  flowmeter  readings. 

The  test  work  was  concluded  by  a  brief  visual  study  of  the  lean 
extinction  process  utilizing  a  quartz  combustor. 

Quartz  Combustor  Tests 

The  installation  of  the  quartz  combustor  is  shown  in  Fig  30.  High 
speed  movies  (500  frames  per  second)  were  taken  of  the  combustion  process 
with  the  aim  of  determining  the  mechanism  of  lean  extinction.  In  gen¬ 
eral  ,  the  precise  manner  of  blow-out  could  not  be  discerned,  primarily 
because  of  the  rapid  occurrence  of  extinction  which  could  not  be  cap¬ 
tured  at  the  frame  speed  used.  Another  difficulty  was  the  oscillatory 
nature  of  the  combustion  process.  Under  conditions  which  were 


o 

o 


00 

o 

o 


o 

o 


o 

o 


C\J 

o 

o 


£  I 

o  u 

0J 

Li-  </) 

OC.  r—  {/) 
<D  JO 
3  »— 


Finure  45.  Primary  and  RZ  Fuel  Flows  at  Lean  Exti 
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Figure  46.  Primary  and  RZ  Fuel  Flows  at  Lean  Extinction: 
4.5"  Dump 


apparently  quite  steady  to  the  naked  eye,  a  frame-by-frame  examination 
of  the  film  showed  that  a  bright  reaction  zone  oscillated  to  and  fro 
along  the  combustor.  This  phenomenon  may  have  been  related  to  the  fluc¬ 
tuating  nature  of  the  reattachment  process.  In  any  event,  the  sequence 
of  the  extinction  process  was  not  clearly  discernible,  and  extinction 
occurred  quite  abruptly. 

Discussion  of  Flame  Stability  Data 

It  will  be  recalled  from  the  discussion  of  dump  combustor  stability 
given  in  Chapter  II  that  a  commonly  accepted  parameter  for  correlating  sta¬ 
bility  data  is  the  X  parameter  given  in  Eq(46).  Attempts  to  correlate 
both  the  ethylene  and  JP-4  data  obtained  during  this  investigation  are 
shown  in  Figs  47  to  50.  It  is  apparent  that  neither  the  effects  of 
dump  height  nor  of  inlet  temperature  were  successfully  correlated. 

Attention  was  therefore  focused  on  the  expression  given  in  Eq(71)  for 
the  equivalence  ratio  at  lean  blow-out.  For  the  assumptions  made, 
namely  that  both  the  ethylene  and  JP-4  combustion  processes  were  des¬ 
cribed  by  a  single-step  reaction  with  an  apparent  activation  energy, 

E  =  42  Kcal/mole,  and  reaction  order  n  =  1.75,  the  expression  given 
in  Eq ( 71 )  is,  with  T  put  equal  to  Tt  , 

0  n 


Turning  first  to  the  ethylene  data,  it  was  found  that  for  the  range  of 
inlet  conditions  tested,  the  data  could  be  correlated  by  an  expression 
of  the  form  similar  to  that  of  Eq ( 71 )  but,  as  previously  anticipated, 
with  modified  exponents.  Thus,  for  ethylene,  one  obtains 
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Figure  47.  Trial  Correlation  by  X  Parameter:  Effect  of  Step  Size 
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Figure  48.  Trial  Correlation  by  X  Parameter: 
Effect  of  Inlet  Temperature 


Figure  49.  Trial  Correlation  by  X  Parameter: 
Effect  of  Inlet  Temperature 
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Lean  stability  limits  have  been  expressed  in  the  "reduced"  form  of  these 
variables,  namely 
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where  =  0.16.  The  ethylene  data  from  the  2.5",  3.5",  and  4.5"  dump 
diameter  combustors  are  shown  in  Figs  51,  52,  and  53,  respectively.  The 
effect  of  step  size  h  is  included  in  the  overall  data  correlation  shown 
in  Fig  54:  this  correlation  is  valid  over  the  range  of  inlet  conditions 
tested;  namely, 

10  <  PQ  <  40  psia 

1000  <  Tt  <  1250°R 

0 

0.2  <  Mq  <  0.9 

A  similar  approach  was  taken  with  the  JP-4  test  data,  and  successful 
correlations  of  the  form  given  in  Eq(71)  were  developed.  In  this  case, 
the  appropriate  exponents  were 
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and  in  terms  of  reduced  variables,  this  becomes 
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The  overall  correlation  using  "reduced  variables"  is  shown  in  Fig  55,  and 
this  is  valid  for  the  same  range  of  conditions  previously  noted  for  the 
ethylene  test  series.  Once  again,  a  very  satisfactory  data  correlation 
is  obtained. 
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Figure  52.  Lean  Stability  Limit  Correlation:  3.5"  Dump  Diameter; 
Ethylene  Fuel 


132 


izes;  OP-4  Fue 


Thus  the  data  correlations  obtained  are  in  general  agreement  with  the 
form  anticipated  by  stirred  reactor  theory.  Indeed  the  exponents  for 
h,  PQ,  and  Tt  are  quite  close  to  the  anticipated  values.  However, 
an  interesting  feature  of  the  correlations  given  in  Eq(73)  and  Eq ( 75 ) 
is  the  relative  insensitivity  of  the  weak  limit  equivalence  ratio  to 
velocity  variation.  Although  this  dependence  is  normally  low  for  bluff 
bodies  (Ref  31),  it  appears  to  be  even  less  in  the  case  of  coaxial  con¬ 
figurations.  Moreover,  inspection  of  Figs  17  and  18  confirms  this  weak 
dependence  for  both  two-dimensional  steps  and  for  coaxial  geometries. 
This  phenomenon  muse  therefore  be  associated  with  the  flow  features 
which  are  unique  to  step  geometries.  The  features  which  are  specific  to 
the  step  compared  to  a  bluff  body  stabilizer  include:  the  presence  of  a 
wall  bounding  the  recirculation  zone,  the  existence  of  a  well -developed 
boundary  layer,  and  a  shear  layer  which  reattaches  to  the  wall.  From 
Fig  17,  it  is  observed  that  the  use  of  different  wall  materials  which 
affect  the  heat  loss  to  the  wall  did  not  appear  to  change  the  velocity 
dependence.  It  is  also  noted  that  the  velocity  enters  the  analysis  by 
way  of  the  mr  term  as  shown  in  Eq(47);  it  is  therefore  presumed  that 
the  specific  aerodynamic  features  of  the  recirculation  zone  flow  which 
affect  entrainment  are  responsible  for  the  low  velocity  dependence. 
Future  investigations  are  needed  to  clarify  this  phenomenon.  For  a 
given  inlet  stagnation  temperature,  an  upper  limit  is  set  on  the  flow 
velocity  at  the  dump  plane  by  the  occurrence  of  choking  at  that  plane. 
Thus  the  term  U  ^  is  bounded.  In  many  cases  it  follows  that 
extinction  will  occur  due  primarily  to  the  existence  of  low  temperature 

(Tt  )  or  low  pressure  (P.)  conditions  rather  than  due  to  blow-off  due 
^o  o 

to  excessive  velocity. 


to 


136 


In  view  of  the  satisfactory  agreement  of  the  data  with  the 
prediction  of  the  stirred  reactor  model,  and  with  regard  to  the  relative 
insensitivity  of  the  lean  extinction  limit  to  velocity,  no  application 
of  alternate  stabilization  criteria,  i.e.,  time-ratio  models,  is  made. 

It  is  appropriate  to  note,  however,  that  many  flame  stability  data  can 
be  explained  equally  well  by  either  stirred  reactor  or  time-ratio  models. 
This  dual  approach  would  also  be  anticipated  to  apply  to  the  flame  sta¬ 
bility  d  a  associated  with  step  geometries.  In  this  case,  however,  a 
time-ratio  model  would  have  to  explain  the  relative  insensitivity  of  the 
lean  limit  to  inlet  velocity  variation.  This,  in  turn,  would  require  a 
characteristic  length  scale  whose  magnitude  was  approximately  propor¬ 
tional  to  inlet  velocity.  Thus,  further  investigations  of  the  shape 
and  structure  of  the  RZ,  and  the  associated  shear  layer,  under  actual 
operating  conditions,  are  of  considerable  interest.  Furthermore,  it 
could  be  anticipated  that  the  parametric  expression  given  in 

Eq ( 24)  would  not  correlate  the  data  due  to  the  demonstrated  weak  effect 
of  the  velocity  term;  such  proved  to  be  the  case.  Indeed  the  most 
useful  parameter  for  collapsing  both  weak  and  rich  limit  data  was  the 
correlating  group  given  in  Eq(75).  Thus  all  the  data  obtained  from  the 
JP-4  tests  are  displayed  in  terms  of  this  correlating  group  in  Fig  56. 
This  correlation  is  much  more  satisfactory  than  that  shown  in  Fig  16. 
Discussion  of  Split-Flow  Injection  and  Quartz  Combustor  Tests 

Useful  data  were  obtained  for  the  2.5"  and  4.5"  diameter  inlets. 

The  recirculation  fraction  obtained  from  the  2.5"  inlet  is  shown  in 
Fig  57.  Also  plotted  are  the  maximum  recirculation  fractions  predicted 
by  Wingfield's  expression  given  in  Eq(49).  Also,  a  correction  to  allow 
for  the  change  in  density  between  hot  and  cold  conditions  is  displayed 
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in  Fig  57.  This  correction  for  density  ratio  is  given  in  Appendix  A, 
Eq(A-30).  Curves  corresponding  to  density  ratios  of  2  and  3  are  given 
in  Fig  57.  It  will  be  seen  that  the  experimental ly  determined  recircu¬ 
lation  fraction  is  in  good  agreement  with  the  expression  of  Wingfield 
if  a  density  ratio  of  approximately  2.5  is  assumed.  A  cross  check  on 
the  density  ratio  was  made  based  on  the  reading  of  the  thermocouple 
installed  in  the  recirculation  zone.  Although  the  purpose  of  this 
thermocouple  was  to  signal  flame  extinction,  it  does  provide  an  approxi¬ 
mate  indication  of  the  RZ  temperature.  Examination  of  the  temperature 
data  shows  that  the  temperature  ratio  at  the  point  of  extinction  lies 
generally  in  the  range  2. 5-3.0.  This  agreement  between  Wingfield's 
expression  and  the  experimental  data  is  gratifying.  Of  course,  the  2.5" 
inlet  corresponds  to  the  smallest  value  of  diameter  ratio  (d/D)  tested. 
Thus  it  is  expected  to  yield  recirculation  fraction  data  approaching 
that  predicted  by  theories  devised  for  industrial  furnaces  of  low 
diameter  ratio.  One  may  then  conclude  that  all  the  air  entering  into 
the  recirculation  orocess  is  in  fact  taking  part  in  the  combustion  pro¬ 
cess.  However,  it  is  considered  that  more  data  from  similar  low  diameter 
ratio  combustors  is  required  to  confirm  this  interesting  result. 

The  4.5"  diameter  inlet  data  are  also  shown  in  Fig  57.  In  this 
case  the  experimental  point  lies  on  the  cold  flow  line  predicted  by 
Wingfield's  expression.  In  this  case  the  combustor  diameter  ratio 
(0.75)  is  much  too  high  for  the  similarity  theories  discussed  in 
Appendix  A  to  apply.  Nevertheless,  the  data  are  of  the  right  order  of 
magnitude,  and  it  is  interesting  to  speculate  whether  the  agreement  is 
fortuitous,  or  if  indeed  the  presence  of  a  hot  shear  layer  close  to  a 
small  step  (0.75")  invalidates  the  density  correction  applicable  to 
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geometries  of  lower  diameter  ratio.  Once  again  further  data  are 
required  to  supplement  these  exploratory  investigations. 

The  experimentally  determined  values  of  mrc/m0  were  not  employed, 
however,  in  the  evaluation  of  reactor  loading  expression  given  in  Eq(50). 
This  is  due  to  the  fact  that  Eq(50)  includes  the  term  A0/V.  The 
magnitude  of  this  term  is  dependent  on  the  length  ratio  (k)  of  the  hot 
recirculation  zone,  and  the  precise  value  of  k  is  uncertain.  However, 
the  applicability  of  the  split-fuel  flow  technique  to  dump  combustors 
has  been  demonstrated.  Also,  the  recirculated  mass  flow  values  derived 
from  these  tests  will  be  of  distinct  value  to  those  engaged  in  creating 
reactor-network  models  of  the  combustor  process. 

Turning  to  the  tests  of  the  transparent  quartz  combustor,  the  prime 
value  of  these  tests  was  twofold.  Firstly,  they  demonstrated  that  the 
combustion  process  in  the  coaxial  combustor  can  be  observed  under  actual 
operating  conditions.  Further  work  is  definitely  required  with  higher 
frame  speeds  to  study  the  actual  mechanism  of  flame  extinction. 

Secondly,  visual  inspection  of  the  movies  on  a  frame-by-frame  basis 
revealed  a  bright  reaction  zone  which  oscillated  rapidly  up  and  down 
the  chamber.  The  recirculation  region  appeared  as  an  annular,  dark 
zone  with  flame  originating  at  the  circumferential  lip  of  the  step. 

The  actual  flame  surface  appeared  as  a  turbulent  brush  intermittently 
anchored  at  the  circumference  of  the  dump  inlet.  Thus  the  basically 
unsteady  nature  of  combustion  in  a  dump  combustor  was  clearly 
visual ized. 

As  a  result  of  these  various  test  series,  several  conclusions 
may  be  drawn.  These  are  itemized  in  the  next  chapter. 
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VI  CONCLUSIONS 


The  research  described  in  this  document  was  conducted  to  explore 
the  stability  of  flames  located  in  the  sudden  expansion  region  of  a 
coaxial  dump  combustor.  The  associated  goals  of  the  investigation  were 
to  determine  the  flame  stability  characteristics  over  a  range  of 
geometries  and  inlet  conditions;  to  model  the  stabilizing  process  and 
derive  a  correlating  parameter;  and  to  explore  a  split  fuel  injection 
technique  to  determine  the  recirculated  mass  flow  fraction.  It  was 
also  desired  to  review  the  literature  concerning  the  flame  stability 
and  aerodynamics  of  coaxial  dump  combustors. 

An  experimental  investigation  was  therefore  performed  to  determine 
the  flame  stability  associated  with  a  6"  diameter  combustor  with  step 
heights  varying  from  0.5"  to  1.75".  Inlet  conditions  were  varied  over 
the  approximate  ranges:  pressure,  10-40  psia;  inlet  total  temperature, 
1000-1250°R;  and  Mach  number,  0.2-0. 9. 

From  this  study  the  following  conclusions  were  drawn: 

1.  The  mechanism  of  flame  stability  in  a  sudden  expansion 
region  is  similar  to  that  observed  with  bluff  body  stabilizers,  and  the 
presence  of  a  wall  bounding  the  recirculation  zone  does  not  produce  any 
marked  differences  in  overall  behavior.  Thus,  as  in  the  case  of  a 
bluff  body,  the  lean  stability  limit  is  enhanced  by  increased  inlet 
pressure,  inlet  temperature,  and  step  height,  and  degraded  by  increased 
inlet  velocity. 

2.  Over  the  range  of  this  investigation,  the  effect  of 
velocity  on  the  weak  extinction  limit  is  very  small  (^gQ  a  U0^). 
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This  weak  dependence  is  a  conmon  feature  of  two-dimensional  step 
geometries  and  coaxial  dump  systems. 


3.  The  dominant  influence  in  flame  stability  is  the  inlet 
temperature;  the  subsidiary  variables  are  step  height,  inlet  pressure, 
and  velocity. 


4.  The  recirculation  zone  associated  with  a  sudden  expansion 
can  be  successfully  modelled  as  an  adiabatic  stirred  reactor  using  a 
one-step  chemical  reaction.  From  this  model  a  correlating  group  for 
the  lean  stability  limit  may  be  derived  in  the  form 
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Spec  fic  correlating  groups  were  derived  for  both  ethylene  and  JP-4 
fueled  systems  and  are  given  by  the  following  equations: 

Ethylene 


4>LBo  0-16 


1 


U 


0.10 


1 


(h/12)0*2  (Pq/14. 7)°* 16  (Tto/i000) 


0.8 


(74) 


JP-4 


*LB0  =  °'22 


1 


U 


0.10 


1 


(h/12)0-2  (P0/14.7)0'15  (Tto/1000)0-8 


(76) 


These  expressions  apply  in  the  following  range  of  conditions: 

a.  Inlet  Pressure: 

b.  Inlet  Temperature:  1000  <  T. 


10  <  p0  <  40  psia 


to 


<  1 250  R 


c.  Mach  Number:  0.2  <  M  <0.9 

o 

5.  For  JP-4  fuel,  the  existence  of  complete  stability  loops  for 
a  coaxial  step  was  demonstrated  for  the  first  time  by  operating  at  very 
low  combustor  pressures.  Rich  extinctions  were  not  studied  with  ethylene 
due  to  the  hazardous  nature  of  operating  at  these  conditions.  Both  the 
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lean  and  rich  limit  data  for  JP-4  fuel  were  collapsed,  using  a 
correlating  group  of  the  form 


6.  Neither  the  ethylene  nor  JP-4  data  were  collapsed  by  the 
correlating  parameter  {X)  previously  employed  in  dump  combustor  studies. 

7.  The  effect  of  vitiation  is  to  reduce  the  stable  operating 
range  of  a  given  combustor. 

8.  A  variation  in  combustor  length  to  diameter  ratio  from  3 
to  6  produced  little  variation  in  flame  stability  characteristics. 

9.  The  split  fuel-flow  technique  to  determine  the  air  mass 
flow  fraction  consumed  in  the  recirculation  zone  can  be  successfully 
applied  to  coaxial  dump  combustors.  For  small  diameter  ratios  (d/D), 
the  result  obtained  was  compatible  with  estimates  drawn  from  industrial 
furnace  theory. 

10.  For  the  conditions  investigated  in  this  study,  it  was 
observed  that  direct  injection  of  fuel  into  the  recirculation  zone 
resulted  in  a  reduction  in  combustor  noise  and  much  smoother  operation. 

11.  It  was  demonstrated  that  a  transparent  quartz  combustor 
can  be  successfully  employed  to  observe  the  combustion  process  under 
realistic  operating  conditions.  It  was  shown  that  combustion  is,  in 
fact,  an  oscillatory  process  with  the  reaction  zone  moving  to  and  fro 
along  the  combustor;  such  combustion  appeared  quite  steady  to  the 
naked  eye. 

12.  A  comprehensive  review  of  the  literature  concerning  the 
flame  stability  and  aerodynamics  of  coaxial  dump  combustors  has  been 

presented. 
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VII  RECOMMENDATIONS  FOR  FUTURE  WORK 


The  experimental  investigation  of  flame  stability,  together  with 
the  state-of-the-art  reviews  undertaken  in  this  study,  have  uncovered 
many  important  research  areas  relevant  to  the  development  of  both 
coaxial  ana  multi-inlet  dump  burners.  These  research  areas  are  dealt 
with  in  this  section  under  the  following  headings:  flame  stability, 
combustion,  and  aerodynamics. 

Flame  Stability  of  Coaxial  Dump  Burners 

The  current  investigation  was  limited  to  one  size  of  combustion 
chamber  (D  =  6").  However,  the  correlating  group  for  lean  extinction 
limits  would  be  expected  to  depend  on  combustor  size  as  shown  by  Eq(63). 
Thus  it  is  important  that  the  effect  of  varying  the  combustor  diameter 
(D)  be  explored.  This  is  an  essential  extension  of  the  present  test 
series,  both  to  demonstrate  and  expand  the  application  of  the  correlat¬ 
ing  groups  given  in  Eq(74)  and  Eq ( 76) .  Another  recommended  investigation 
would  be  to  utilize  a  fluid-cooled  combustor  which  would  enable  the 
effect  of  wall  temperature  on  flame  stability  to  be  systematically 
explored.  It  should  also  be  noted  that  the  use  of  ablative  liners  on 
practical  combustors  will  create  a  hot  wall  condition.  If  wall  tempera¬ 
tures  are  maintained  at  high  levels,  then  significant  modification  of 
the  flame  stabilizing  mechanism  will  ensue.  Thus  both  hot  and  cold 
wall  effects  need  to  be  studied.  The  basic  stirred  reactor  analysis 
may  readily  be  extended  to  non-adiabatic  boundary  conditions. 

More  basic  research  investigations  of  the  flame  stabilizing  process 
associated  with  both  sudden  expansion  and  step  configurations  are  needed. 
No  data  are  apparently  available  on  the  structure  of  the  recirculation 


zone  or  of  the  associated  shear  layer.  A  first  task  is  to  undertake 
detailed  flow  visualization  of  a  flame  stabilized  on  a  two-dimensional 
step  so  as  to  observe  the  mechanism  of  extinction.  The  current  investi¬ 
gation  has  also  shown  that  a  transparent  quartz  combustor  may  be 
successfully  utilized  under  "smooth"  burning  conditions.  Investigation 
of  such  a  combustor  with  non-intrusive  instrumentation  would  yield  data 
on  the  flow  field,  including  the  geometry  of  the  recirculation  zone, 
and  on  the  flame  stabilizing  mechanism  of  the  sudden  expansion 
configuration. 

Combustion  in  Coaxial  Dump  Combustors 

Two  effects  which  were  observed  in  this  current  investigation  are 
important  to  combustor  development.  The  first  is  the  basically  unsteady 
nature  of  a  combustion  process  associated  with  the  fluctuating  reattach¬ 
ment  of  the  flow  to  a  solid  surface.  This  effect  was  observed  in  the 
movies  made  of  the  combustor  under  burning  conditions.  This  basically 
unsteady  flow  appears  to  couple  with  various  acoustic  and  instability 
mechanisms  producing  sudden  excursions  of  efficiency,  heat  transfer, 
and  noise  level.  It  is  therefore  important  in  all  tests  of  dump  com¬ 
bustors  that  appropriate  quick  response  instrumentation  be  employed  to 
characterize  the  unsteady  modes  of  operation,  and  thereby  relate  the 
changes  in  performance  to  the  driving  phenomena.  A  second  observation 
concerns  the  noticeable  reduction  in  combustion  noise  level  which 
occurred  when  using  direct  injection  into  the  recirculation  zone. 

This  injection  technique  is  not  generally  employed  in  current  dump  com¬ 
bustors.  However,  it  offers  the  immediate  advantage  of  direct  control 
of  the  local  fuel-air  ratio  in  the  recirculation  zone  together,  appar¬ 
ently,  with  the  possibility  of  suppressing  unsteady  modes  of  operation. 
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Aerodynamics  of  Coaxial  Dump  Burners 

As  will  be  noted  in  Appendix  A,  much  activity  has  been  focused  on 
the  one-dimensional  analysis  of  sudden-expansion  flows.  However, 
little  has  been  done  to  explore  the  structure  of  the  recirculation  zone. 
Thus  there  are  little  data  even  on  the  basic  size  and  shape  of  the 
recirculation  zone  as  a  function  of  Reynolds  number,  Mach  number,  and 
the  geometric  variables  such  as  the  ratios  Lc/D,  d/D,  and  d*/D  (where  d* 
is  the  exit  nozzle  throat  diameter).  In  particular,  existing  flow 

3 

studies  do  not  appear  to  embrace  the  range  of  Reynolds  number  from  10 
5 

to  10  .  Again,  few  investigations  have  examined  the  unsteady  aspects 
of  these  flows;  most  data  refer  to  time-mean  quantities.  Also,  the 
effects  of  initial  boundary  layer  thickness,  inlet  flow  profile,  and 
initial  turbulence  characteristics  on  the  flow  structure  downstream  of 
the  dump  plane  require  investigation.  Finally,  the  effect  of  combustion 
on  the  aerodynamics  of  the  flow  is  an  area  which  appears  to  have  received 
little  attention;  in  this  case  the  detailed  study  of  the  aerodynamics  of 
reacting  flows  associated  with  a  two-dimensional  step  would  be  a  good 
starting  point.  The  associated  heat  and  mass  transfer  phenomenon  are, 
of  course,  of  considerable  relevance  to  combustor  design.  Finally,  the 
use  of  residence-time  distribution  (RTD)  techniques  to  observe  the  gross 
mixing  characteristics  of  dump  combustors  should  be  explored.  A  few 
investigations  of  the  coaxial  configuration  have  been  made  and  are 
referenced  in  this  document.  However,  much  work  is  needed  to  charac¬ 
terize  the  coaxial  geometries  relevant  to  ramjet  combustors.  Such 
techniques  should  be  explored  to  permit  a  relatively  simple  approach 
to  the  modelling  of  the  complex  flow  structures  observed  in  multi¬ 
inlet  combustors. 
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Summary 


It  is  clear  from  the  above  discussion  that  the  coaxial  combustor 
offers  many  opportunities  for  research  relevant  to  integral-rocket- 
ramjet  combustors.  Such  problems  may  conveniently  be  addressed  with 
relatively  simple  test  geometries  but  will  require  sophisticated 
instrumentation.  This  area  of  endeavor  should  prove  a  fertile  field 
of  work.  However,  such  work  needs  to  be  regarded  as  a  stepping  stone 
to  the  understanding  and  modelling  of  the  complex  combustion  processes 
associated  with  multi-inlet  combustor  systems  which  are  currently  under 
development. 
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APPENDIX  A 


AERODYNAMICS  OF  THE  COAXIAL  DUMP  COMBUSTOR 

It  will  be  apparent  from  the  discussion  given  in  Chapter  II  that  in 
estimating  flame  stability  limits,  both  the  chemical  and  aerodynamic 
aspects  of  the  flameholder  flow  field  must  be  modelled.  It  is  there¬ 
fore  appropriate  to  review  the  aerodynamic  phenomena  associated  with 
the  dump  combustor. 

A  number  of  studies  deal  with  the  one-dimensional  analysis  of  the 
overall  stagnation  pressure  losses  incurred  by  sudden-expansion  systems 
(Refs  96  and  97).  Such  analyses  are  based  on  the  equations  of  conser¬ 
vation  of  mass  and  momentum.  A  key  assumption  is  that  for  subsonic 
flow  conditions,  the  pressure  at  the  plane  of  the  dump  face  is  uniform 
and  equal  to  the  corresponding  pressure  of  the  dump  plane  inlet  flow. 

Good  agreement  with  cold  flow  experimental  pressure  losses  is  obtained, 
providing  that  the  effects  of  inlet  velocity  profile  and  turbulence  are 
included  (Refs  98,  99,  and  100).  However,  for  the  purposes  of  this 
investigation,  the  specific  details  of  the  flow  field  are  of  importance 
rather  than  the  overall  description  provided  by  a  one-dimensional 
analysis. 

The  axisymmetric  sudden-expansion  configuration  is  a  ubiquitous 
component  of  various  flow  systems.  Analyses  of  the  flow  characteris¬ 
tics  of  this  device  appear  in  the  literature  of  fluid  mechanics,  heat 
transfer,  industrial  furnace  practice,  chemical  engineering,  fluidics, 
and  even  biological  engineering.  However,  much  of  this  literature 
deals  with  laminar  flow  conditions  and  with  geometries  of  low  diameter 
ratio,  i.e.,  d/D  less  than  about  0.1.  Unfortunately,  for  ramjet  systems, 
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data  are  required  for  highly  turbulent  flows  and  for  high  diameter 
ratios,  i.e.,  d/D  usually  greater  than  0.5. 


The  key  feature  of  the  sudden  expansion  system  is  the  existence 

of  the  recirculation  zone  generated  in  conjunction  with  the  separation- 

reattachment  process.  As  shown  in  Fig  A-l,  the  entering  flow  separates 

from  the  wall  and  entrains  material  from  the  surroundings,  thereby 

increasing  the  size  of  the  jet  flow  and  progressi vely  reducing  jet 

velocity.  The  jet  flow  rate  ultimately  reaches  a  maximum  value.  After 

this  point,  a  region  of  fluid  loss  from  the  jet  to  the  recirculation 

region  occurs  so  that  subsequently  the  flow  rate  of  the  reattached 

flow  becomes  equal  once  again  to  the  initial  flow  rate.  The  variation 

★ 

of  the  recirculated  mass  flow  fraction  (^rc/^0)  with  distance  is  shown 
schematically  in  Fig  A-2.  The  process  of  entrainment  and  fluid  loss 
(de-entrainment)  will  be  considered  in  more  detail  later  in  this 
Appendix. 

The  flow  field  in  a  coaxial  dump  combustor  may  be  sub-divided  into 
three  regions:  the  recirculation  zone,  the  central  non-recirculating 
flow  field,  and  the  shear  layer  which  separates  the  first  two  regions. 

Obviously,  the  volume  of  the  recirculation  zone  and  the  mass  flow 
entrained  into  it  are  of  interest  if  this  zone  is  modelled  as  a  stirred 
reactor.  The  details  of  the  shear  layer  would  also  be  of  interest  if  a 
characteristic-time  model  of  the  flame  stabilizing  process  were  used. 

The  central  flow  region  in  which  flame  is  propagated  is  of  interest  in 
the  estimation  of  overall  combustion  efficiency.  The  structure  of  the 
recirculation  zone  will  now  be  discussed  in  more  detail. 

*  mrc  is  used  to  denote  the  measured  recirculated  mass  flow. 

1 
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Figure  A-l.  Flow  in  a  Sudden  Expansion 
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Characteristics  of  the  Recirculation  Zone 

Turning  first  to  the  area  of  fluid  mechanics,  the  literature  on 
axi symmetric  sudden  expansion  flows  is  large.  A  number  of  significant 
theoretical  and  experimental  investigations  is  listed  in  Table  A-I. 
Attention  is  drawn  first  to  the  pioneering  work  of  Kalinske  (Ref  101), 
who  studied  the  mean  and  fluctuating  water  velocities  in  a  series  of 
divergent  pipes  including  a  sudden  expansion  configuration.  This  inves¬ 
tigation  utilized  a  streak  line  photographic  approach,  and  thus  the 
reduction  of  data  was  a  laborious  undertaking.  An  investigation  of 
similar  pipe  geometries  to  those  of  Kalinske  was  carried  out  by 
Chaturvedi  (Ref  89),  using  air  as  a  working  fluid.  In  addition  to  the 
mean  flow  velocity  profile,  the  turbulence  quantities  u',  v‘ ,  and  w' 
were  measured,  using  hot-wire  instrumentation.  Chaturvedi's  data  are 
limited  to  only  one  diameter  ratio  (d/D  =  0.5),  but  still  represent  the 
most  comprehensive  data  set  available,  and  indeed  have  served  as  a 
standard  against  which  theoretical  analyses  (Refs  87  and  102)  can  be 
assessed.  Some  descriptive  features  of  the  flow  field  (Ref  89)  are 
shown  in  Fig  A-3.  The  reverse  flow  velocity  in  the  recirculation  zone 
is  clearly  seen.  The  variation  of  the  axial  turbulence  intensity  with 
distance  is  also  of  considerable  interest  in  that  the  turbulence  is 
observed  to  diffuse  relatively  slowly,  and  the  high  intensity  level 
remains  closely  confined  to  the  shear  layer  for  several  diameters 
downstream. 

As  noted  earlier,  a  quantity  of  some  importance  for  combustor 
modelling  is  the  length  from  the  plane  of  the  dump  to  the  point  where 
the  flow  attaches  to  the  wall.  From  a  study  of  water  flow  in  a  sudden 
expansion.  Back  and  Roschke  (Ref  103)  demonstrated  that  the  reattachment 


TABLE  A- I  INVESTIGATIONS  OF  AXISYMMETRIC  SUDDEN  EXPANSIONS 


Rouse,  Experimental  Water 

Jezd insky 


Data  From  Ref  89 


distance  L,  expressed  in  terms  of  step  height  h,  varied  strongly 
with  the  upstream  Reynolds  number  (R  )  in  the  range  20  <  R£  <  4200. 
This  variation  is  illustrated  in  Fig  A-4.  Many  other  investigators 
(Refs  33-86  and  89)  have  determined  this  reattachment  length,  and  their 
data  are  also  shown  in  Fig  A-4.  It  will  be  observed  that  in  the 
laminar  flow  regime,  the  reattachment  length  L  increases  with  R  . 
However,  a  regime  is  reached  where  L/h  attains  a  maximum;  further 
increases  in  Rg  result  in  a  decrease  in  L/h  to  a  minimum  followed 
by  a  final  rise  to  an  approximately  constant  value  at  a  higher  Reynolds 
number  (Re  >  3.10^).  This  constant  value  of  L/h  is  denoted  by  k, 
so  that  k  is  approximately  equal  to  8.5.  Two  recent  investigations 
of  sudden  expansion  burners  are  reported  by  Pennucci  (Ref  84)  and 
Drewry  (Ref  85),  respectively.  Both  investigators  used  surface  flow 
visualization  to  locate  the  reattachment  point;  of  course,  the  reattach¬ 
ment  process  really  occurs  over  a  zone,  and  the  choice  of  a  distinct 
line  of  reattachment  cannot  be  precise. 

In  a  recent  application  of  the  laser-Doppler  velocimeter,  Moon  and 

Rudinger  (Ref  106)  measured  the  value  of  k,  using  a  coaxial  lucite 

5  6 

configuration  for  Reynolds  numbers  in  the  region  of  10  to  10  .  Their 
results  are  also  shown  in  Fig  A-4  and  are  generally  consistent  with 
those  of  other  investigators. 

It  may  be  inferred  from  Fig  A-4  that,  for  high  values  of  Reynolds 
number  (Re  >  3.105),  the  value  of  k  lies  between  8  and  9  for  cold 
flow  and  is  independent  of  Re.  However,  there  are  indications  in  the 
literature  that  two  geometric  variables,  namely  the  inlet  diameter 
ratio  d/D  and  the  exit  diameter  ratio  d*/D,  may  influence  k.  Regarding 
the  latter  ratio  (d*/D),  both  Pennucci  and  Drewry  used  geometries  which 
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Reynolds  Number 

Figure  A-4.  Recirculation  Zone  Length-Ratio  as  a  Function  of  Reynolds  Number 


incorporated  a  terminal  convergent  nozzle,  thus  imposing  a  downstream 
boundary  condition.  Pennucci  observed  that  for  one  diameter  ratio 
(d/D  =  0.5),  decreasing  d*/D  from  0.63  to  0.45  led  to  a  small  decrease 
in  the  reattachment  distance  corresponding  to  a  reduction  in  k  from 
9  to  8.  A  similar  effect  was  noted  by  Drewry.  This  result  may  imply 
that  k  exhibits  a  Mach  number  dependency.  Another  intriguing  inves¬ 
tigation  is  that  of  Goldshtik  and  Silent'ev  (Ref  88),  who  studied  the 
effect  of  blockage  in  a  two-dimensional  opposed-step  configuration. 
Blockage  is  defined  as  the  ratio  of  flow  area  before  the  step  to  flow 
area  after  the  step.  In  this  two-dimensional  case,  the  corresponding 
value  of  k  varied  from  6  at  a  blockage  of  0.4  to  a  maximum  of  7  at  a 
blockage  of  0.68,  and  finally  fell  to  a  minimum  value  of  about  3  as  the 
blockage  approached  unity.  Although  this  investigation  was  performed  on 
a  two-dimensional  configuration,  similar  results  may  apply  to  axisymmetric 
systems.  Indeed,  Lipstein  (Ref  107)  indicates  such  an  effect.  Thus  the 
assumption  that  k  lies  between  8  and  9  step  heights  (Fig  A-4)  must  be 
treated  with  caution  pending  further  investigation  of  geometric 
variables. 

Both  Pennucci  and  Drewry  also  recorded  wall  static  pressure 
variation  in  the  chamber  and  measured  the  total  pressure  losses  asso¬ 
ciated  with  sudden  expansion  flows.  The  typical  variation  of  wall 
pressure  is  shown  in  Fig  A-5;  it  will  be  noted  that  pressure  continues 
to  rise  for  some  distance  after  the  reattachment  zone.  Neither 
Pennucci  nor  Drewry  was  able  to  correlate  the  location  of  the  flow 
reattachment  point  with  any  characteristic  feature  of  the  pressure  rise 
curve.  However,  a  re-examination  of  their  data  shows  that  the  reattach¬ 
ment  point  corresponds  closely  to  the  point  where  the  static  pressure 
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rise  reaches  half  its  maximum  value;  i.e.,  1/2  (Pmax  -  Pn)^ n) -  This  is 
entirely  an  empirical  observation.  It  is  not  intended  to  explore  this 
observation  any  further  in  this  discussion;  the  reader  is  referred  to 
other  investigations  (Refs  104  and  105)  which  discuss  correlations 
already  obtained  in  related  geometries.  Another  interesting  observa¬ 
tion  noted  by  both  Pennucci  and  Drewry  was  the  existence  of  a  rotating 
flow  in  the  RZ  region;  looking  upstream,  this  flow  rotated  in  a  counter¬ 
clockwise  direction.  This  phenomenon  was  also  reported  by  Iribarne 
et  al •  (Ref  86) . 

In  connection  with  the  reattachment  process,  it  is  important  to 
stress  that  this  reattachment  point  fluctuates  rapidly  about  the  loca¬ 
tion  indicated  by  flow  visualization  techniques.  An  interesting 
investigation  was  conducted  by  DeRossett  and  Przirembel  (Ref  83),  who 
studied  the  unsteady  nature  of  the  reattachment  process  behind  an 
axisymmetric  step.  They  utilized  a  dual-element  hot  film  probe  to 
measure  the  flow  velocity  close  to  the  wall  and  demonstrated  that  in 
the  reattachment  region,  flow  oscillations  existed  at  a  frequency  of 
the  order  of  several  hundred  Hertz.  Thus  reattachment  occurred  over 
a  region  rather  than  on  a  line;  furthermore,  the  width  of  this  region 
was  found  to  increase  with  both  step  height  and  the  thickness  of  the 
approaching  boundary  layer. 

Two  other  facets  of  this  investigation  are  worth  noting.  First, 
in  order  to  describe  the  reattachment  process,  a  "forward  flow  proba¬ 
bility  (FFP)"  was  defined.  At  a  given  location,  this  was  described  as 
the  percentage  of  time  that  the  flow  was  directed  upstream;  this  was 
observed  over  a  sufficiently  long  period  to  obtain  a  valid  statistic. 

A  typical  variation  of  this  FFP  with  distance  from  the  step  is  shown 
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in  Fig  A-6,  and  the  length  of  the  reattachment  zone  may  be  defined  in 
terms  of  the  level  of  FFP  (typically  at  the  zero  or  100%  level). 

Second,  the  shape  of  the  curve  confirms  the  existence  of  a  small  vortex- 
type  flow  in  the  corner  of  the  step  in  addition  to  the  main  recirculating 
flow.  The  existence  of  this  type  of  vortex  was  detected  in  nominally 
two-dimensional  flows  both  by  Abbott  and  Kline  (Ref  108)  and  other 
researchers.  In  the  basic  dissertation  by  DeRossett  (Ref  109),  infor¬ 
mation  is  also  given  on  the  fluctuating  pressure  levels  which  accompany 
the  reattachment  process.  The  maximum  RMS  values  of  the  fluctuation 
associated  with  various  flow  velocities  appear  to  be  located  just 
upstream  of  the  point  of  complete  reattachment. 

This  latter  result  confirms  earlier  work  by  Rouse  and  Jezdinsky 
(Ref  110)  and  that  of  Narayanan  and  Reynolds  (Ref  111)  concerning  two- 
dimensional  flows.  Attention  is  drawn  to  the  fluctuating  nature  of 
the  reattachment  process  because  high  speed  movies  of  combustion  in  a 
coaxial  dump  combustor,  made  during  this  current  investigation,  clearly 
revealed  the  very  unsteady  nature  of  the  process.  Furthermore,  this 
unsteadiness  is  of  considerable  technological  interest  in  relation  to 
practical  combustors  where  severe  heat  transfer  and  coupled  combustion 
instability  effects  may  result  from  this  basic  fluid  dynamic  phenomenon. 

The  above  discussion  has  yielded  some  insight  into  the  nature  of 
the  reattachment  process  and  has  highlighted  the  relatively  constant 
value  of  k  (8-9)  for  highly  turbulent  flow  under  non-burning  conditions. 
Currently,  no  detailed  information  can  be  found  on  the  flow  structure 
of  the  simple  axlsymmetrlc  sudden  expansion  under  combustion  conditions, 
although  it  is  anticipated  that  a  number  of  studies  will  eventually  be 
made  with  non-interference  instrumentation.  An  investigation  was  made 


173 


Forward 

Flow 

Probability 

(Percent) 


! 

I 


Figure  A-5.  Illustration  of  Forward  Flow  Probability 
(Ref  83) 
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by  Baev  and  Tret'yakov  (Ref  58),  under  combustion  conditions,  of  a 
two-dimensional  opposed-step  combustor.  In  this  case,  the  length  of 
the  recirculation  zone  was  not  affected  by  fuel-air  ratio  changes,  but 
»  increased  substantially  with  the  exit  velocity  from  the  combustion 

chamber.  No  exit  nozzle  appeared  to  be  used  in  this  combustor. 

^  A  similar  discussion  applies  to  subsonic  flow  over  two-dimensional 

steps.  Here  again  a  substantial  literature  exists  concerning  plane 
flows  (Refs  111-117)  on  the  mean  flow  structure,  the  turbulence  levels, 
and  on  the  unsteadiness  of  the  reattachment  process,  However,  only 
limited  observations  exist  on  the  flow  structure  associated  with  actual 
combustion.  Thus,  Kawamura  (Ref  76)  found  that  the  length  of  the 
recirculation  zone  was  about  half  that  of  the  corresponding  cold  flow 
(3.5  step  heights  compared  to  7)  and  exhibited  little  variation  with 
either  fuel/air  ratio  or  velocity. 

In  closing  this  discussion  of  the  length  of  the  recirculation  zone, 
it  may  be  concluded  that  there  are  many  factors  which  influence  this 
length,  and  much  research  needs  to  be  done,  particularly  under  actual 
combustion  conditions.  Nevertheless,  as  illustrated  in  Fig  A-4,  a 
reasonable  assumption  is  that  k  »  8,5  for  high  Reynolds  number  flows. 

«  This  conclusion  is  also  supported  by  studies  of  the  heat  and  mass 

transfer  in  sudden-expansion  systems.  However,  such  secondary  evidence 
,  will  not  be  introduced  in  this  discussion. 

>  The  volume  and  surface  area  of  the  recirculation  zone  are  also  of 

interest.  A  typical  recirculation  zone  shape  is  shown  in  Fig  A-7. 

This  shape  may  be  approximated,  as  shown,  by  an  equation  of  the  form 


175 


The  surface  area  (Sp)  corresponding  to  this  parabolic  approximation 
is; 
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The  corresponding  volume  is  given  by 
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A  characteristic  dimension  (hp)  of  the  recirculation  zone  may  be 
defined  by  writing 
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A  curve  of  hp/D  as  a  function  of  p-  for  various  values  of  k  is 
given  in  Fig  A-8.  For  k  =  8.5,  a  good  representation  of  this 
relationship  is 

it  1.095 
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In  addition  to  the  reattachment  distance  (L),  another  important 
quantity  is  the  mass  flow  entrained  into  the  recirculation  zone.  With 
few  exceptions,  information  on  this  quantity  is  largely  found  in  the 
literature  concerning  industrial  furnaces  and  burner  systems.  Rele¬ 
vant  aspects  of  this  literature  will  now  be  reviewed. 

Entrainment  in  Confined  Jet  Systems 

A  typical  arrangement  of  an  industrial  furnace  is  shown 
schematically  in  Fig  A-9.  Unlike  the  dump  combustor,  such  furnaces 
are  generally  fed  by  two  streams.  The  central  or  primary  jet  is  typi¬ 
cally  a  high  velocity  fuel,  or  fuel-laden,  jet.  The  secondary  flow 
which  surrounds  the  primary  provides  air  for  combustion.  Two  situa¬ 
tions  can  arise:  if  the  secondary  air  flow  is  large,  then  entrainment 
of  air  by  the  primary  jet  proceeds  steadily  until  restricted  by  the 
duct  wall.  However,  if  the  secondary  air  flow  is  below  that  which  the 
jet  can  entrain,  some  downstream  fluid  will  be  recirculated  upstream 
and  a  recirculating  flow  will  be  established.  Two  quantities  are  then 
of  prime  interest;  namely,  the  axial  location  of  the  recirculation  zone 
and  the  maximum  quantity  of  fluid  being  recirculated.  Two  theories 
have  been  widely  used  to  correlate  these  quantities:  the  early  work 
of  Thring-Newby  (Ref  118)  and  the  more  rigorous  analysis  of  Craya-Curtet 
(Ref  119).  These  two  theories  have  each  yielded  a  similarity  criterion 
describing  confined  jet  flows  commonly  denoted  by  the  symbols  6 
(Thring-Newby  criterion)  and  Ct  (Craya-Curtet  criterion).  An 
appreciable  volume  of  literature  (e.g..  Refs  120  and  121)  refining  and 
extending  these  theories  has  been  published  in  recent  years. 

The  relevance  of  the  above  work  to  the  present  investigation  lies 


in  the  fact  that  the  sudden  expansion  combustor  corresponds  to  an 


industrial  furnace  with  zero  secondary  air  flow;  this  special  case  is 
documented  in  the  literature.  Unfortunately,  the  area  ratios  (d/D)  of 
industrial  furnaces  are  much  lower  than  those  corresponding  to  ramjet 
combustors.  This  fact  limits  the  direct  application  of  available 
furnace  data  to  ramjets. 

Thring  and  Newby  (Ref  113)  analyzed  the  fow  within  the  furnace 
on  the  assumption  that  the  primary  jet  entrained  flow  in  a  manner 
similar  to  a  turbulent  free-jet.  Field  et  al.  (Ref  122)  modified  the 
analysis  utilizing  the  classical  expression  for  free-jet  entrainment 
derived  by  Ricou  and  Spalding  (Ref  123)  from  direct  flow  measurements. 


This  expression  is: 
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where 


m  =  mass  flow  rate  at  distance  x^  downstream 

m  =  initial  mass  flow  rate  of  jet 
o 

p,  =  ambient  density 

a 

Pq  =  jet  fluid  density 
d  =  jet  diameter 
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Eq(A-5)  is  valid  for  values  of  xL/d  greater  than  about  13.  For 
shorter  distances.  Hill  (Ref  124)  has  shown  that  lower  rates  of 
entrainment  exist. 

The  corresponding  entrained  mass  flow  mg  will  be  given  by: 
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The  basis  of  Field's  theory  is  that  emerging  jet  is  assumed  to 
entrain  as  a  free-jet  up  to  a  given  point  Q  where  entrainment  stops 
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(Fig  A-10);  disentrainment  then  takes  place  until  the  jet  strikes  the 
wall  at  point  R.  A  jet  spreading  angle  of  9.7  degrees  is  assumed  by 
Field  (which  incidentally  appears  to  be  twice  as  high  as  Field  ini¬ 
tially  intended),  leading  to 
xR  =  2.925D 


where  D  is  the  combustor  diameter. 


m_ 


A  point  P  is  defined  as  the  place  where  —  defined  by 

mg 

Eq(A-6)  falls  to  zero,  i.e., 
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Finally,  it  is  assumed  that  point  Q  is  located  midway  between 
points  P  and  R;  thus: 
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Finally,  from  Eq(A-6)  and  Eq(A-8),  and  equating  mrc  to  m( 
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where  Q  is  termed  the  Thring-Newby  parameter. 

In  some  cases  all  the  secondary  flow  (m.)  will  be  entrained  before 
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recirculation  starts.  A  modified  parameter  is  then  used,  defined  as: 
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The  corresponding  recirculating  flow  equation  is 
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Unfortunately,  the  above  theory  is  only  strictly  valid  for  values 
of  jy  <  0.05  whereas  in  the  present  investigation,  d/D  lies  in  the 
approximate  range  0.40  to  0.80.  Nevertheless,  it  is  necessary  to  dis¬ 
cuss  the  6  parameter  so  that  relevant  experimental  data  may  be 
presented.  More  exact  analyses  of  furnace-type  flows  have  led  to 
improved  similarity  parameters. 

The  first  rigorous  analysis  of  confined  jet  flow  was  reported 

initially  by  Craya  and  Curtet  (Ref  119),  and  subsequently  by  Curtet 
(Ref  125).  This  work  considered  the  equations  of  motions  for  iso¬ 
thermal,  incompressible  flow  and  led  to  a  similarity  parameter  denoted 


by  ms,  which  can  be  expressed  as: 
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where  is  termed  the  "kinematic  mean  velocity"  and  is  expressed  as: 
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d  denotes  primary  jet  velocity  and  UQ  corresponds  to 


where  U 
secondary  flow  velocity 


The  term  Udyn  is  called  the  dynamic  mean  velocity  and  is  defined  as: 
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In  yet  another  analysis  of  the  general  problem  of  confined  jet  systems, 
Becker  et  al.  (Ref  126)  defined  a  parameter  which  they  termed  the  Craya- 
Curtet  number  C^;  this  is  related  to  the  m$  parameter  given  in 
Eq (A-l 3)  by  the  expression: 
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Corresponding  comprehensive  experimental  studies  have  shown  Ct  is  in 
general  an  appropriate  parameter  to  characterize  the  mixing  of  isothermal 
confined  jets.  For  large  duct  to  jet  area  ratios,  it  may  be  shown  that 
Ct  «  0.  Another  similarity  parameter  is  derived  by  Magnussen  (Ref  127), 
whose  theory  for  recirculation  was  based  on  the  pressure  increase  which 
occurs  downstream  of  the  jet  region.  His  similarity  parameter  is 
denoted  by  B,  where  B  is  defined  as: 

U+  (pA)c 


B  = 


(A-17) 


m 

where  U+  is  a  characteristic  entrance  velocity  and  the  subscript  c 
denotes  conditions  at  the  plane  of  maximum  recirculation. 

For  a  constant  area  duct  with  constant  density  flow,  it  can  be  shown 
that  B  is  related  to  the  Craya-Curtet  parameter  by: 
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Also  for  a  sudden  dump  configuration,  the  following  relation  between 
the  Thring-Newby  parameter  (0)  and  B  exists  (Ref  127). 

0  =  —  (A-19) 

»/eT 

There  are  thus  four  similarity  parameters  in  use  to  describe 
coaxial  confined  jet  mixing;  0,  Ct,  m$,  and  B.  Unfortunately, 
their  use  has  been  predicated  on  the  assumptions  of  isothermal, 
incompressible  flow,  and  on  large  duct-to-jet  areas.  As  noted  earlier, 
investigations  have  been  made  to  extend  the  analysis  of  jet  mixing;  in 
particular,  to  include  the  case  of  non-uniform  density,  and  a  non- 
isothermal  Craya-Curtet  number  has  been  defined  (Refs  120  and  121). 
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Also,  another  method  of  analyzing  cocurrent  flows  using  flow 
non-uniformity  coefficients  has  been  reported  (Ref  128). 

The  above  similarity  parameters  have  been  briefly  introduced 
because  some  data  from  furnace  practice  are  relevant  to  ramjet  combus¬ 
tors  despite  the  significant  differences  in  duct  area  ratios.  Further¬ 
more,  the  analysis  of  the  more  complicated  multi-inlet  ramjet  combus¬ 
tors  (Fig  2)  will  require  an  appeal  to  similarity  parameters  in  order 
to  establish  a  rational  approach  to  combustor  analysis.  The  similarity 
parameters  for  the  coaxial  dump  combustor  will  now  be  described. 

For  the  sudden  expansion  combustor,  there  is  no  secondary  air 
flow  so  that  Eq(A-14)  yields 

UK  ’  Ud  (jf)  2 

and 

Udyn  =  Ud  (f ) 

From  Eq(A-13) 
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ms  =  (J)  -  0.5  (A-20) 

and  from  Eq(A-16) 
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Also,  from  Eq(A-9) 

9  =  J  (A-22) 

and  finally  from  Eq(A-19) 
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The  values  of  M,  Ct, 


0  ,  and  B  for  the  combustor 


geometries  utilized  in  the  present  investigation  are  listed  in  Table 
A-II.  The  relevant  data  from  furnace  practice  will  now  be  reviewed. 

A  number  of  investigations  have  been  performed  to  measure  the 
recirculating  mass  flow  fractions  in  both  cold  and  combusting  flows. 

A  most  useful  source  of  data  is  the  work  of  Hubbard  (Ref  90)  because 
several  tests  were  performed  with  zero  secondary  air  mass  flow.  Thus 
Hubbard  gives  data  on  axial  velocity  decay,  recirculation  zone  velo¬ 
cities,  and  recirculated  mass  flow  fraction.  In  Fig  A- 11 ,  it  is 
shown  that  the  recirculated  mass  flow  fraction,  corresponding  to  zero 
secondary  flow,  falls  on  the  curve  which  relates  that  fraction  to  the 
0  parameter.  This  suggests  that  data  obtained  for  two  stream  systems 
can  be  utilized  for  sudden  dump  combustors.  Specifically,  it  suggests 
that  for  constant  density  mixing,  the  quantity 


can  be  regarded  as  the  equivalent  diameter  ratio  for  a  sudden  expansion 
system.  In  physical  terms,  this  requires  that  the  dump  plane  corre¬ 
spond  to  the  point  where  all  the  secondary  flow  is  entrained  prior  to 
the  commencement  of  recirculation.  Another  unique  feature  of  Hubbard's 
data  is  that  they  include  relatively  large  values  of  d/D;  his  data  is 
plotted  in  Fig  A-12,  together  with  the  recirculation  fractions  rec¬ 
ommended  by  other  sources.  The  following  expressions,  suitably  modi¬ 
fied  to  apply  to  the  coaxial  dump  combustor,  are  recotimended  by  the 
various  sources: 
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Recirculated  Mass  Flow  Ratio 
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Figure  A-ll.  Variation  of  Recirculation  Mass  Flow  Fraction  withym$ 
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Figure  A-12.  Estimates  of  Recirculation  Mass  Flow  Fraction  from 
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Further  inspection  of  Hubbard's  data  reveals  that  several  of  the 
larger  d/D  points  correspond  to  sudden  dump  operation.  These  points 
have  been  flagged  in  Fig  A-12.  On  inspection,  it  appears  that  the 
limited  sudden  expansion  data  arising  from  furnace  studies  are  reasona¬ 
bly  well  correlated  by  Wingfield's  expression: 

sf  =  °-398  (370  -  K059)  <A-29> 

It  should  be  recalled  that  the  above  discussion  related  to  cold  flow 
models  and  that  appropriate  allowance  for  change  of  density  must  be 
made  before  applying  Eq(A-29)  to  hot  flow  situations.  In  this  regard, 
Beer  et  al.  (Ref  129)  have  demonstrated  similarity  of  hot  and  cold 
mixing  patterns  using  a  modified  Thring-Newby  parameter.  The  correc¬ 
tion  for  hot  flow  conditions  is  simply  to  use  an  equivalent  primary 
nozzle  diameter  (d^)  through  which  passes  the  same  mass  flow  and 
momentum  flux  as  the  actual  nozzle,  but  with  a  density  Pfr  instead 
of  P0.  Thus 
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A  rigorous  derivation  of  this  correction  term  is  given  by  Guruz  et  al. 
(Ref  121).  The  effect  of  density  ratio  on  the  recirculated  mass  flow 
fraction  given  by  Eq(A-29)  is  shown  in  Fig  A-13. 


In  this  review  of  furnace  data,  it  has  been  stressed  repeatedly 
that  the  results  apply  to  combustors  of  lower  diameter  ratio  (d/D ) 
than  ramjet  systems.  In  these  low  diameter  ratio  systems,  the  rate  of 
entrainment  may  be  derived  from  the  free-jet  expression  of  Ricou  and 
Spalding  (Ref  123)  given  in  Eq(A-5).  However,  during  the  course  of 
this  review,  it  was  noted  that  Hill  (Ref  124)  had  obtained  entrainment 
data  for  the  initial  region  of  a  free-jet.  This  data  is  more  appro¬ 
priate  to  coaxial  combustors  of  high  diameter  ratio  where  the  primary 
jet  is  not  fully  developed.  In  this  regard,  Dealy  (Refs  130  and  131) 
has  shown  that  similar  velocity  profiles  are  developed  at  ^  =  0.25, 
but  not  for  ^  =  0.50.  Therefore,  the  results  of  Hill  were  reviewed 
and  an  analytical  expression  derived  for  his  data. 

Entrainment  in  the  Initial  Region  of  a  Turbulent  Jet 

The  experimental  investigation  of  Hill  (Ref  124)  utilized  the 
same  technique  as  Ricou  and  Spalding,  but  obtained  data  for  entrainment 
at  x^/d  ratios  in  the  region  1  <  xL/d  <  21.  Hill's  data  is 
reproduced  in  Fig  A- 14  This  shows  the  variation  of  the  entrainment 
coefficient  C2  with  xL/d.  Note  that  the  entrainment  rate  is 


where  P  ,  is  the  ratio  of  secondary  stream  to  primary  stream 
re  1 

density  and  by  comparison  with  Eq(A-5),  it  is  seen  that  C2  =  0.32  in 
Ricou  and  Spalding's  work. 


The  data  given  in  Fig  A- 1 4  may  be  approximated  by  the  expression: 


C2  =  ~  (b,  +  Ln  L 


-  (b,  +  Ln  \  0  <  xL/d  <  12 


yielding 


“e^L1  .  ft  l. 

d 


xL  xL/d 

ft 2  dxL  =  V^rel  fC 2  d  ^ 


and  Eq ( A-32 )  finally  yields 


h  "  rel 


0.033  +  0.083  Ln  ^ 


for  entrainment  fraction,  where  the  range  of  applicability  is 
XL 

o  <  -r  <  12 


(A-32) 


(A-33) 


For  xL/d  >  12,  Ricou  and  Spalding's  (Ref  123)  expression  applies: 


”  °'32 


(A- 34) 


A  comparison  of  the  entrainment  fractions  given  by  Eqs(A-33)  and 
(A-36)  is  given  in  Fig  A- 1 5 . 

Residence  Time  Distribution  Diagnostics 

As  noted  in  Chapter  II,  the  residence  time  distribution  (RTD)  of  a 
combustor  will  yield  information  about  the  type  of  flow  pattern  and 
thus  has  potential  for  use  as  a  diagnostic  tool.  Two  types  of  RTD  are 
of  interest;  namely,  the  RTD  for  the  recirculation  zone  itself,  and  the 
RTD  for  the  whole  combustor.  The  RTD  for  the  recirculation  zone  itself 
can  yield  information  on  the  mean  residence  time  (tl  of  fluid  in  the 
zone  and,  given  an  effective  volume  (V),  the  associated  mass  flow  rate 
can  be  deduced.  Thus  the  mass  flow  is  given  by 
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Figure  A-15.  Comparison  of  Entrainment  Fractions  Given  by  Two  Theories 
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Many  studies  have  been  made  of  the  RTDs  of  the  recirculation  zone 
behind  bluff  bodies  such  as  discs  (Ref  35).  However,  the  RTD  of  a 
simple  sudden  expansion  RZ  does  not  appear  to  have  received  attention. 

The  overall  RTD  of  the  combustor  may  be  used  to  estimate  the 
associated  combustion  efficiency  as  demonstrated  by  Beer  and  Lee  (Ref 
129). 

A  few  investigations  of  the  RTDs  associated  with  coaxial  dump 
combustors  have  been  reported  in  the  literature  (Refs  131-138).  Once 
again,  these  have  been  generally  concerned  with  combustors  of  low 
diameter  ratio  (d/D),  and  little  data  of  immediate  application  have 
been  generated.  However,  the  work  of  Moeller  and  Dealy  (Ref  134)  did 
yield  an  interesting  correlation  between  mixing  effects  and  the 
Craya-Curtet  parameter  (Ct).  Thus  the  RTD  approach  has  considerable 
potential  and  may  prove  to  be  of  distinct  value  in  determining  the 
significant  overall  characteristics  of  the  complex  flow  fields  found 
in  multi -inlet  combustors. 

Summary 

The  purpose  of  this  appendix  was  to  review  aerodynamic  phenomena 
associated  with  the  coaxial  dump  combustor  with  particular  regard  to 
obtaining  data  on  the  recirculated  mass  flow  fraction.  In  accomplish¬ 
ing  this  review,  it  is  apparent  that  much  further  work  remains  to  be 
done  to  determine  adequate  similarity  parameters  for  combustors  of  high 
diameter-ratio  (d/D)  and  for  multi-inlet  combustors. 
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APPENDIX  B 


MASS  TRANSFER  IN  RECIRCULATION  ZONES 


In  order  to  evaluate  the  air  loading  of  a  homogeneous  reaction 
zone,  it  is  necessary  to  evaluate  the  term  mr/V,  where  mr  is  the 
air  mass  flow  participating  in  the  reaction  and  V  is  the  associated 
reactor  volume.  The  evaluation  of  these  terms  is  not  straightforward. 
For  bluff  bodies  installed  in  ducts,  a  number  of  cold  flow  studies  have 
been  made  to  establish  both  the  volume  of  the  RZ  and  the  associated 
mass  transfer  phenomena  (Refs  35,  49,  and  56).  However,  under  actual 
combustor  conditions,  significant  changes  in  both  the  V  and  mr 
terms  occur,  and  few  experimental  investigations  other  than  the  work  of 
Winterfeld  (Ref  56),  Bovina  (Ref  49),  and  Lefebvre  et  al.  (Ref  93)  have 
been  made.  Also,  it  should  be  noted  that  diagnostic  explorations  of 
combusting  recirculatory  flows  are  usually  carried  out  at  conditions 
well  removed  from  blow-off  in  order  to  provide  quasi-steady  conditions. 
Thus,  the  fluid  dynamics  of  the  actual  unsteady  blow-off  process  are 
not  documented.  Thus,  even  in  the  case  of  bluff  body  stabilizers, 
there  is  little  experimental  data:  the  sudden-expansion  configuration 
appears  to  have  received  no  attention. 

For  engineering  purposes,  it  is  desirable  to  express  m^  in 
terms  of  the  upstream  flow  condition.  Thus  for  a  bluff  body  stabilizer 


where  m0f  is  the  hypothetical  mass  flow  which  would  pass  through  an 


area  equal  to  the  flameholder  area  Af.  Winterfeld  (Ref  56)  has  called 
Qn  the  exchange  coefficient,  and  Ernst  (Ref  139)  terms  Cm  the  recircu¬ 
lation  efficiency.  In  general,  there  have  been  two  approaches  to  the 
determination  of  mr.  The  first  is  to  apply  a  mass  transfer  analysis 
(Ref  140).  Typically,  such  investigations  (Refs  21  and  32)  led  to  the 
result  that  Cm  was  approximately  constant,  and  thus  the  fraction  of 
the  approach  flow  which  entered  into  the  recirculation  process  was 
constant.  With  this  result  one  may  write 


mr  a  Po^Q^f 


(B-2) 


The  second  approach  to  estimating  the  mass  transfer  is  by  the  use  of 
residence  time  measurements.  This  approach  has  been  pioneered  by 
Quick  (Ref  141),  Winterfeld  (Ref  56),  and  Bovina  (Ref  49). 

For  a  recirculation  zone,  the  transported  mass  flow  mr  may  be 
expressed  as 

V 


Pr_ 

t 


(B-3) 


where  is  the  mean  density  of  the  recirculation  zone  and  t^  is 

the  mean  residence  time  of  a  particle  in  the  zone.  Winterfeld  has 
shown  that  for  a  number  of  bluff  body  shapes,  and  for  a  given  blockage, 
under  burning  conditions 

t  <*  1- 


U, 


And  from  Eq(B-3)  one  obtains 


"V  a  PrUoV 


(B-4) 


(B-5) 
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Winterfeld  also  defined  an  exchange  velocity  (Ey)  to  characterize  the 
turbulent  exchange  process  across  the  boundary  of  the  recirculation 
zone.  Thus 


mr  =  PrSEy  (B-6) 

The  exchange  velocity  was  found  to  vary  directly  with  the  freestream 
velocity  so  that 

"V  «  PrU0S  (B-7) 

For  a  given  shape  of  bluff  body  flameholder  installed  in  a  fixed 
diameter  pipe,  the  ratio  V/S  will  not  vary  significantly  during  opera 
tion.  Thus  Eq(B-5)  and  Eq(B-7)  are  quite  consistent. 

Shchetinkov  (Ref  45)  has  observed  from  some  earlier  data  of 
Winterfeld  (Ref  142)  that  the  dimensionless  volumetric  flow  rate 
through  the  RZ  is  independent  of  the  presence  of  combustion  and  varies 
little  with  blockage.  Thus 

"V 

- * —  =  constant 

UoDf  pr 
or 


mr  a  prUQS  (B-8) 

which  is  a  similar  relation  to  that  given  in  Eq(B-7). 

In  the  case  of  two-dimensional  bodies,  the  phenomenon  of  vortex 
shedding  occurs  and  introduces  another  variable,  i.e.,  the  Strouhal 
number,  affecting  residence  time  measurements  (Ref  143).  With  this 
cautionary  remark,  attention  is  drawn  to  the  work  of  Silant'yev 
(Ref  144),  who  investigated  transfer  phenomena  at  the  boundary  of  a 


recirculation  zone  generated  by  the  flow  over  two-dimensional  steps. 
In  this  investigation  the  mean  residence  time  was  found  to  vary  as 
f ol 1 ows : 


t  « 


h_ 

IL 


(B-9) 


where  h  =  step  height. 

The  final  result  quoted  is 


t  = 


20  h 

u 


(B-10) 


The  definition  of  the  free  stream  velocity  (UQ)  is  a  little  ambiguous 
in  this  investigation  (Ref  144),  but  Eq(B-lO)  is  considered  the  correct 
interpretation.  No  comment  appears  in  the  work  of  Silant'yev  on  the 
possible  impact  of  vortex  shedding  on  the  mean  residence  time. 

Returning  to  the  discussion  of  flow  around  axisyrmietric  bluff 
bodies,  a  significant  series  of  cold-flow  studies  has  recently  been 
reported  by  Humphries  and  Vincent  (Refs  35,  145-148).  They  have  shown 
that  the  dimensionless  group  H  s  UT/D^  may  be  expressed  as: 


H 


1 


*f  kJ_  ,  \ 
>’  Df’  u2  *  u*  HJ 


(B-ll) 


where  =  integral  length  scale  of  free  stream  turbulence 
=  energy  of  free  stream  turbulence 
Dg  =  molecular  diffusion  coefficient 
u  =  molecular  kinematic  viscosity 
=  Strouhal  number 


If  molecular  diffusion  and  vortex  shedding  are  neglected,  then 

"  ‘  {"*’  5?  u7} 
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(B-12) 


Under  "smooth"  free  stream  conditions,  H  is  found  to  be  constant 
irrespective  of  Rg.  Under  turbulent  free  stream  conditions,  the 
corresponding  expression  for  H  becomes 


where  A 


(A) 

is  a  free  stream  turbulence  parameter. 


(B-13) 


An  interesting  experimental  result  for  discs  is  that  H  decreases  from 
a  value  of  about  7.5  at  A  =0  to  H=6at  A  =.04:  a  corre¬ 


sponding  reduction  in  the  L/D  of  the  recirculation  zone  from  2.25  to 
about  1.5  was  noted.  The  work  of  Humphries  and  Vincent  has  placed  the 
study  of  cold  flow  residence  times  on  an  improved  theoretical  basis  and 
has  laid  a  firm  foundation  for  similar  studies  for  reacting  flows. 

Although  residence  times  have  been  measured  on  other  geometries, 
e.g.,  blunt  vehicle  bases  (Ref  149),  no  studies  of  the  sudden  expansion 
recirculation  zone  appear  to  have  been  made.  Indeed,  the  only  investi¬ 
gation  related  to  the  coaxial  combustor  is  that  performed  by  Silant'yev 
(Ref  144)  concerning  flow  over  a  step:  his  results  are  given  by 
Eq(B-10) . 

For  a  coaxial  combustor  with  reacting  flow,  it  may  be  postulated 
that  the  entrainment  process  is  similar  to  that  of  a  free  jet.  As 
noted  in  Appendix  A,  this  is  a  reasonable  assumption  for  combustor 
geometries  of  low  d/D  ratio.  In  this  case,  mass  entrainment  from  the 
recirculation  zone  will  take  place  in  the  initial  region  downstream  of 

t  I 

the  dump  plane  up  to  a  point  say  x^.  Subsequently,  downstream  of  x^ 
mass  will  be  de-entrained  from  the  jet  to  supply  the  recirculation  zone. 
The  corresponding  entrained  mass  flow  can  then  be  expressed  as 
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and  a  mean  density  term  appears  in  this  equation. 

To  recapitulate  at  this  point,  it  will  be  recalled  that  the  aim 
of  this  Appendix  was  to  evaluate  mr  in  terms  of  upstream  flow  condi¬ 
tions.  As  a  result  of  this  review,  three  methods  of  expressing  mr 
are  seen  to  be  available.  Firstly,  one  may  express  iry  from  Eq(B-l) 
as  follows: 
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Of  course,  if  the  entrainment  relation  given  in  Eq(B-14)  is  truly 
descriptive  of  the  coaxial  combustor,  then  the  error  in  assuming 


Pr  =  pQ  is  relatively  smaller  since  the  square  root  of  the  temperature 


ratio  Jj-  is  involved  rather  than  the  actual  temperature  ratio. 
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APPENDIX  C 


CHEMICAL  REACTION  MODEL 


In  order  to  model  the  recirculation  zone  of  the  coaxial  dump 
combustor  as  a  stirred  reactor,  two  tasks  must  be  completed.  The  first 
is  to  perform  the  necessary  calculations  to  determine  the  maximum 
reactor  loading  *  as  a  function  of  equivalence  ratio  and  inlet 
temperature.  However,  4>mx  is  a  complex  function  of  <t>  and  T,  and 
thus  the  second  task  is  to  express  0mgx  in  a  simplified  form  suitable 
for  engineering  use.  Typically,  in  this  study  is  expressed  over 

maX 

a  limited  range  of  <t>  and  T  in  the  form 


<b  a 

max 


Ts  <Pt 


(C-l) 


Thus  this  Appendix  outlines  the  numerical  calculations  performed  to 
determine  the  exponents  in  expressions  of  the  form  given  in  Eq(C-l). 
These  calculations  were  performed  for  both  ethylene  and  JP-4  fuels  over 
a  wide  range  of  assumed  activation  energies  (E)  and  reaction  orders  (n). 

It  was  shown  in  Chapter  II,  Eq(2),  that  the  loading  of  a  stirred 
reactor  may  be  expressed,  for  an  overall  reaction  order  of  two,  as 


r 

k. 


VP. 


(i  +  m)  n 

R2 


e-E/RTr 


v.  Tr 
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iVVf 


(C-2) 


For  a  reaction  of  order  n  the  corresponding  expression  becomes 


mr  _  (1  +  m)  fk  1  e  E/RTr\  J_  n/2Y  n/2 
Vp  n  Rn  T  (n-0.5)j  «t>  Ao  Af 

assuming  that  XQ  and  X^  share  equal  reaction  orders, 
for  XQ  and  Xf  from  Eq (4a)  and  Eq(4b),  one  obtains 


(C-3) 

Substituting 
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0 


n/2 


"V_  .  h.  u*, )  1  y/4)} _ 

VPrn  Rn  t"'0'5  60  |(m+l)(x  +  y/4)  +  ♦  +  €♦  (y/4-l)|n 

(C-4) 

If  attention  is  confined  to  a  given  hydrocarbon  fuel,  then  x  and  y 
are  determined.  An  activation  energy  E,  a  reaction  order  n,  and  a 
reaction  rate  constant  may  also  be  assigned  for  the  reaction. 

The  reaction  temperature  is  calculated  using  the  expression 

Tr  =  l  +  rjc  AT 

where  AT  is  the  adiabatic  temperature  rise  corresponding  to  the  given 
values  of  4>  and  T,  r]c  is  the  temperature  rise  combustion  effi¬ 
ciency,  and  it  is  assumed  that  r)c  -  e.  As  previously  illustrated  in 
Fig  6,  for  a  given  equivalence  ratio  and  initial  temperature,  0 
attains  a  local  peak  value,  4>mx,  at  a  given  value  of  £  denoted  by 
€opt'  As  P°1nted  out  Previously,  the  value  0max  is  of  significance 
in  that  it  determines  the  maximum  mass  flow  the  reactor  can  sustain 
before  extinction. 

For  the  purposes  of  this  investigation,  it  was  initially  desired 
to  express  0max  in  the  following  forms: 

*mx  «  eT/t  <c-5> 

and 

0mav  «  Ts  <t> 1  (C-6) 

max 

The  expression  given  in  Eq(C-5)  has  previously  been  used 
successfully  by  Ballal  and  Lefebvre  (Ref  31).  However,  Eq(C-6)  yields 
a  stability  correlating  group  conveniently  expressed  in  exponent  form. 
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Absolute  values  of  are  not  required  for  the  evaluation  of  the 

^max 

exponents  occurring  in  either  Eq(C-5)  or  Eq(C-6),  so  values  of 

'  1  1 
^max  =  ^^nax  were  aciecluate'  Thus  values  of  cJ>max  were  calculated 

both  for  ethylene  (CgH^)  and  JP-4  fuels  burning  with  air,  over  a  range 

of  inlet  temperatures  varying  from  600-1200°R  and  for  equivalence 

ratios  ranging  from  0.3  to  1.0.  These  calculations  were  made  for 

ethylene  for  values  of  E  ranging  from  30  to  54  KCal/Mole  and  for 

values  of  n  from  1.25  to  2.00,  A  reduced  number  of  calculations  were 

made  for  JP-4  using  the  formula  Cg  gH^g  g  given  in  Ref  151.  In  Fig  C-l , 

a  typical  variation  of  4>max  W1tfl  *  anc*  T  is  shown  for  ethylene 

with  E  =  42  KCal/Mole  and  n  =  2.  This  latter  figure  is  of  particular 

interest  in  that  it  shows  that  Eq(C-5),  expressed  in  logarithmic  form, 

is  a  suitable  form  for  approximating  the  variation  of  <t>max  with  T 

and  <t>.  As  noted  previously  in  Eq(9),  the  term  b  is,  rigorously,  a 

function  of  the  equivalence  ratio.  Now,  for  each  combination  of  E 

and  n,  a  data  set  of  the  form  shown  in  Fig  C-l  was  obtained.  For 

each  set  of  data,  expressions  of  the  form  given  in  Eq(C-5)  and  Eq(C-6) 

I 

were  used  to  approximate  $max  over  the  range  of  <t>  and  T  corre¬ 
sponding  to  lean  extinction  conditions.  The  ranges  used  were  <t> ;  0.4 

to  0.7,  and  T;  800-1200°R.  The  goodness  of  fit  generally  obtained 
is  illustrated  in  Fig  C-2  and  Fig  C-3.  A  multiple  linear  regression 
program  was  used  to  approximate  $>max  and  values  of  the  exponents  b, 
s,  and  t  were  obtained.  In  Fig  C-4,  the  variation  of  the  exponents 
s  and  t  with  the  assumed  activation  energy  and  reaction  rate  is 
shown  for  ethylene:  similar  results  are  obtained  for  JP-4  as  shown  in 
Fig  C-5. 
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The  exponents  s  and  t  appear  in  Eq(70)  of  Chapter  III.  These 


exponents,  together  with  the  reaction  order  n,  are  used  to  obtain 
the  quantities  a, 0,7, a  of  Eq(71),  which  are,  of  course,  exponents  of 
velocity,  pressure,  step  height,  and  temperature,  respectively.  A 
typical  variation  of  these  latter  exponents  with  E  and  n  is  shown 
in  Fig  C-6. 

A  few  values  of  the  exponents  a, 0,7,  and  a  were  also  calculated 
for  a  reactor  loading  expression  based  on  the  hypothesis  that  the  mass 
flow  mr  was  given  by  Eq(B-17).  In  this  case  the  air  loading  term  was 
given  by 


A  chemical  reaction  model  based  on  Eq(C-l)  was  utilized  for  <J>max  and 
the  corresponding  temperature  ratio  (T/Tr)  appearing  on  the  right-hand 
side  of  Eq(C-7)  was  also  evaluated.  The  values  obtained  for  a, 0,7 
and  6  did  not  correlate  the  experimental  data  and  this  line  of  inquiry 
was  not  further  pursued. 
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APPENDIX  D 


INLET  FLOW  DATA 

As  pointed  out  in  Chapter  IV,  one  of  the  experimental  goals  was  to 
obtain  an  approximately  constant  velocity  profile  at  the  dump  inlet 
plane.  Previous  investigators  had  used  straight  pipe  entries,  and  little 
information  on  the  associated  velocity  profiles  was  available.  For  this 
investigation,  a  series  of  convergent  inlet  sections  of  cubic  wall  pro¬ 
file  were  designed  and  fabricated.  Cold  flow  calibrations  of  three  dump 
sizes,  namely  2.5",  3.5",  and  4.5",  were  undertaken.  For  these  tests, 
the  convergent  section  was  directly  coupled  to  the  combustor,  and  vari¬ 
ous  exit  nozzles  were  used  to  vary  the  Mach  number  at  the  dump  plane. 
These  inlets  were  calibrated  over  a  range  of  dump  plane  Mach  numbers 
and  mass  flow  levels  which  bracketed  the  anticipated  operating  condi¬ 
tions.  All  calibrations  were  carried  out  at  a  total  temperature  of 
about  530°R.  At  each  condition,  the  dump  plane  Mach  number  profile  was 
measured,  utilizing  a  traversing  pitot  static  tube.  The  inlet  Reynolds 
number  based  on  dump  plane  diameter  was  calculated  from  the  expression 


The  experimental  Mach  number  profiles  typically  obtained  are  shown 
in  Figs  D-l ,  D-2,  and  D-3.  It  will  be  observed  that  profiles  of  sub¬ 
stantially  constant  Mach  number  were  obtained.  The  cold  flow  Reynolds 
numbers  ranged  from  4  x  105  to  1.5  x  106.  For  hot  flow  experimental 
conditions,  the  corresponding  Reynolds  number  was  always  in  excess  of 
1  x  105. 
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Figure  D-l.  Flow  Profile  for  2.5"  Inlet 


Figure  D-3.  Flow  Profile  for  4.5"  Inlet 
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